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PROTON SPIN-LATTICE RELAXATION IN 
POLYATOMIC GASES! 


M. Broom, M. Lipsicas,? AND B. H. MULLER? 


ABSTRACT 


The proton spin-lattice relaxation time 7, has been measured in gaseous 
samples of methane, ethylene, and ethane as a function of pressure at room 
temperature and also at 193°K for methane. In the pure gases 7; is propor- 
tional to density, p, at low densities indicating that intramolecular interactions 
couple the spin systems to the lattice, as is the case in hydrogen gas. T;/p at low 
densities gives information on the mean square angle through which the mole- 
cules are rotated per collision. Relaxation due to paramagnetic O2 is observed 
at higher densities when oxygen gas is added as an impurity. The relaxation 
probability per collision with an oxygen molecule is about 5 times larger for the 
ethylene—oxygen system than for the other two systems studied. This anomaly is 
discussed in terms of the theory of Oppenheim and Bloom. It is shown that a 
study of the temperature dependence of 7; due to O2 impurities provides a 
new way of obtaining detailed information on the Lennard-Jones parameters 
for the interaction between O2 and the solvent molecules. 


1. INTRODUCTION 

It has been shown in an earlier paper (Lipsicas and Bloom 1961) that 
measurements of the proton spin-lattice relaxation time, 7), in He gas are 
capable of giving unique information about the anisotropic part of the H.-H2 
interaction. This is so because, unlike monatomic spin systems such as He’, 
the principal mechanism enabling the nuclear spin system in H2 gas to exchange 
energy with the translational degrees of freedom of the gas, and hence establish 
thermal equilibrium at the temperature 7 of the gas, involves those inter- 
molecular interactions which cause transitions between the rotational sub- 
states of the H2 molecules. The experimental results reported here of 7; in 
CHy,, C2H,4, and C2H¢ gases as functions of temperature and pressure demon- 
strate that the main spin-lattice relaxation mechanism in polyatomic gases is 
basically of the same sort as in He. Just as in He, there is an intramolecular 
dipole-dipole interaction between the magnetic moments of the protons on 
each molecule. Collisions between molecules cause molecular reorientation 
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and a consequent modulation of the intramolecular interaction. If the intra- 
molecular interactions are all known, one can, in principle, relate 7, to those 
anisotropic intermolecular interactions which produce transitions between the 
rotational substates of the molecules. 

The mechanism for 7; in polyatomic molecules mentioned above was first 
introduced to interpret 7, in liquids (Bloembergen, Purcell, and Pound 1948). 
The time dependence of the intramolecular interactions was calculated in 
terms of a model of the molecule as a sphere undergoing a rotational random 
walk in a viscous medium and 7) was related to the viscosity of the liquid. 
More rigorous treatments of the rotational random walk model have recently 
been carried out by Furry (1957) and Favro (1960) for spherical and non- 
spherical molecules respectively. Hubbard (1958) has incorporated some of 
these recent results in a careful treatment of the polyatomic spin system using 
Redfield’s (1957) density matrix formulation of relaxation theory. Hubbard 
relates 7, to the rotational diffusion coefficient D which appears in Furry’s 
theory. The rotational random walk model is thus on a fairly firm theoretical 
foundation. It is difficult, however, within the scope of the model itself to 
answer questions regarding its validity for different classes of liquids. Also, 
even if it is valid, the relationship of D to intermolecular interactions is not 
at all obvious. 

In the experiments reported here 7; is studied in polyatomic gases both at 
low densities, where only two-body collisions are important, and as a function 
of density up to liquidlike densities. It may, therefore, be possible eventually 
to answer the questions formulated above by a detailed analysis of our results 
and further study of relaxation in gases. We make a start in this direction 
for the low density gases, where the mean time between successive steps in 
the rotational random walk is simply the mean time between collisions as 
calculated from the transport properties of the gas. We are able to show 
from our 7; values that although the rotational random walk model may be 
applicable in CH, gas it is inapplicable to C2H, or C2He gas, because our results 
can only be understood if each collision between pairs of these latter molecules 
causes a very large molecular reorientation. 

Before proceeding to a detailed discussion of our results it will be useful 
to outline some of those basic differences between H» and the polyatomic 
molecules studied here which are important from the point of view of nuclear 
spin relaxation theory. 

The intramolecular interactions which contribute to 7, in Hy» are of two 
types, a dipole-dipole interaction between the magnetic moments of the two 
protons and a spin rotational interaction between the proton magnetic 
moments and the rotational magnetic moment of the molecule. The orders of 
magnitude of these two interactions are comparable with one another in He, 
each of them producing a magnetic field at each proton of order 30 gauss 
(Kellogg et al. 1939). Very few direct measurements of the intramolecular 
interactions have been made in polyatomic molecules, but it is expected that 
the spin rotational interaction will be much smaller than for He, and should 
produce negligible relaxation effects as compared with the dipole-dipole inter- 
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action. Some recent measurements on C,H» (Anderson and Ramsey 1960), 
where the rotational contribution to the magnetic field at the protons is of 
the order of 1 gauss, seem to bear this out. 

Hydrogen has two modifications, ortho and para, corresponding to total 
nuclear spins J = 1 and 0 respectively. Resonance is only observed due to 
the ortho-H2 molecules. Only odd rotational states (J = 1, 3, .. .) are allowed 
for ortho-H»z molecules, and because of the small moment of inertia of the He 
molecule, the excited rotational states (J > 1) are not sufficiently populated to 
produce appreciable 7, contributions even at temperatures up to room 
temperature (Lipsicas and Bloom 1961; Needler and Opechowski 1961). 
Measurements of 7; in He may be related directly to mean lifetimes of the 
molecules in each of the three magnetic substates of the J = 1 state. 

The differences between He and the polyatomic molecules are illustrated 
by the case of CH4, which has three modifications, two of which have non- 
vanishing nuclear spin. They are meta (I = 2), ortho (I = 1), and para 
(I = 0), which have statistical weights 5:9:2 respectively. Since the rotational 
states allowed for each modification are determined by symmetry considera- 
tions, and since transitions between different modifications are highly for- 
bidden, one might expect, in general, to observe a relaxation curve described 
by the sum of two exponentials, one each from the ortho and meta molecules. 
In all the experimental results to be described, we always report, within 
experimental error, only a single exponential recovery of the magnetization 
towards equilibrium. This is explained by examination of the expressions for 
the relative populations (Q2, Q1) of the J = 2 and J = 1 modifications, re- 
spectively, at a temperature 7 (Maue 1937). 


Q. = 5(1+7 exp(—12a)+9 exp(—20a)+ ...], 


aaa 


9 exp(—2a)+15 exp(—6a) +42 exp(— 12a) 


+54 exp(—20a)+99 exp(—30a)+..., 


where a = h?/2IokT = 7.68/T, 

Io = moment of inertia of the CH, molecule. 
Thus, the rotational states allowed for meta- and ortho-CHy, are 

metas f= (Ona... 

Ortho? <= 4 25.4. Ges esc 

At room temperature the value of J which we obtain by equating the 

rotational energy of the J state, J/(J +1)h?/2Jo, to kT is approximately J = 6. 
Thus, very high values of J are appreciably populated at the temperatures 
where our experiments have been done. Examination of (1) shows that the 
ortho and meta modifications share the same higher rotational states. Although 
there will be two time constants the difference between them will be small at 
high temperatures and difficult to resolve. As we shall discuss later, an addi- 
tional reason for anticipating more than one exponential arises from Hubbard’s 
(1958) analvsis, but there, too, it turns out to be difficult to resolve the two 
exponentials experimentally. 
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Because the high rotational states are appreciably populated by the poly- 
atomic gas molecules, they may be pictured as approximately classical mole- 
cules as far as their angular momentum properties are concerned. Thus, the 
rotational random walk model often used to describe molecular reorientation 
due to collisions with other molecules has some chance of success, since the 
molecules can, in principle, undergo very small reorientations. 


2. EXPERIMENTAL METHODS 

All the 7 measurements were made using the pulse technique. The pro- 
cedure used for most of the measurements here was a slight modification of 
the usual pulse technique, convenient for the measurement of long relaxation 
times. A succession or train of pulses was applied in a time short compared 
with 7;. This train ensured that the magnetization at the end of the train 
was close to zero. At a time ¢ later a similar train of pulses was applied and the 
height A (¢) of the first echo due to this train was measured. 7; was deduced 
from the usual relation 


(2) A(t) = A()[l—exp(—t/7))] 


where (©) is the value of A (¢) for t > 7. Usually we waited a time 2107, 
to obtain A(). In some experiments ¢ was measured using a stop watch, 
while in others we used a system of decade counters. 

For details of the electronics, high pressure system and low temperature 
system, the reader is referred to the paper by Lipsicas and Bloom (1961). 

The samples of methane, ethane, and ethylene gas used here were Research 
Grade (299.6% pure) or Pure Grade (299% pure) obtained from Phillips 
Petroleum Co. The impurities listed by the supplier were other hydrocarbons 
whose influence on 7; is negligible in these quantities. 

In some of the experiments to be described later oxygen gas was added 
purposely as an impurity. This was done originally to ensure that the experi- 
mental results, especially in the low density region, really were characteristic 
of the pure gas and not affected by paramagnetic impurities. Later we found 
that the influence of paramagnetic impurities on 7, was of interest for its own 
sake. The actual measurements of impurity concentration were done very 
crudely, however, and it is difficult to assign experimental errors to this 
important quantity. In practice we simply filled our high pressure bomb to a 
certain measured pressure of O2. We then added the gas being studied to the 
highest available pressure, usually less than 100 atmospheres. The ratio of 
Os concentration to solvent concentration was then calculated from the 
measured pressures and the known compressibility factors for the gases. The 
pressure dependence of 7; was measured by releasing small quantities of the 


mixture. 
Only in the case of CH, did we obtain an independent estimate of the O2 


concentration. This was obtained by measuring 7, in the liquid where the 
relationship between 7, and O- concentration is accurately known (Sandhu, 
Lees, and Bloom 1960). 
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3. THEORY 


7, may be expressed in terms of Fourier transforms of the correlation 
functions of those interactions which enable energy to be exchanged between 
the nuclear spin system and the other degrees of freedom (lattice) of the 
material in which the spins are located. As we have pointed out in the Intro- 
duction, the nuclear spin-molecular rotational interaction may be neglected 
for polyatomic molecules, so that only correlation functions associated with 
dipole-dipole interactions need be considered. Then, for a system of identical 
nuclei of spin J and gyromagnetic ratio y, Bloembergen, Purcell, and Pound 
(1948) and Kubo and Tomita (1954) have derived the formula 


(3) 7 = SG? 1(1+1)[Ju(wo) + Je(2e0)], 
(4) Jx(w) = r (F(t) Fi(t-+1)) exp(iwr)dr, k= 1,2, 
(5a) Fy(t) = = 3, Y21(0i;, bij) les ’ 


32 . ne 
(5b) F3(t) = f= p ® V22(0ij, Gis) / ris ‘ 


where the external magnetic field Ho is oriented along the z-axis and the 
vector Tf, joining the 7th and jth nuclei has polar and azimuthal angles 4; 
and ¢;; respectively. Y2:(8@,¢) and Y22(0,o) are spherical harmonics. wo = yHo 
is the Larmor frequency of the nuclei. The bracket ( ) indicates an average 
over the equilibrium ensemble for the lattice. 

The sum over j includes all those nuclei on the same molecule as the ith 
nucleus as well as those on neighboring nuclei. Contributions to 7; due to 
interactions between nuclear spins on different molecules are inversely propor- 
tional to density and hence negligible in low density gases. At liquidlike 
densities their contributions may be comparable to those of the intramolecular 
interactions. The evaluation of the intermolecular dipole-dipole terms in 
equation (3) in terms of molecular properties of dense gases or liquids has 
been discussed elsewhere (Bloembergen, Purcell, and Pound 1948; Torrey 
1953; Oppenheim and Bloom 1959, 1961). Here, we shall only discuss intra- 
molecular interactions so that the sum over j in (5) includes only other nuclear 
spins on a single molecule. In the following paragraphs we shall first consider 
the case of rigid molecules and then attempt to allow for the effect of internal 


degrees of freedom. 


A. Rigid Molecules 
Usually correlation functions in (4) may be approximated by 


(6) (F.()Flt+7)) = (FOF: ())e77" 


° e ° ° . 3.2 
where 7,, the correlation time, usually satisfies the inequality wor, K 1 for 
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gases and liquids. If this assumption and equation (6) are satisfied then 7 
for rigid molecules can be written down using equation (3) and the facts that 
r;j is time independent and (| Yim|?) = 1/4x 


7 1 2 6 
(7) T= 2y*h I (I+1)t. » Ta. 

For a molecule containing N protons (J = 4) equidistant from one an- 
other so that 7) = ro 


, a ee 
(8) +” 3(N—1) ~* 3(N—1) y's. 

Hubbard (1958) has recalculated the relaxation of the nuclear spin magneti- 
zation for rigid molecules having three or four identical nuclei of spin } 
equidistant from one another using Redfield’s theory (1957) which takes into 
account off-diagonal terms in the density matrix for nuclear spins. He finds 
that the relaxation is governed by two exponentials in each case so that, 
strictly speaking, one cannot define a unique 7}. For the 4-spin case, CHa, 
the z-component of magnetization M(t) was found to have the following time 
dependence as it relaxed towards its equilibrium value Mo. 


—_ l -_ ay 1/2 A 
¢ eis oh FE oe a eee cai J 
(9) Mo—M(t)={Mo—M(0)} \ 184 exp] —345 (73+2(46) Tr, 


ce i e as a 
+( 184 mp] gg 2) Ff - 


The coefficients of the two exponential terms are in the ratio 1/50. It is thus 
extremely difficult to separate the two exponentials experimentally. Our 
signal/noise was not sufficiently great to resolve two such terms and we always 
observed only a single time constant. If we assume that one measures a single 
time constant given by the time constant of the exponential with the larger 
coefficient, then 


T, = 0.225 JT) for CHa. 


This differs by less than 1% from 7; = eT obtained from equation (8), 
putting V = 4. 

Similar arguments lead to 7; = 0.332 T> for a 3-spin system using Hubbard’s 
theory while the usual (diagonal density matrix) theory gives 7; = 37». 

Hubbard (1958) has pointed out that if the off-diagonal terms in the density 
matrix are neglected, the result for 7, depends on the representation used. 
Here we have found close agreement between his result and 7, calculated 
assuming a diagonal density matrix in a representation in which J; for each 
proton is diagonalized. Hubbard’s calculation becomes complicated for a 
4-spin system such as CoH, where the protons are not equidistant; thus we 
shall take the above result as justification for relating 7, to 7. for CeHa, 
and also for more complicated polyatomic molecules, using equations (3)—(7). 

For CH4, where 79 = 1.78 X10-* cm, we have 
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(10) ' Tyre = 1.24X10-" sec? (CH,) 


while for C2H,4, where the known structure of the C2H, molecule (American 
Institute of Physics Handbook 1957) gives 1.86 A, 2.42 A, and 3.05 A for the 
interproton distances in equation (7) 


(11) Tit, = 3.85 X10—"! sec? (C2H,). 


In treating C2H¢ as a rigid molecule, one overestimates 7, by less than 10% 
by considering it as a composite of two identical 3-spin systems, since the 
protons on different methyl groups are comparatively far from each other. 
Then, one obtains 


(12) Tr, = 1.86X10—'! sec? (C2He). 


Thus, within the scope of the above theory, measurements of 7; may be 
used to evaluate 7, as defined by equation (6). Furry (1957) has shown that 
if the molecule is assumed to undergoa rotational random walk, then r, = 1/6D, 
where D is the “rotational diffusion coefficient”. In the model of a molecule 
behaving as a classical sphere of radius 6 in a medium of viscosity 


(13) t= “SRT 


It is not expected that equation (13) would be applicable to dilute gases, 
but it will be interesting to apply it to the dense gas. For the dilute gas, one 
can to a reasonable approximation define the mean time between collisions 
rt from the measured transport properties of the gas. For a hard-sphere gas 
(Hirschfelder 1954) 


: cin Gay" 
(14) = = 4a'p\—— 


where a is the diameter of the molecule, m is the mass of the molecule, and p 
is the number of molecules per cm*. 

Since the anisotropic interactions which produce molecular reorientation 
are short range one can picture the orientation of a molecule as changing only 
as a result of collisions whose frequency is defined by the usual kinetic theory 
cross sections. One can then, quite precisely, define two opposite limiting 
cases. 

Case 1. Weak Coupling Approximation 

The mean square angle of rotation per collision (a?) is assumed small. Then, 
the theory of Furry (1957) would hold and with the additional assumption 
that 7 is the mean time between successive rotations in the rotational random 
walk, one obtains 


(15) = = (a") K1. 


ce 
Thus, to test whether the rotational random walk theory holds, for the 
dilute gas one only has to examine the order of magnitude of 7/7,. If it does 
hold +/7,. gives a measure of (a?). 
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Case 2. Strong Coupling Approximation 
Each collision leaves the molecule randomly oriented. In this case it is 
easily shown that rt = 7,¢. 


B. Non-rigid Molecules 

In a non-rigid molecule such as CeHg, there is the possibility of each of the 
two CH; groups undergoing rotations about the C—C axis of the molecule in 
addition to the rotation of the molecule as a whole. The activation energy 
for CH; rotation in C2H¢ is approximately 3000 cal/mole (Kemp and Pitzer 
1937) so that at sufficiently high temperatures this rotation would be taking 
place to an appreciable degree even between collisions. There is also the 
possibility of collisions occurring which cause an appreciable reorientation of 
a CH; group but not of the molecule as a whole. A complete theory of 7; in 
gases should, therefore, be formulated in terms of a complete description of 
the internal degrees of freedom of the molecules. Such a theory might follow 
the lines of the detailed treatment of the Hz molecule by Needler and Ope- 
chowski (1961), for example. 

In the absence of such a theory we can at least draw some conclusions 
about the limiting case that the motions associated with internal degrees of 
freedom are much faster than the reorientational motions of the molecule as 
a whole. If vi, is a characteristic time associated with the internal motions, 
Muller (1954) has pointed out that for the case Tin, <K7,~, and also wor, K 1, 
the Hamiltonian describing the intramolecular dipole-dipole interactions may 
be replaced by its average over the internal motions and 7; calculated in the 
usual way using this new effective Hamiltonian. One always obtains a longer 
7, than for the rigid molecule. For rapid rotation of a CH; group around its 
symmetry axis, it is found that 7; is multiplied by a factor 4 (Muller 1954) 
so that for ethane we find, using equation (12) 


(16) Tit. = 7.4X10-"! sec? (CoH¢ with rapid internal CH; rotations). 


Relaxation by Paramagnetic Impurities 

The presence of O» or other paramagnetic impurities often provides a potent 
mechanism for 7), since the magnetic moments of paramagnetic atoms or 
molecules are about 1000 times as large as those of nuclei. In polyatomic gases 
at sufficiently low densities, however, the effect of a small amount of para- 
magnetic impurity should be negligible. The reason for this is that, for relaxa- 
tion by paramagnetic impurities, 7; should be inversely proportional to 
density, while for the pure gas 7; « 1/7, « density. 

For the gas, the contribution of a paramagnetic impurity to 7; may be 
shown to be (Bloembergen et al. 1948; Oppenheim and Bloom 1961) 


1/2 
(17) nz Talia Ne pore) +m) | i F(y)dy gas. 


The corresponding formula for the dense gas or liquid is 
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1 _ 16m _yalp)Np (°F) 4, 
Ds 3 (D,+Du)a . ¥ = 


(18) = liquid, 
where 7p is the nuclear gyromagnetic ratio, (u;) is the mean square value of 
the magnetic moment operator for the paramagnetic molecule, N, is the 
number of paramagnetic molecules per cm’, M/, and m are the masses of the 
paramagnetic molecule and the hydrocarbon molecule respectively, D, and 
Dy are the diffusion constants of the paramagnetic molecule and the hydro- 
carbon molecule respectively, and a is a characteristic length for the molecular 
interaction such as the Lennard—Jones parameter. 


(19) Fly) = | fee eae) ax | 


g(x) is the radial distribution function of the paramagnetic molecule and 
hydrocarbon molecule system, where x is the separation of a pair of molecules 
expressed in units of the characteristic length a. 





(20) g(x) = exp | “ Vie) | 


where V(x) is the potential energy of a pair of molecules separated by a 
distance xa. For a liquid or dense gas, g(x) is more complicated (Oppenheilm 
and Bloom 1961). 

Evaluation of the integral in equation (18) using a realistic intermolecular 
potential such as the Lennard-Jones potential must be done numerically. 
Results of such calculations are not yet available. We give here the result for 
J, F(y)dy assuming a square well potential for the intermolecular potentia 
and using equation (20) for the radial distribution function of the dilute gas. 
For best agreement with the Lennard-Jones potential in fitting second virial 
coefficients of gases, the parameters of the square well are as follows (Hirsch- 
felder et al. 1954): 


Vie) =o So 
(21) V(x) = Vo = 0.56¢ L<*< Ls 
Vie) =0 6 > 18 


where ¢ is the Lennard—Jones energy parameter. 
For this case, equation (17) becomes 





‘a A _ 4 ynlup)Np fies ay (¥2) 
ws fy & @ 2k7 “\R7 
where 

(23) f(x) = e+0.31 (e7/?—1)?—0.72 e*/*(e?/27—1); 


f(x) is plotted in Fig. 1. 
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kT 
Fic. 1. Plot of f(Vo/RkT) versus Vo/RT. The function f(Vo/kT) is defined by equation (23). 


4. EXPERIMENTAL RESULTS AND DISCUSSION 


Methane (CH,) Gas 

Measurements of 7; in methane were made at two temperatures, 193° K 
and 293° K. The former temperature is only 2 Kelvin degrees above the 
critical temperature. The results for 7; are shown in Figs. 2 and 4 as functions 
of density. Figure 3 shows the density as a function of pressure for the lower 
temperature to indicate the effect of three-body and higher-order collisions 
on the equation of state for the gas. The density measurements were made 
using the amplitudes of the nuclear magnetic resonance signal as described 
elsewhere (Lipsicas, Bloom, and Muller 1961). 

The main qualitative features of our results may be summarized as follows: 

(a) 7, is proportional to the density p in the low density gas. This confirms 
the over-all picture we have presented earlier of the relaxation being due to 
intramolecular interactions which are modulated by those collisions between 
pairs of molecules which cause molecular reorientation. 

(6) T:/p for low density CH, gas is about 200 times as large as for H». This 
partially reflects the fact that 7i7,, which is about 1.24X10-" sec? for CH, 
according to equation (10), is about 35 times as large as the corresponding 
value for He (Bloembergen ef al. 1948). As has been discussed in the Intro- 
duction, however, 7, has a somewhat different significance for the two cases. 

(c) The addition of paramagnetic impurities in the form of O2 gas does not 
influence the very low density results at all but does influence the higher 
density measurements. Thus, in the samples containing a relatively large 
amount of Os, 7; goes through a maximum value when plotted as a function 
of density. It is likely that, at the highest densities in the measurements at 
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Fic. 2. Plot of 7; versus density for CH, gas at 193° K; closed circles are measurements 


for samples containing 0.1% Ox. 
Fic. 3. Plot of density versus pressure for CH, gas at 193° K. 
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Fic. 4. Plot of 7; versus density for CH, gas at 293° K. 
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193° K, paramagnetic impurities are playing a role even in our purest sample. 

The slopes of the straight lines drawn through the low density 7; values in 
Figs. 2 and 4 give: 

T,/p = 3X10-? sec/Amagat, T = 193° K, 

(24) 
T;/p = 2.05X10-? sec/Amagat, T = 203° K. 

Using equation (10) to obtain 7, from the results of equation (24) and 
evaluating 7 from equation (14) using a = 3.8X10-* cm, the Lennard-Jones 
parameter for CH, (Hirschfelder et al. 1954), we obtain 


t/t. = 0.27, T = 193° k, 


(25) ¥ 
t/t. = 0.15, T = 293° K. 


Thus the inequality (15) is not really satisfied, especially at the lower 
temperature, and it is not clear whether one can interpret the values of 7/7, 
in equations (25) as giving the mean square angle of rotation per collision 
(a?). However, the fact that r/7, at 193° K is approximately twice as large as 
7/T. at 293° K is consistent with what one would anticipate under the approxi- 
mation (a?) <1. For this case, the angle of rotation of a molecule due toa 
collision with another having relative velocity V is proportional to 1/V. Thus 
(a?) « ((1/V)?) « T-!; (293/193) is equal to 1.52 as compared with (a?)193° «/ 
(a?) 093° x = 1.8. It would be interesting, therefore, to generalize Furry’s random 
walk treatment of the molecular rotations to include the case (a?) not neces- 
sarily much smaller than unity. As we shall see below, the CoH, and C:He¢ 
gases most certainly do not satisfy the condition (a?) < 1. 

Measurements of rotational dispersion by means of ultrasonic absorption 
experiments (Kelly 1957) indicate that ‘the number of collisions required to 
excite the rotational mode in methane lies between 14 and 17’. The fact 
that fewer collisions are required to reorient the CH, molecules according to 
our measurements indicates, as one might expect from energy conservation, 
that collisions which reorient the molecules without causing a change in the 
J quantum number are (about 2 times) more frequent at room temperature 
than collisions which cause a change in J. 

Turning now to relaxation of CH, by O: impurities, the results for a sample 
containing 0.1% O2 impurity are given in Fig. 2. As mentioned earlier, the 
amount of O2 was determined by measuring 7) in the liquid where the relation- 
ship between 7; and .\, is known experimentally (Sandhu, Lees, and Bloom 
1960). At the higher densities where the impurity relaxation predominates, 
Tp is approximately constant and equal to 350 Amagat-sec. Since there are 
0.1% Oz molecules this gives 7).V, = 0.35 Amagat-sec for O2 in CH, gas. 

The theoretical value for 7,;V, may be found by evaluating equation (22). 
In fact, if (u;) was known accurately, this experiment could be used to evaluate 
the parameters a, Vo associated with the O.-CH, intermolecular interaction. 
If one takes (up) = y(O2)*h*?S(S+1) = 0.67X10-* e.m.u. (American Institute 
of Physics Handbook 1957) and using y, = 2.67104, T = 193° K, 7,N, = 
0.35 Amagat-sec and the known masses of O2 and CH, molecules, equation 
(22) gives 
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f(Vo/kT) _ 


A 0.24. 
a(A)? 


(26) 
The Lennard-Jones force constants for O2 and CH, from second virial 
coefficients are (Hirschfelder ef al. 1954) 


° 


e/k = 148.2°K, a=3.817A, CH, 
e/k = 118° K, a = 3.46A, Oo. 


We approximate the force constants of the O.-CH, interaction using the 
empirical combining laws. 


(27) €12 = (€1€2)'/%, G32 = 3(a:+4a2), 
which gives 
(28) e/k = 182°K, a=3.64A, CH,-On. 


Using equations (21) and (28) we get Vo/kT = 0.38 at 193° K corresponding 
to f(x) = 1.29. Thus, we estimate 
f(Vo/kT) _ 1.29 


——~; = 0.098 


(29) a(A)? ~ @.64): 


in good agreement with the result obtained from 7 measurements and given 
in equation (26) considering the crudeness of the square well model and the 
use of the empirical combining laws. 

It will be interesting to see whether the agreement is improved after the 
integration in equation (17) has been carried out using a more realistic 
potential. Also of interest will be a comparison of the theory with 7; measure- 
ments carried out over a wide range of temperature. 


Ethylene (C2H.) Gas 
Measurements of 7; in ethylene gas at 293° K are shown in Fig. 5 for the 
purest sample available and for a sample with 0.35% O2 added. The main 
features of the results are that 7; is again proportional to density in the low 
density region with a slope about 10 times greater than for CH,, and that 
oxygen seems to be much more effective in shortening 7; than for CH, gas. 
From the linear region in the 7; versus p curve for the pure gas we obtain 


(30) T;/p = 0.25 sec/Amagat, T = 298° K. 


We follow the same procedure as used to interpret the methane results, using 
equation (11) to obtain 7, from the 7; measurements and equation (14) to 
obtain 7. In equation (14) we take the hard-sphere diameter of the C,H, 
molecules to be 4.2 A. We then obtain 

(31) - = 0.66, 


¢ 


3 | 


which indicates that the molecules are rotated through large angles in each 
collision. 
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Fic. 5. Plot of 7; versus density for C,H, gas at 293° K. 


When 1/7; is plotted against p for the sample containing 0.35% Oz it 
gives a straight line with p7,; = 27 Amagat-sec. Thus 7,V, = 0.096 Amagat- 
sec, which is about 1/4 the value obtained for methane at 193° K. Following 
the same procedure described in detail in the discussion of methane, we obtain 


F(Vo/kT) _ og. 


i a(A)? 


If we use the empirical combining laws (27) to evaluate the constants for 
the C2H,.-Oz interaction, we obtain from the Lennard-Jones parameters for 
CoH, (Hirschfelder et al. 1954) 

e/k = 199° K, 52 A, C:H,, 


a = 4, 
ese = bs KK, a = 3.99 A, CoH 4-Osv. 


(33) 


This corresponds to Vo/kT = 0.29 for T = 293° K and hence, using Fig. 1, 


ik tae as 

f(VofkT) _ 122, _ gon 
a(A) (3.99) 

which is a factor 12 smaller than that given by the interpretation of 7, in 

equation (32). 


(34) 
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A bold conjecture at this point is to attribute a good part of the factor of 12 
difference between (32) and (34) to a drastic breakdown of the empirical 
combining laws given by equation (27) for the case of the C2H4—O2 interaction 
because C2H, is an unsaturated hydrocarbon. 

It may be that the discrepancy by a factor of approximately 2.5 between 
the f/a? values obtained for the CH,—Oz interaction from 7, measurements 
(26) and the use of the empirical combining laws in equations (27), (28), and 
(29) can be attributed to a systematic error arising from the use of an un- 
realistic square well potential. Even if that is the case, we are still left with a 
factor of approximately 5 in f/a? to be explained for the C2:H4—Oz interaction. 

If the range of the interaction given by a = 3.99 A in (33) is assumed to be 
approximately correct, the factor of 5 mentioned above would indicate that 
the correct Vo for the CH,-O; interaction is given by the solution to 


(35) f(Vo/RT) = 6, 
which corresponds to Vo/RT = 2.2 and hence, 
(36) Vo/k = 644° K, e/k = 1150° K. 


On the other hand there is reason to believe, because of the planar structure 
of CeH,, that many collisions between C2H, and O2 molecules involve a 
distance of closest approach considerably smaller than the value of a given in 
equation (33). This would be consistent with the anomalous results obtained. 

It does not seem worth while to analyze our results further at this time for 
several reasons, some of which have been previously discussed. (1) The 
temperature dependence of So F(y)dy in equation (17) has not been carried 
out for a realistic potential. (2) In the derivation of equation (17), the inter- 
molecular potential was assumed to be isotropic (Oppenheim and Bloom 1961) 
and this is certainly not the case for C2H,. (3) (uj) may be temperature depen- 
dent. (4) The result would indicate a strong temperature dependence for 
T,N, in CoH, with oxygen impurities and no experimental study of this has 
yet been made. 

Ethane (C2E,) Gas 

The results of 7, measurements in ethane gas as a function of density at 
314° K are shown in Fig. 6. Again, 7; is proportional to density in the low 
density region and 


(37) T;/p = 0.23 sec/Amagat. 


The rigid molecule model (equation 12), taking the hard-sphere diameter 
of the CH, molecule as 3.95 A, leads to 


(38) * = 1.42. 
Te 
On the basis of a model of C2H¢s having many rotations of the CH; groups 
between collisions, equation (16) gives 


(39) = 0.35. 


r 
Te 








1108 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 








ss 
ot 
e 
64 
54 OXYGENATED (0.5% 03) 
°/ 
e/ 
// 
4 / 
| f° 
f 
a, 
J 
4 
/ ° 
2 
o 
ow 
” 
KI 
T ay T ae ae 
10 20 30 40 50 60 


P AMAGATS 


Fic. 6. Plot of 7; versus density for C2He gas at 314° Kk. 


It is impossible to make precise statements about whether equation (39) 
is valid without further experiments, since we may be dealing with the case 
in which a fraction of the CH; groups are rotating between collisions. How- 
ever, equation (38) would indicate that at least some internal motion is 
present, since we must always have 7, >r. 

The influence of Os impurities on 7; seem to be of the same order of magni- 
tude as for CHg, as indicated by the results for a sample containing 0.5% O2 
shown in Fig. 6. 
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NOTE ADDED IN PROOF 


Since this paper was written, molecular beam experiments on CH, have 
been performed which indicate that our assumption that the spin rotational 
interaction in methane is negligible as compared with the dipole-dipole 
interaction may not be justified. This point is discussed by C. S. Johnson, Jr., 
and J. S. Waugh in a paper closely related to the work reported here, to be 
published in the Journal of Chemical Physics. We wish to thank Professor 
Waugh for his interesting correspondence on this point. 








OBSERVATION OF THE (b'2}-a'A,) TRANSITION IN O.! 


J. F. Noxon? 
ABSTRACT 


The Q branch of the (0,0) band of the electric quadrupole (b'Z$-a'A,) tran- 
sition in O»2 has been observed at 1.908 uw in the emission spectrum of a discharge 
through O2 and He. By a comparison with the (0,0) atmospheric O2, band 
(b'E$-X*%7), the absolute transition probability for the (b-a) system has been 
found to be 2.51073 sec“, with an uncertainty of a factor of 2. The (0,0) band 
of the infrared atmospheric (a'A,-X*Z7) system of O2 has also been observed 
in emission. Using the observed intensity of the (0,1) atmospheric O2 band in 
the aurora and airglow one may predict that the (0,0) (b-a) band should be 
detectable in a strong aurora if observations are made from high altitude. 


1. INTRODUCTION 


The lowest excited electronic states of O2 are a'A, and b'} (see Herzberg 
1950, p. 446) which lie respectively 0.98 ev and 1.63 ev above the X*27 
ground state. Since transitions between any two of these three states are of 
the type g-g, they are forbidden by the selection rules for electric dipole 
radiation and both excited states are consequently metastable. 

Spin-orbit interaction leads to a magnetic dipole character for the (b'}{- 
X*Y;) and (a'A,-X*Z7) transitions (Van Vleck 1934), and they give rise to 
the so-called atmospheric and infrared atmospheric band systems of O2; the 
(0,0) bands fall respectively at 7600 A and 1.27 uw. These systems, as their 
names imply, appear in the absorption spectrum of the earth’s atmosphere 
and have there been investigated in detail: the infrared atmospheric bands 
by Herzberg and Herzberg (1947) and the atmospheric bands by Babcock and 
Herzberg (1948), as well as earlier workers referred to in their paper. The 
atmospheric O» bands are present in the emission from both the aurora and 
the night airglow (see Chamberlain and Meinel 1954) and have been observed 
as well in laboratory sources (Kaplan 1947; Kvifte 1951, and references in 
the latter paper). Quite recently the first observation in emission of the 
infrared atmospheric system was made by Vallance Jones and Harrison (1958), 
who found the (0,1) band at 1.58 uw in the twilight airglow spectrum. Owing 
to reabsorption by O» at lower heights the strong (0,0) bands of both systems 
do not show up in the emission from the upper atmosphere. 

The absolute transition probabilities for these two band systems have been 
determined by absorption studies. Childs and Mecke (1931) found the Einstein 
A coefficient for spontaneous radiation of the atmospheric bands to be 0.14 
sec~!, while Vallance Jones and Harrison (1958) have obtained an A coefficient 
of 1.9X10-‘ sec—! for the infrared atmospheric system. Both A coefficients 
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refer to emission from v’ = 0 in the upper electronic state. Since both of 
these systems correspond to a singlet-triplet intercombination and depend 
upon spin-orbit interaction, the transition probabilities are considerably lower 
than for ordinary magnetic dipole radiation for which the probability is 
around 10? or 10*sec—! (Van Vleck 1934). 

The remaining transition, between the two excited states, b'S+ and a'Ay, is 
forbidden as electric dipole by its g-g nature and is also forbidden as magnetic 
dipole since AA = 2. It does not violate the selection rules for electric quadru- 
pole radiation, however, and is therefore presumably quadrupole in nature. 
Bands of this system have not hitherto been observed; this is understandable 
for several reasons. The transition is evidently highly forbidden. The strongest 
bands all lie at long wavelengths near 2 yp, as will be shown. Finally, the 
transition does not involve the ground state of O» and therefore cannot be 
conveniently studied in absorption. 

The transition probability for emission of the (b-a) system has been estimated 
by Seaton (1958) to be of the order of 1 sec~!. This estimate, however, was 
largely conditioned by a requirement that the total probability for radiation 
from the O.(b'2}) state exceed the probability for emission of the atmospheric 
bands (0—X) alone, 0.14 sec™'. The requirement arose from other considera- 
tions involving the deactivation of metastable oxygen atoms in the upper 
atmosphere. Seaton’s argument has been criticized by Wallace and Chamber- 
lain (1959), who suggested that the transition probability for (6-a) might well 
be smaller than 1 sec. 

There appears to be no observation of any other electronic transition in a 
diatomic molecule which is primarily quadrupole in nature, although such 
transitions are well known in atoms: for example, the green auroral line of 
atomic oxygen (see Herzberg 1944). On the other hand, many molecular 
electronic transitions are permitted both by magnetic dipole and electric 
quadrupole selection rules; the quadrupole contribution to such transitions 
is expected to be much weaker than the magnetic dipole contribution, how- 
ever, and can be detected only by the appearance of weak satellite branches 
with AJ = +2. The first observation of such quadrupole satellite branches 
has recently been reported by Wilkinson and Mulliken (1957), who studied 
the Lyman-—Birge—Hopfield (a'II,-X'Z}) bands of Nz; they deduce the 
quadrupole contribution to the total LBH intensity to be about 0.15. In a 
pure electric quadrupole transition the branches with AJ = +2 (S and O 
branches) are expected to be comparable in intensity with branches for 
which AJ = +1, 0 (R, P, and Q branches). 

It follows that an attempt to observe the (b-a) transition in Os is of interest, 
not only because of the nature of the transition but also because a simul- 
taneous measurement of the intensity of the (b-a) and atmospheric (b—X) 
bands in emission will permit a determination of the absolute transition 
probability of the former. This paper reports work along these lines, using a 
laboratory discharge source, in which the (-a) transition has been observed 
and studied, together with bands of both the atmospheric and infrared 
atmospheric systems of Oo». 
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EXPERIMENTAL 

The discharge was quite similar to one used by Kvifte (1951) in order to 
study the atmospheric O, bands in emission. In order to achieve a high 
concentration of metastables the cylindrical pyrex discharge tube was of large 
dimensions, about 80 cm long and 8 cm in diameter, and was filled to about 
150 mm Hg pressure with pure helium containing about 0.5% Oe. Aluminum 
electrodes were placed in side tubes at either end of the large tube and the 
discharge was operated by a 5000-volt sign transformer at about 60 ma. It 
was found that a higher oxygen concentration actually led to a reduction in 
the atmospheric O» band intensity, in agreement with the findings of Kvifte 
(private communication). With a prolonged operation of the discharge the 
atmospheric bands became weaker, presumably due to cleanup of the Oz. Some 
residual nitrogen in the helium was found to contaminate the near infrared 
spectrum with first positive bands. By admitting the oxygen first the nitrogen 
could be oxidized by applying a Tesla discharge to the helium filling line; the 
oxides formed were then frozen out in a liquid air trap on the way to the 
discharge tube. 

Most of the observations were made using an infrared scanning spec- 
trometer, equipped with a cooled lead sulphide detector, which had been 
designed for studies of the spectrum of the night airglow (Noxon, Harrison, 
and Vallance Jones 1959). With a Corning 7-56 filter to remove higher orders, 
the 1- to 2-y~ region could be examined in the first order; by using a Corning 
7-69 filter the 7000-A to 1-u range could be studied in the second order, 
without any other alterations. It was found desirable to use an instrument 
of greater sensitivity in the near infrared to study the atmospheric O2 bands 
in more detail; for this purpose Dr. D. M. Hunten kindly made available 
his scanning spectrometer equipped with a DuMont K1292 photomultiplier 
as the detector. Preliminary studies to determine the emission spectrum of 
the discharge in the 3500-A to 9000-A range were made photographically. 

In order to measure the relative intensities of the several emission features 
observed, two low brightness sources, originally designed for auroral and night 
sky work, were employed (see Harrison and Vallance Jones 1957). One low 
brightness source had been calibrated against a black body over the range 
4000 A to 1.2 u while the second had been calibrated over the partially over- 
lapping range 1.0 to 2.0 yu; at 1.2 the two sources were found to agree to 
better than 20%. Wavelengths were determined by using second-, third-, and 
fourth-order Hg lines. 


OBSERVATIONS 

A preliminary calculation of the expected positions of bands belonging to 
the (b-a) system of O» showed that the origins of the (0,0), (1,1), and (2,2) 
bands fall at 1.908 uw, 1.938 u, and 1.969 » respectively. Since the potential 
curves for the two electronic states lie directly above one another, the Franck— 
Condon principle predicts that bands with Av = 0 should be considerably 
stronger than those with Av # 0. Consequently, attention was directed to 
the 1.9- region, and the spectrum of the discharge was scanned from 1.85 
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to 1.95 u at a spectral slit width of 70 A. The scanning was carried out very 
slowly, with an output time constant of 30 seconds, in order to achieve the 
highest practical signal-to-noise ratio. 

Figure 1 shows the resulting spectrum which was consistently found to 
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Fic. 1. Solid line: Observed emission spectrum of the discharge through 0.5% Oz in helium 
at 150 mm total pressure, showing the peak at 1.908 uw. The spectral slit width was 70 A. 

Dotted line: Illustrating the appearance of the spectrum observed using pure helium at 
150 mm pressure. 

Dashed line: Assumed contribution of the Q branch of the (0,0) (b-a) O»2 band. 

The predicted positions of the rotational lines of the (0,0) band between 1.89 and 1.93 u 
are indicated. The numbering corresponds to the quantum number of the lower rotational 
level (J’’). Note the great intensity of the background continuum. 


contain the definite peak at 1.908 1+0.003 u. This will later be identified as 
the narrow and relatively intense Q branch of the (0,0) band of the (b-a) 
system. It will be observed that the peak is superposed upon a strong back- 
ground. The scans showing the peak at 1.908» were all obtained with the 
O.-He mixture at high pressure. When the discharge tube was filled to the 
same total pressure with pure helium the peak was absent, as indicated by 
the dotted line in Fig. 1; similarly, the peak was also absent in the He-Oz 
spectrum when the total pressure was reduced to a few millimeters. The 
background emission, in particular the maxima near 1.925 y and 1.935 u, 
remained in all cases. It is presumably to be associated with helium. 

The peak at 1.908 uw has the shape which would be expected from a narrow 
emission feature scanned with a spectral slit width of 70 A, although the long 
wavelength side of the peak does not stand out in its entirety. In Fig. 1 the 
emission present only in the high pressure He—Oz mixture is delineated by a 
dashed line; this was done by matching the spectrum of pure helium at high 
pressure with that of the He-O, mixture at 1.90 and 1.92 » and subtracting 
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the first from the second. When the spectral slit width was reduced to 35 A 
the 1.908-u peak, although weak, appeared to have a width at half height of 
about 35 A; the true width of the emission feature responsible can therefore 
scarcely exceed about 20 A. 

Under conditions where the 1.908-1 peak was strong, a scan in the second 
order revealed the very much more intense (0,0) band of the (b—X) atmos- 
pheric O2 system, and the neighboring permitted atomic oxygen triplet at 
7774 A. Such a scan is shown in Fig. 2. The same physical slit width was used 


0, (b-X) (0,0) 
R 


O 77744 





7600 7700 7800 


Wavelength A 


Fic. 2. Spectrum showing the (0,0) atmospheric O2 band and the allowed atomic oxygen 
triplet at 7774 A. 0.5% Oz in helium at 150 mm total pressure. The spectral slit width was 


35 A. 


as in Fig. 1, but the spectral slit width is naturally reduced by a factor of 
2. While adequate for a determination of the total intensity of the band, the 
resolution is not sufficient to separate clearly the P and R branches of the 
(0,0) band. Both the atomic oxygen line and the atmospheric band were 
absent in the high pressure discharge through pure helium, as expected; in 
the low pressure He—O2 mixture (0.5% O2) the atomic line was very strong, 
but the atmospheric O2 band was either weak or absent entirely. This is to 
be expected as well, since the band arises from a metastable state which will 
be effectively quenched by rapid diffusion to the wall at low pressure. 

The total intensity of the 1.908-u peak could be estimated to within about 
30% by appropriate subtraction of the helium background, and so could be 
compared with the total intensity of the (0,0) atmospheric band. It was 
found that the intensity ratio of the two emissions remained constant to 
within +30% over a range of experimental conditions which gave rise to a 
factor-of-10 variation in the absolute intensity of both emissions. About 20 
separate determinations of the ratio were made. The variation of absolute 
intensity was brought about by changes in total pressure, O» concentration, 
and discharge current. With the help of the calibrated low brightness sources 
the mean ratio, in photon units, of the intensity of the 1.908-4 peak to the 
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(0,0) atmospheric band was found to be 0.004; bearing in mind all possible 
errors involved, the uncertainty in this value is considered to be no more 
than a factor of 2 either way. 

Attention was next directed to the 1.27-y region in the hope of observing 
the (0,0) infrared atmospheric band in emission. A strong feature at this 
wavelength was observed and is shown as the solid line in Fig. 3. Some varia- 
tion between the intensity of this feature and the (0,0) atmospheric band was 


12,500 12,600 12,700 12,800 12,900 
e 
Wavelength A 


Fic. 3. Solid line: Observed emission at 1.27 « using a 45-A spectral slit width; the dis- 
charge conditions were similar to those in Figs. 1 and 2. 
Dashed line: Synthetic profile of the (0,0) infrared atmospheric O2 band constructed for 


500° K and a 45-A spectral slit width. 


found, but the emission at 1.27 u was absent both in pure helium at high 
pressure and in the He—O, mixture at low total pressure. The total intensity 
of the 1.27-u emission was generally about 5% of the (0,0) atmospheric band. 

Finally, the atmospheric O2 bands were examined at higher resolution using 
the Hunten scanning spectrometer. With this instrument it was possible to 
work with a 2-A spectral slit width and so to resolve the rotational structure 
in the P branch of the (0,0) atmospheric band. From the value of J’ corre- 
sponding to the most intense rotational line a rotational temperature of 500° K 
+100° K was estimated; because of the long life of the b'2} upper state, it 
is to be expected that this rotational temperature is a good measure of the 
kinetic temperature in the discharge. 

A search was made for other bands belonging to the atmospheric O2 system. 
This revealed only the (0,1) band at 8640 A. From the absence of other bands 
with v’ > 0 one may conclude that less than 1 or 2% of the b'S} molecules 
were in v’ > 0. Kvifte, on the other hand, did observe the (1,1) band in his 
laboratory discharge, and found it to be about 0.5% of the intensity of the 
(0,0) band; this indicates a concentration in v’ = 1 of about 0.6%. These 
estimates involve the use of the vibrational transition probabilities calculated 
for the atmospheric bands by Fraser, Jarmain, and Nicholls (1954). In this 
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connection it may be noted that the present observations give an intensity 
ratio for the (0,1) to (0,0) band of 0.05+0.01. This agrees well with the 
ratio of 0.047 predicted by the calculation of Fraser et al., ignoring any possible 
variation in the electronic transition moment. It is somewhat larger than the 
value of 0.02 found by Kvifte in his laboratory study of the atmospheric 
bands. 

Apart from the atmospheric Oz bands, the features at 1.27 and 1.91 u, 
which will be identified as the (0,0) band of the (a—-X) and (b-a) systems 
of Os, and a number of strong permitted atomic oxygen lines, the only other 
emissions due to oxygen noted between 3500 A and 20,000 A were several 
weak bands of the first negative system of O32 in the 5000-A to 6000-A region. 
No trace appeared of the forbidden atomic oxygen lines at 5577 A and 6300 A, 
nor were any bands observed belonging to the forbidden Herzberg system of Oo». 
As noted earlier, there is a strong quasi-continuum present at long wavelengths; 
this became noticeable at about 1.2 4 and appeared to be characteristic of a 
discharge through helium. Spectroscopic traces of N2 in the helium were 
evident by virtue of weak first positive band emission in the visible and 


near infrared. 
DISCUSSION 


The (b'}-a'A,) System of O» 
Figure 4 illustrates the rotational levels in the b'S} and a'A, states and 
shows the transitions to be expected; because of the apparently quadrupole 


+0 2 4 J=K 





aAg 


Fic. 4. Rotational levels of the B'S} and a'A, states of Os, showing transitions for 
AJ =0, +1, +2. 


2 3 4 5 JK 


nature of the electronic transition one may expect the O and S branches, 
with AJ = +2, to be comparable in intensity with the P, Q, and R branches. 
The wavelengths of the rotational lines in the (0,0) band were calculated 
from the constants for the two electronic states given by Herzberg (1950); 
for a temperature of 500° K most of the intensity will come from rotational 
lines with J’ < 20. Figure 1 shows the predicted positions of a number of the 
rotational lines expected; the numbering is by J” in the lower state.The 
calculation shows that lines in the Q branch with J’ < 20 all lie between 19,084 
and 19,103 A, a spread of less than 20 A. Hence, with the resolution employed 
in the present work, the Q branch will be indistinguishable from a very narrow 
atomic line; this is a consequence of the nearly identical values for Bo in the 
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two electronic states. A second consequence is that the P and R branches lie 
more or less symmetrically on either side of the Q branch, as do the S and 
O branches; the spacing of the rotational lines changes very slowly with in- 
creasing J. As a result, the lines of the R branch up to J’ = 20 extend over 
about 120A, those in the P branch over about 200A, and those in the S 
and O branches over about 300A and 400A respectively. Finally, Fig. 1 
shows that within a 70-A range of 1.908 uv there fall, at shorter wavelengths, 
only the first two R lines, and, at longer wavelengths, only the first two P 
and the first O lines. Consequently, the Q branch should not be blended to 
any extent at a spectral slit width of 70A or less. 

The appearance of the peak at 1.908 u, together with the constant ratio 
of intensity to the (0,0) atmospheric band, makes it highly probable that 
the feature is in fact the Q branch of the (0,0) band of the (b-a) system. 
Because the other four branches are broad and unresolvable in the present 
work, no identification of them can be made in the presence of the strong 
helium background emission. A spectral slit width of 20A might resolve the 
stronger O branch lines near 1.93 uw; a resolution at least 3 or 4 times greater 
will be required to untangle the P, R, and S lines. An identification of S or 
O lines having an intensity comparable with lines of the P, Q, and R branches 
would be of great interest, and should definitely confirm the electric quadrupole 
nature of the (d-a) transition. Intensity considerations appear to place such 
a high resolution beyond the capability of the present spectrometer and dis- 
charge. 

Since other branches, as well as the helium background, are present in the 
region on either side of the Q branch, it is clear that a subtraction technique 
can only approximate the true intensity of the Q branch. This is not expected 
to yield more than a 30% error, and this has been included in the estimated 
error of a factor of 2 sane with the ratio of the Q branch to the (0,0) 
atmospheric band. It will be recalled that this ratio was found to be about 
0.004. In order to estimate the absolute transition probability for the (0,0) 
band of the (d-a) system the assumption is made that the Q branch repre- 
sents one fifth of the total band intensity. This is approximate, for no attempt 
has been made to calculate individual line intensities. Recalling that the 
Einstein A coefficient for the (0,0) atmospheric band is 0.14 sec, one finds 
that the A coefficient for the (0,0) band in the (b-a) system is about 2.5 107% 
sec-!, with an uncertainty of a factor of 2. 

This probability is far lower than the value of 0.88 sec! assumed by Seaton 
(1958), but it is not an unreasonable value for electric quadrupole radiation 
near 2.0 uw. According to Van Vleck (1934), pure quadrupole emission at about 
7500 A has a transition probability of about 0.1 sec— (to within a factor of 
10 either way). Since the expression for the quadrupole transition probability 
contains a A~> term, a factor of 100 reduction is expected for quadrupole 
radiation near 2.0 u. The predicted probability is thus of the order of magni- 
tude 10-*sec~!, not too dissimilar to the experimental value found here. 
Finally, because transitions with Av ¥ 0 are expected to be of low probability, 
a negligible error is introduced in equating the A coefficient for the (0,0) 
band to that for all radiation out of v’ = 0. 
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The (0,0) Infrared Atmospheric O2 Band in Emission 

The feature at 1.27 » shown in Fig. 3, falls at the correct wavelength to be 
identified with the (0,0) band of the (a—X) infrared atmospheric system of 
O.; this band has not previously been observed in emission. In order to 
confirm the identification, a synthetic profile of the (0,0) band was constructed 
for a temperature of 500° K, using the wavelengths given by Herzberg and 
Herzberg (1947) and the line strength factors given by Van Vleck (1934). A 
triangular slit function was employed to correspond with the 45-A spectral 
slit width used in obtaining the spectrum in Fig. 3. Because of the long life- 
time of the a'A, state, one expects the rotational levels to be in a Boltzmann 
distribution at about 500° K, similar to those in the b'Z} state. The resulting 
synthetic profile is shown in Fig. 3 as a dashed line and appears to agree well 
with the observed spectrum. This agreement confirms the identification of 
the 1.27-u feature as the (0,0) band of the infrared atmospheric system of 
Os. Knowing the A coefficients for both the (0,0) (a—X) and (b—X) bands, 
one may deduce from the relative intensity of the two bands that the con- 
centration of O.(a'A,) in the discharge was about 35 times that of O2(b'Z+). 

A very crude estimate of the absolute concentration of both metastable states 
may also be made, since the absolute sensitivity of the spectrometer is known. 
Assuming that the metastables fill the volume of the discharge uniformly, one 
finds that the equilibrium concentration of a'A, metastables was about 1% 
of the total amount of O2 present, and the concentration of b'Z¢ about a 
factor of 35 smaller. The estimate is probably good to within a factor of 10. 


Predicted Emission of (b-a) Bands in the Aurora and Night Airglow 

The strong Av = 0 bands of the (b-a) system all lie between 1.9 and 2.0 u. 
In this region atmospheric absorption is complex and nearly complete, due 
to the presence of HO and CO: in the earth’s atmosphere (Mohler et ai. 
1950). The possibility of observing (b-a) emission from the aurora or night 
airglow is therefore remote, so long as the observations are made from the 
ground. On the other hand, observations from a balloon or a high flying 
aircraft should be little affected by such absorption. 

The intensity of (b-a) emission from upper atmospheric sources may be 
predicted on the basis of the intensity observed for the atmospheric bands, 
which arise from the same upper electronic state. All but a few per cent of 
the O.(b'=}) molecules which emit radiation in the aurora and airglow are 
observed to be in v’ = 0 (Wallace and Chamberlain 1959). Most of the (b-a) 
emission should therefore be in the (0,0) band at 1.908 uw. The intensity to 
be expected must be predicted from the intensity observed for the (0,1) 
atmospheric band, for the (0,0) atmospheric band is completely reabsorbed 
by O» in the lower atmosphere. From the present work one may deduce that 
the (0,0) (d-a) band has an intensity equal to 0.02/0.05 = 0.4 times that of 
the (0,1) atmospheric band, measured in photon units. 

Chamberlain (1954) has shown that a sizeable fraction of the (0,1) intensity 
observed from the ground may arise from Oz: lying well below the emitting 
level, as a result of absorption of the (0,0) band and subsequent re-emission 
as the (0,1) band. But more recent discussion by Wallace and Chamberlain 
(1959) indicates that the effect is small, presumably because the O2 at lower 
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altitude is collisionally deactivated before reradiating. Consequently it will 
be assumed that the intensity of the (0,1) band as observed from the ground 
is an accurate measure of the intensity of the band as originally emitted from 
the aurora or airglow. 

The intensity of the (0,1) band in the airglow is observed to be about 
1.5kR (kilorayleigh) in the zenith direction (Dufay 1959) and the corre- 
sponding intensity of the (0,0) (d-a) band will be 0.6 kR. This is nearly 100 
times smaller than the predicted intensity of the OH vibration-rotation bands 
in the night airglow at 1.94 (Chamberlain and Smith 1959). Consequently, 
(b-a) emission is not likely to be detected in the night airglow. 

In the aurora, on the other hand, the (0,1) atmospheric band has been 
observed with an intensity as high as 500kR (Hunten 1958). The corre- 
sponding intensity of the (0,0) band in the (b-a) system will be about 200 kR, 
which should be easily detectable above the OH background if observed near 
the zenith and from well above the earth’s surface. 
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A BALLOON-BORNE SPECTROMETER FOR STUDY OF 
THE AIRGLOW BEYOND 2.0 yz! 


J. F. Noxon anp A. VALLANCE JONES 


ABSTRACT 


An account is given of the design, construction, and flight of a balloon-borne 
infrared spectrometer intended for the study of the night airglow spectrum in 
the 2.2- to 3.5-u region. This region is not accessible from the ground because 
of the strong thermal emission from the lower atmosphere. The spectrometer 
employed a 64X64 mm plane diffraction grating And a liquid-oxygen-cooled PbS 
detector. The spectral slit width was 700 A and the noise-equivalent differential 
brightness was about 2X10~' watt. cm~*. sterad~. micron. The observations 
permitted an upper limit of about 700 kR (kilorayleighs) to be placed on the 
zenith brightness of the 1-0 OH band; this value is very close to that predicted 
by Chamberlain and Smith. No emission was observed between 2.3 and 2.7 u. 
This establishes an upper limit on the zenith brightness of emission features in 
this region at 200 kR, for features narrow compared to the spectral slit width, and 
at about 3000kR/micron for the differential brightness of any continuous 
emission. The general design of spectrometers for this wavelength region is con- 
sidered and a discussion is given of improvements suggested by the results of 
this experiment. 


1. INTRODUCTION 


The exploration, from the ground, of the night airglow spectrum in the 
region beyond 2.0 u has been attempted by Noxon, Harrison, and Vallance 
Jones (1959) and by Moroz (1960). In both these investigations it was found 
that, for wavelengths greater than about 2.5 wu, the thermal emission of the 
infrared-active molecules of the lower atmosphere completely obscured any 
emission bands arising from above 50 kilometers; in particular the fundamental 
sequence of OH vibration-rotation bands in the 3.0-y region were unobserv- 
able. Observations from the ground are also restricted to certain narrow trans- 
mission ‘‘windows’’ in the absorption spectrum of the lower atmosphere. 
Consequently, in 1957, we proposed making observations with an infrared 
spectrometer flown by balloon to a height above 15 kilometers. This program 
was carried out as a joint project by the University of Saskatchewan and 
the Canadian Armament Research and Development Establishment of the 
Defence Research Board of Canada (CARDE). The University of Saskatche- 
wan was responsible for the design and construction of the spectrometer and 
its electronics while CARDE was responsible for flying the spectrometer, for 
the telemetry of the results, for the power supplies of all the air-borne equip- 
ment, and for the provision of facilities for preflight testing of the instrument. 

The spectrometer, designed for observations between 2.0 and 4.0 4, was 
flown on April 27, 1960. Good spectra were obtained up to a height of 7.5 km, 
but above this haieh signals received on all the telemetry channels became 

1Manuscript received April 24, 1961. 
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noisy. While the spectra obtained were of considerable interest, and will be 
discussed below, the main purpose of this paper will be to describe the spectro- 
meter and the special problems encountered in this wavelength region. The 
investigation is being continued using a larger spectrometer constructed at 
CARDE, 

2. DESIGN OF SPECTROMETER 


In this section a short summary is given of the basic principles controlling 
the design of infrared spectrometers using photoconductive detectors. When 
the spectrometer employed in the present work was designed, it was assumed 
that the weight and size should be held to a minimum, so that it was necessary 
to consider rather carefully the minimum size of instrument required to detect 
the Av = 1 OH bands. Now that it is possible to fly larger instruments, the 
relations developed in this section are useful for predicting their performance 
on the basis of the results obtained in this investigation. 








4 





seed | 


Fic. 1. Scale drawing of balloon-borne spectrometer. The components are labelled in 
accordance with the description given in Section 3 of the text. 


For a spectrometer of the type shown in Fig. 1 the power reaching the 
detector is given by 
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(1) P = a(By-Ad).w.t 


where a is the area of the detector, B, is the differential brightness of the 
source, w is the solid angle over which the detector receives radiation from 
the condensing optics, and ¢ is the transmission of the optical system. 
(B,:t- Ad is the effective brightness of the surface seen by the detector through 
the spectrometer optics.) 

The noise-equivalent power for many detectors is given by the relationship 


(2) P,, a S.a'2, (Af/f)!2 = k.qi? 


- 


where S is the index of detector performance introduced by R. C. Jones,* Af 
is the band width of the amplifying system, f is the modulation frequency of 
the chopping system, and k = S.(Af/f)!”. 

In describing the performance of a spectrometer it is useful to define the 
quantity, “noise-equivalent brightness’’, Byy. This is the differential bright- 
ness of a surface which would give a signal at the output of the spectrometer 
equal to the r.m.s. noise in the output. From equations (1) and (2) 


(3) Byy = k/(a!?.Ad.w.#). 


However, a must be chosen so that the desired spectral slit width is obtained. 
The width of the exit slit, s, is given by the equation 


s = F.(d0/dd). AX 


where F is the focal length of the condensing lens (Fig. 1) and d@/dd is the 
angular dispersion of the grating. 

The length of the detector, /, is equal to aF, where a is the angular length 
of the detector. The value of a@ is limited by aberrations in the optical system. 
From these relations one obtains the result that 


(4) a = F*.a.(d6/dX).AX. 


Combining equations (3) and (4) one obtains the final result that the noise- 
equivalent brightness is given by the equation 


(5) Byy = S(Af/f)!2.[F. foc. (d8/dd)}!. cot(Ad)32]-1. 


Equation (5) is useful both in the design of spectrometers and in making 
comparisons between existing instruments and new design proposals. 

The design adopted for the balloon-borne spectrometer was checked with 
the help of equation (5) from which it was predicted that a signal-to-noise 
ratio of about 15 should be obtained for the OH bands between 2.8 and 
2.95 uw. The quantities substituted into equation (5) are listed for convenience 
in Table I and correspond to the design described in the next section. The 
apparent zenith differential brightness of the sky in OH emission was calcu- 
lated by constructing synthetic spectra for a rotational temperature of 220° K 
at spectral slit widths of 250 A and 700 A. These synthetic spectra, shown 


*See Smith, R. A., Jones, F. E., and Chasmar, R. R. 1957. The detection and measure- 
ment of infra-red radiation (Oxford), p. 271. 
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TABLE I 
Spectrometer constants and calculation of noise-equivalent brightness 

Spectral slit width, Ad 650 A 
Condenser focal length, F 6.0 cm 
Condenser solid angle, w 1.14 steradian 
Angular length of detector, a 0.1 radian 
Jones’ factor of detector for monochromatic 

radiation near 3.0 pu, S 2.89 X 107° watt. cm7 
Transmission of optics, ¢ 0.2 
Dispersion at detector, d0/dd 3.07 X 108 radian. cm 
Amplifier band width, Af 0.8 sec™! 
Radiation chopping frequency, f 90 sec 
Noise-equivalent brightness, Bw 6.9X10-* watt. cm~?. sterad~. micron=! 
Predicted sky brightness (AA = 700 A, 

dX = 2.95 uw for zenith angle 75°), By 1.0X 1077 watt. cm~*. sterad~!. micron7! 
.”. Signal/noise ratio expected 14.5 








in Fig. 2, are similar to the one calculated by Noxon, Harrison, and Vallance 
Jones (1959); in the present case, however, no correction for atmospheric 
absorption was made. The spectra of Fig. 2 were calculated using the relative 
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Fic. 2. Synthetic OH spectra at spectral slit widths of 700 A (dotted line) and 250 A 
(solid line). The relative brightnesses of the bands in the Av = 1 sequence are approximately 
those given by Chamberlain and Smith (1959) and the brightness of the 1-0 band is chosen 
to be 700kR (4.41078 watt. cm~?. column™), 


band intensities of Chamberlain and Smith (1959), on the basis of an absolute 
brightness of 60 kR for the 4-2 OH band in the zenith direction. This value 
for the absolute brightness of the 4—2 band represents the latest measurements 
of the absolute brightness carried out by Shemansky and Vallance Jones 
(1961) and Moroz (1960) and is lower than the value used by Noxon, Harrison, 
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and Vallance Jones. The absolute brightnesses used for the Av = 1 OH bands 
in calculating the spectra of Fig. 2 were about 30% lower than the values 
predicted by Chamberlain and Smith. The value of the Jones factor, S, was 
calculated from the test data supplied by Infrared Industries with the detector 
used. While there may be some uncertainty with regard to some of the values 
listed in Table I the results of the calculation suggested that a spectrometer 
of the type considered in the next section should be suitable for initial studies 
in the 3.0-u region. 


3. DESCRIPTION OF SPECTROMETER 


The main optical and mechanical features of the spectrometer are shown 
in Fig. 1. Radiation from the sky enters the spectrometer through a hole in 
the side of the case and falls on the mirror chopper, C, from which it is reflected 
on to the entrance slit, S. It then falls on the spherical collimator mirror, M, 
from which it is reflected to the grating, G. The spectrum is focused by an 
arsenic trisulphide lens, L, at the detector, D, which serves also as an exit 
slit. The reflecting chopper, C, provides a means of comparing the radiation 
from the sky with the radiation from the black-body cavity, B, which is 
focused on the entrance slit by the spherical mirror, N. The black-body 
cavity is cooled by dry ice and consequently the signal obtained at the detector 
would, in the ideal case, become a measure of the emission from the sky. 
However, in practice the mirror chopper and the mirror, N, have emissivities 
which are between .02 and .05 and consequently they will radiate as gray 
bodies. To the extent that the emissivities and temperatures of the mirror 
chopper and the mirror, N, differ, a false signal will be contributed. 

A plane, 300 line/mm, diffraction grating, blazed at 3.0 4, was employed. 
This grating was rotated by the cam and lever arrangement, W, so as to scan 
a wavelength range between 1.6 and 3.7 with a repetition period of 80 
seconds. Wavelengths shorter than 2.0 » were removed by a germanium filter 
placed over the entrance slit. The arsenic trisulphide lens, L, had a diameter 
of 75 mm and a focal length of 60 mm. This lens, produced by the Fraser 
Glass Company, had four spherical zones giving the best approximation to 
an ideal plano-convex aspherical lens. The entrance slit, S, and the size of 
the detector, D, were chosen to provide a spectral slit width of 650 A at 
3.0 u. Laboratory observation of the infrared helium line at 2.06 u showed 
that a spectral slit width of 700 A was attained. The detector was an Infrared 
Industries extended-response PbS cell. It was cooled by a hollow finger, P, 
projecting into the cell from a tank of liquid oxygen, T. This tank, which was 
insulated by 2 inches of styrofoam, contained sufficient oxygen to last for 
4 hours. Since the pressure at 27 km is in the neighborhood of 14 mm Hg, 
an intolerable variation in the temperature of the detector would occur if 
the liquid-oxygen container were open to the atmosphere. Consequently the 
oxygen container was provided with a simple spring-loaded valve which 
maintained a pressure differential of about one atmosphere between the out- 
side and the inside pressures. In this way, the variation of the boiling point 
of the liquid oxygen was confined to a small temperature range. The rotating 
chopper disk and the cam of the wavelength scanning system were driven by 
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miniature 400 c.p.s. synchronous geared motors. In order to keep moisture 
from condensing on the window of the detector the inner chamber was filled 
with nitrogen gas which was allowed to boil off from a second insulated 
container filled with liquid nitrogen. 

The output from the detector was amplified by an electron tube preamplifier 
followed by a tuned, 90 c.p.s., transistor amplifier. Three outputs from the 
transistor amplifier, having relative gains of 1, 9, and 80, were fed to non- 
synchronous diode detectors; the outputs from the two most sensitive channels 
were also applied to diode-bridge phase-sensitive detectors for which the 
reference signals were derived from an optical pickup on the chopper, C. 
The outputs from these five detectors were telemetered to the ground together 
with information on the position of the grating, the detector temperature, and 
the temperature of the inside of the spectrometer case. Values of the ambient 
air temperature and pressure were also telemetered to the ground during 
the flight. 

The spectrometer was mounted on the gondola so that its line of sight was 
inclined at a zenith distance of 75°. The azimuth angle of the line of sight 
was free to vary with the rotation of the balloon. It is known that variations 
in the inclination of the gondola during flight do not exceed 5°; the period 
of rotation of the gondola may be as short as 4 seconds. 

The sensitivity and performance of the spectrometer were tested before 
flight by operating the spectrometer with the entrance hole in the outer case 
closed. Under these conditions the spectrometer sees radiation corresponding 
to an ideal black body at the temperature of the inside of the outer chamber. 
In this way the minimum detectable differential brightness at 3.04 was 
estimated to be 2X 10-8 watt. cm~?. sterad~!. micron~! at a spectral slit width 
of 700 A. The difference between this observed value and the predicted 
design value of 0.7X10-° watt. cm~*. sterad~!. micron™! is reasonable when 
additional noise caused by the telemetry transmitter is considered. A check 
was made by inserting temporarily a thermal shield around the motor to 
insure that radiation from the somewhat warmer chopper motor did not 
contribute noticeably to the signal during this calibration. 


4. RESULTS AND DISCUSSION 

The spectrometer was flown from Quebec, Que., on April 27, 1960; a 
polyethylene balloon filled with helium was used to lift the gondola. The 
balloon was launched at 0245 hours E.S.T., reached a height of 7.5 km at 
0305 hours, and attained a maximum height of 27.5km at 0445 hours. It 
was cut down a little later and was recovered in good condition in November, 
1960. Spectra were obtained up to a height of 7.5 km; above this height the 
spectral telemetry channels became noisy. 

Figure 3a shows two sample scans on the intermediate sensitivity channel, 
one on the ground (labelled 1) and the other (labelled 2) at the maximum 
height of 7.5 km. The scan at 7.5 km as obtained on the highest sensitivity 
channel is shown as curve 1 in Fig. 3b. There is also shown as curve 2 on 
Fig. 3b the predicted OH spectrum scaled from Fig. 2 for a zenith angle of 
75° and the measured spectrometer sensitivity. 
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Fic. 3. Specimen spectra as recorded by telemetry. 
(a) (1) Spectrum at ground level. (2) Spectrum from 7.5 km. The zenith angle of observa- 
tion was 75° and the records were obtained on the intermediate sensitivity channel. 
(6) (1) Spectrum at 7.5 km recorded on the highest sensitivity channel. (2) Synthetic OH 
spectrum scaled to match observed spectrum between 2.7 and 2.9 yu. 


Most of the data required for discussing the results have been brought 
together in Fig. 4, which shows the spectra obtained from the ground (B) and 
from 7.5 km (E) plotted on a logarithmic scale. The predicted synthetic OH 
spectrum has been plotted on the same figure as curve F while curves A and 
C are calculated black-body spectra for temperatures of 280° and 250° K 
respectively. Curves D and G show the spectra of a gray body of emissivity 
0.04 at temperatures of 290° and 250° K respectively. The relation of these 
curves will be discussed in the following paragraphs. 

The spectrum obtained on the ground is similar to that observed by Noxon, 
Harrison, and Vallance Jones (1959) and explained by them as thermal 
emission from the lower atmosphere. The brightness at 2.7 uw is expected to 
be close to that from a perfect black body at the ambient air temperature; 
Fig. 4 shows that this is true. During the ascent, the emission at wavelengths 
shorter than 3.2 » decreased in brightness, falling by a factor of 10 by the 
time the spectrometer reached 7.5km. At longer wavelengths the decrease 
in brightness with height became progressively less until at 3.6 u the signal 
was observed to be nearly as strong at 7.5km as on the ground. It must 
therefore be concluded that a false thermal emission from the mirror chopper 
made a significant contribution to the brightness observed at these longer 
wavelengths. Indeed there is a possibility that a large part of the signal 
observed at 7.5km may have arisen from the spurious emission from the 
mirror chopper since, as is shown by curve D in Fig. 4, a gray body having 
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Fic. 4. Comparison between observed, synthetic, and thermal spectra. (A) Spectrum of 
black body at 280° K. (B) Observed night airglow spectrum from ground level. (C) Spectrum 
of black body at 250° K. (D) Spectrum of gray body of emissivity 0.04 at 290° K. (E) Ob- 
served night airglow spectrum at 7.5 km. (F) Predicted synthetic OH spectrum at a spectral 
slit width of 700 A, replotted from Fig. 2. (G) Spectrum of gray body of emissivity of 0.04 at 
250° K. The positions of the Q branches of the first four OH bands are shown. 


an emissivity of .04 and a temperature of 290° K would have a spectral 
distribution of brightness matching the observed spectrum at 2.8 u and from 
3.3 to 3.64. (The chopper was known to have a reflectivity of about 0.95 
and hence an emissivity of about 0.05.) 

The most interesting conclusion results from the comparison in Fig. 4 
between the observed spectrum at maximum altitude, curve E, and a syn- 
thetic OH spectrum, curve F. This synthetic spectrum was scaled from the 
700 A curve in Fig. 2 by multiplying by a van Rhijn factor of 4. As will be 
seen from Fig. 4 the measured brightness at 7.5 km has fallen very close to 
the expected values of brightness for the OH emission and, consequently, 
one may conclude that the actual brightness of the OH bands during the 
flight did not in fact exceed the calculated values. Thus one may put an 
upper limit of approximately 4.4X10-° watt. cm-?. column=!, or 700kR, 
on the brightness of the 1-0 band. 

A further result may be deduced immediately from the observed spectra 
concerning upper atmospheric emissions in the region between 2.3 and 2.7 yu. 
Since no emission features were observed in this region it may be concluded 
that there are present no discrete night airglow features with a brightness 
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greater than about 200kR (~10~-° watt. cm~*. sterad—). This unit corresponds 
to the brightness of a monochromatic feature which would give a detectable 
signal. The corresponding limit to the differential brightness of a broad feature 
is about 2X10-* watt. cm~*. sterad~. micron~! (~3000 kR/micron). 

As a result of the information obtained in this flight it is possible to pre- 
scribe with fair confidence the requirements for an instrument capable of 
making good measurements of the Av = 1 OH bands. The observed spectrum 
at 7.5km and the predicted spectrum shown in Fig. 3b indicate that the 
sensitivity of the spectrometer was adequate for the detection of the OH 
bands at their predicted brightness. There are still two directions in which 
improvements may be made. First, it would be desirable to cool the surface 
of the mirror chopper after the balloon has risen above 20 km. Secondly, the 
spectral slit width might be decreased to about 250 A. Both of these changes 
would decrease the unwanted effect of stray thermal radiation from the 
mirror chopper. As may be seen from Fig. 4 a mirror chopper with a residual 
emissivity of 0.04 at 290° K produces an emission of the same order as the 
expected signal. If the surface of the chopper were reduced to a temperature 
of 250° K by means of a stream of cold dry nitrogen, the curve G in Fig. 4 
shows that the thermal radiation level would be reduced well below that of 
the expected OH signal at 3.0 nu. 

The suggested reduction of the spectral slit width would also be beneficial 
in reducing the signal from continuous thermal radiation and would provide 
more certainty in the identification of the OH bands. Figure 2 shows synthetic 
spectra giving the apparent differential brightness of the bands at spectral 
slit widths of 700 A and 250 A. It will be seen that the structure of the band 
system is much more fully developed for the smaller spectral slit width. How- 
ever, if it is desired to maintain the signal-to-noise ratio at the same value 
obtained with the spectrometer described here then the size of the spectro- 
meter must be increased. In equation (5) the only quantity which is altered in 
scaling up the spectrometer to a larger size is the focal length of the condensing 
system, F, since F is proportional to the linear dimensions of the instrument. 
An inspection of Fig. 2 shows that at the peak of the 2-1 band at 2.95 » the 
apparent differential brightness is increased by nearly a factor of 2 in reducing 
the spectral slit width from 700 to 250 A. The corresponding decrease in 
(Ad)? in the equation (5) is therefore partially compensated for by the increase 
in B,. It is thus only necessary to double F, and hence to double the size of 
the instrument, in order to hold constant the signal-to-noise ratio for the 
peak of the 2-1 band. As was mentioned in Section 1 the experiment is soon 
to be repeated at CARDE using a larger spectrometer fulfilling these specifi- 
cations. Consequently it is not planned to make another attempt with the 
instrument described in this paper. 
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MESURE DU TEMPS DE RELAXATION D’UNE VIBRATION 
DE N.O PAR UNE METHODE INFRAROUGE! 


GERMAIN GAUTHIER ET JULES MARCOUX 


RESUME 


Dans une cuve d’absorption infrarouge permettant un débit de gaz, on produit 
un déséquilibre thermique entre les degrés de liberté du gaz. On enregistre les 
changements qui apparaissent dans la courbe de transmission de la bande de 
N.O a 4.55 uw en fonction de la température et du temps, sous une pression de 
10 mm de mercure. On obtient deux distributions de température: une pour la 
vibration et une pour la translation; la température de vibration est continuelle- 
ment en retard sur la température de translation. Des calculs fournissent un 
temps de relaxation variant entre 107? et 107? seconde, a la pression de 
l’expérience. 


INTRODUCTION 


Quand un gaz subit un changement de température, l’usage courant du 
terme ‘‘température’’ veut dire que l’énergie cinétique de translation des 
molécules du gaz a changé. Une nouvelle distribution statistique des vitesses 
s'est établie dans un temps correspondant a quelques collisions. En plus de 
cet échange d’énergie entre les molécules du gaz, il y a aussi échange d’énergie 
entre les différents degrés de liberté de la molécule: translation, vibration et 
rotation. L’étude de ces échanges d’énergie a |’intérieur de la molécule est 
intéressante 4 un double point de vue: compréhension du comportement des 
chaleurs spécifiques et influence sur les problémes de combustion et de 
dynamique des gaz. 

On appelle phénoméne de relaxation, le retard de l’échange d’énergie entre 
les degrés de liberté externes (translation) et internes (vibration et rotation), 
ou le temps qu’une distribution énergétique des niveaux de vibration et de 
rotation met a s’ajuster 4 une nouvelle température de translation. Cet 
ajustement peut s’effectuer par deux processus: par émission ou absorption 
de radiation, ou par échange d’énergie au cours de collisions. La relaxation 
de translation et celle de rotation ne prennent ordinairement que quelques 
collisions (entre 1 et 10) correspondant a un temps de I’ordre de 10~* seconde 
aux conditions normales de température et de pression (Herzfeld 1955). La 
rotation et la translation peuvent donc en pratique étre considérées toutes 
deux comme des degrés de liberté externes. 

Il a été calculé (Lukasik 1956) que le temps de relaxation de vibration pour 
une désactivation par radiation est constant, indépendant de la pression et 
relativement long, de l’ordre de 0.1 seconde. D’autre part, la désactivation 
par chocs demande plusieurs milliers de chocs; le temps de relaxation corres- 
pondant est inversement proportionnel 4 la pression et varie entre 10-* 
et 10-* seconde. 

1Manuscrit recu le 20 février 1961. 

Contribution du Département de Physique du Collége Militaire Royal de Saint-Jean, 


Saint-Jean, Qué. Travail effectué avec l’aide d’un octroi du Conseil de Recherches pour la 
Défense (DRB 9510-22). 
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Une description des méthodes utilisées pour déterminer les temps de relaxa- 
tion peut étre trouvée facilement (par exemple, Herzfeld 1955). La méthode 
la plus connue vient de l’étude de la dispersion et de l’absorption des ultrasons. 
La méthode proposée ici repose sur l’observation des changements qui apparais- 
sent dans le spectre infrarouge d’une vibration moléculaire isolée, quand un 
gaz subit un déséquilibre thermique. Elle permet des variations de pression 
et de température et différe de la plupart des autres méthodes en ce qu’elle 
mesure le temps de relaxation d’une vibration isolée d’une molécule poly- 
atomique. On a pu ainsi déterminer le temps de relaxation de la vibration 
v3 de la molécule N.O a une pression relativement basse, ot la désactivation 
par chocs est encore prépondérante. 


MONTAGE EXPERIMENTAL 


Un diagramme de l'appareil est représenté sur la figure 1. La cuve d’absorp- 
tion infrarouge est constituée de deux cylindres rectangulaires d’aluminium, 
séparés par une bande de téflon de 13.5 mm d’€paisseur. Le cylindre a une 
section intérieure de 2.54 cm par 50.8 cm. On utilise du NO en cylindre 
commercial. On produit une circulation de gaz par pompage a travers la cellule 
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Fic. 1. (a) Cellule d’absorption. (6) Montage expérimental, vu de plan. 
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avec un débit d’environ 1200 litres 4 la minute, au moyen de trois pompes 
mécaniques de type Edwards S 450. On maintient facilement dans la cellule 
une pression de 1 cm de mercure; cette pression est mesurée a un bout de la 
cellule a l’aide d’un manométre anéroide Edwards. La vitesse d’écoulement du 
gaz est obtenue, en mesurant le temps d’écoulement d’une quantité connue 
de gaz dans les conditions de |’expérience. Dans I’expérience décrite ici, la 
vitesse d’écoulement est de 2.90 métres a la seconde. 

La partie B de la cuve est maintenue a la température ambiante par une 
circulation d’eau prés du téflon. Les deux parois horizontales de la partie A 
peuvent étre chauffées jusqu’a 200° C a l’aide de deux éléments électriques 
placés dans une pate d’amiante durcie. Dans |’expérience, les sections A et B 
sont maintenues respectivement a 150° et 30°C, produisant un gradient de 
90 degrés centigrade par centimétre dans la paroi de la cuve, a l’endroit du 
téflon. La température élevée (150° C) est maintenue constante par un ther- 
mostat branché sur les éléments de chauffage et constitué d’un thermocouple 
et d’un relai Pyr-o-Vane Honeywell-Brown, dont la réponse est de +5 degrés 
centigrade. Le gaz est préchauffé dans une chambre annexe remplie d'un 
treillis métallique fin, assurant ainsi une distribution uniforme du gaz. La 
température de la paroi chaude, mesurée a différents endroits, est uniforme 
a +2 degrés centigrade. 

On examine |’état du gaz a différents points de la cuve, a travers cing paires 
de fenétres de sel gemme, une paire dans la partie A pour enregistrer le spectre 
a 150°C et quatre paires dans la partie B pour enregistrer le spectre a des 
distances connues du téflon. La source et le monochromateur d’un spectro- 
métre Perkin-Elmer 12C sont séparés et placés de part et d’autre de la cuve. 
Le faisceau infrarouge est rendu paralléle a l’aide d’un miroir concave. Ce 
faisceau paralléle passe a travers les différentes paires de fenétres de la cuve, 
au moyen de miroirs inclinés 4 45° pouvant se déplacer sur des rails paralléle- 
ment a la cuve. Le faisceau est ensuite repris pour étre focalisé sur la fente 
d’entrée du monochromateur par un miroir sphérique semblable au premier. 
Le faisceau est limité a l’entrée de la cellule par une fente d’environ 3 mm de 


largeur. 
RESULTATS ET DISCUSSION 

Les spectres de la bande de vibration v3 de N2O apparaissant 4 2223 cm7! 
(4.55 w) ont été enregistrés aux positions suivantes: (1) du cété chaud, prés 
de la bande de téflon, donnant le spectre a 150° C; (2) du cété froid a 6.2, 
13.1, 19.5, 27.1 cm du téflon; (3) 4 l'une quelconque de ces positions, avant 
le chauffage, donnant le spectre a 30° C. 

Un certain nombre de spectres de la bande »; ont été enregistrés dans ces 
conditions. La figure 2 représente la transformation des spectres au cours 
d’une série. Dans cette figure, les courbes expérimentales ont été transformées 
pour présenter la transmission (100 J/Jo) en fonction de la longueur d’onde. 
On a mesuré sur chaque spectre la largeur de la bande (Av en unités arbitraires) 
a 75, 85 et 95% de transmission. Ces mesures apparaissent en graphique dans 
la figure 3. On y voit que la largeur de la bande diminue a mesure qu’on 
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Fic. 2. Exemple de courbes de transmission. 


s’éloigne de la séparation entre les sections chaude et froide de l'appareil; a 
quoi il fallait s’attendre puisque le gaz se refroidit continuellement. Rien 
n’indique cependant que la largeur de bande corresponde 4 une température 
oti tous les degrés de liberté de la molécule sont en équilibre thermique, ou a 
une température ot la vibration n’est pas en équilibre thermique avec la 
translation. 

Il faut peut-étre rappeler ici que l’absorption infrarouge est due aux change- 
ments du moment du dipole électrique au cours d’une vibration ou d’une 
rotation. Ces changements varient avec l’énergie ou la température de la 
vibration et sont indépendants de I’énergie ou de la température de la transla- 
tion. Nous pouvons donc supposer (a) qu’un spectrométre infrarouge peut 
étre considéré comme un thermométre pour mesurer les températures de 
vibration; (6) que les changements de la largeur d’une bande infrarouge avec 
la température de vibration ne sont pas affectés par le fait que la vibration n'est 
pas en équilibre thermique avec la translation. 

Afin de déterminer la relation qui existe entre la température de vibration 
et la largeur de la bande, nous avons fait une expérience o0 tous les degrés de 
liberté sont en équilibre a différentes températures connues. Les conditions 
expérimentales sont ici exactement les mémes que celles décrites plus haut. 
La figure 4 représente graphiquement les résultats de cette expérience. Dans 
le domaine des températures étudiées, la largeur de la bande décroit linéaire- 
ment avec la température, 4 75, 85 et 95%. A des absorptions plus basses que 
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Fic. 3. Largeur de la bande a différents endroits de la section froide de la cellule d’absorp- 
tion, lorsque la vibration n’est pas en équilibre thermique avec la translation. 

Fic. 4. Largeur de la bande en fonction de la température du gaz, lorsque tous les degrés 
de liberté sont en équilibre thermique. 
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60% c’est le phénoméne contraire qui se présente: la largeur de bande croit 
avec la température, alors que la largeur prés de 65% d’absorption demeure 
constante. Dans les calculs qui suivent, nous n’avons utilisé que la largeur de 
la bande 4 75%, trés voisine de la largeur 4 mi-hauteur. La hauteur de la 
bande varie avec la température et présente deux branches d’intensités inégales. 
La largeur de bande est ici considérée comme une propriété variable avec la 
température et n’est pas comparée a la “‘largeur 4 mi-hauteur”’ conventionnelle. 

En considérant ces résultats, il est possible de donner une température aux 
largeurs de bandes lorsque la vibration n’est pas en équilibre avec la translation. 
En effet l’expérience demande que le départ 4 0 cm et la fin 4 35 cm de la 
bande isolante de téflon, représentent des conditions d’équilibre thermique. 
La largeur de la bande de déséquilibre doit retrouver, dans ces deux conditions, 
les valeurs de la largeur de la bande d’équilibre. Ce qui permet de comparer 
deux points sur les courbes: le point de départ 4 150° C (largeur de bande 
5.60 unités) et le point d’arrivée 4 environ 30°C (largeur de bande 4.75 
unités). On établit une relation linéaire entre les largeurs de bande et les 
températures de vibration entre ces deux points; ce qui permet d’étalonner en 
températures les largeurs de la bande de déséquilibre. La courbe des tempéra- 
tures de vibration en fonction des distances a |’équilibre initial est tracée sur 
la figure 5 (courbe A). 

La courbe B de la figure 5 représente la variation de la température de 
translation en fonction de la distance au téflon. Cette température de trans- 
lation a été mesurée avec des thermocouples placés 4 différents endroits de 


(A) Température de vibration 


(B) Température de translation 
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Fic. 5. Température de vibration (A) et de translation (B) a différents endroits de la 
cellule d’absorption. 
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l’appareil. Théoriquement, s’il n’y avait pas de convection, on devrait observer 
prés du téflon une variation de température appréciable entre le centre et le 
haut ou le bas de I|’appareil, soit sur une distance verticale de 2.54 cm. Les 
mesures indiquent une différence d’environ un degré centigrade: les courants 
de convection sont donc appréciables. Au cours de l’expérience, nous avons 
constaté que des variations dans le systéme de chauffage produisaient presque 
instantanément des variations dans la température du gaz. 

Les températures mesurées sont limitées dans leur précision par le systéme 
régulateur de température, qui fonctionne 4 +5° C. Cette erreur est indiquée 
sur la courbe B de la figure 5. Cette derniére figure permet de se représenter 
assez bien I’allure du phénoméne. Au temps ¢ = 0, la vibration est en équilibre 
thermique avec la translation. Le gaz se déplace relativement vite (2.90 
m/sec) et la température de vibration 7) ne peut suivre celle de la translation 
T,. La différence 7,—7, des températures est d’abord petite, passe par un 
maximum aprés environ 0.02 seconde et devient nulle aprés 0.12 seconde, 
lorsque l’équilibre s’est rétabli. 

Le temps de relaxation est défini par la relation suivante (Herzfeld et 
Litovitz 1959): 

ie: ad? «ne 

- (7,—T;) 
ot 7, est la température de vibration, 7, la température de translation et 7 le 
temps de relaxation. Pour un phénoméne donné déterminant d7,/dt 4 chaque 
instant, 7 est proportionnel a la différence (7,—7,). La courbe A de la figure 5 
permet de déterminer le rapport 07,/dt, les courbes A et B déterminent la 
différence (7,—T7,). Les différents temps de relaxation 7 ainsi calculés sont 
donnés dans le Tableau I. Pour des pressions assez basses ot la plupart des 


TABLEAU I 


Temps de relaxation calculés en partant des courbes 
de la figure 5 





t (sec) 
TCO) ae te) oT, /at (10 mm Hg) 

119 138 2.61073 0.007 
83 115 2.1 0.016 
67 97 Coe 0.018 
59 82 1.4 0.016 
55 70 1.0 0.015 
51 62 0.7 0.015 
49 56 0.6 0.012 
46 51 0.4 0.011 





collisions sont binaires, le temps de relaxation est inversement proportionnel 
a la pression. Un temps de relaxation de 0.007 seconde 4 10 mm de mercure 
devient donc approximativement 10~* seconde a l’atmosphére. 

Herzfeld et Litovitz (1959) donnent plusieurs valeurs obtenues par dif- 
férents chercheurs. Toutes ces valeurs sont de l’ordre de 10-* seconde, soit 
100 fois plus petites que celles de l’expérience présente. II est difficile de 
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comparer les résultats obtenus ici et ceux donnés par I’absorption des ultrasons. 
Dans ce dernier phénoméne, le temps de relaxation est mesuré 4 une tempéra- 
ture, la température moyenne du gaz, et les variations extrémes de température 
se produisent dans un temps égal a la période de I’onde. Dans la présente 
expérience, les variations extrémes se produisent dans un temps de I’ordre de 
0.02 seconde et l’on peut pointer des températures intermédiaires. 

Nous avons pu vérifier qu’entre 30° et 150° C, la surface sous les courbes, ou 
l’intensité intégrée de la bande d’absorption, ne varie pas lorsque tous les 
degrés de liberté sont en équilibre. Cependant, le déséquilibre thermique des 
degrés de liberté produit, dans l’expérience, une variation appréciable de la 
surface comprise sous les courbes d’absorption. Ces variations devraient 
pouvoir conduire a un calcul de la contribution de la vibration 4a I’intensité 
totale de la bande. Nous espérons continuer ces travaux et les mener dans deux 
directions: (1) en variant la pression, vérifier comment varie le temps de 
relaxation et atteindre si possible le phénoméne de relaxation par émission de 
radiation (Lukasik 1956), (2) déterminer pour une méme molécule, les dif- 
férents temps de relaxation correspondant a différentes fréquences fonda- 
mentales (Massey et Burhop 1952). 
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NEW EMISSION BANDS OF Nj, 21,-A 211, 


Y. TANAKA, T. NAMIOKA, AND A. S. JURSA 


ABSTRACT 


Thirty-eight new bands of nitrogen were observed in the region 2050-3070 A 
in an a-c. condensed discharge of neon mixed with a small amount of nitrogen. 
It was determined that these bands result from the transition *II,-A 7M, of N?. 


INTRODUCTION 


In addition to the two well-known systems of Nt, B?2t-X 2 the first 
negative system and C *2+—X *Z+ the second negative system, Meinel (1950) 
discovered another system of emission bands in the night airglow and assigned 
them to the transition A *II,-X ?2+, the upper state being the one which 
Mulliken (1932) predicted at about 1 ev above the ground state of Nf. 
Shortly after Meinel’s discovery, Worley (1953) studied the absorption spectrum 
of N2in the vacuum ultraviolet region and arranged a Rydberg series (he called 
it the third series) which converged to the A *II, state of Nt. Douglas and 
his co-worker (1951, 1953) have succeeded in producing the Meinel bands in 
the laboratory and have obtained accurate molecular constants for the related 
states. 

Recently, Janin and d’Incan (1958, 1959) reported additional emission 
bands of N¢ in the region 2240-3070 A and assigned them to the transition 
271,-A *II,, the lower state being the upper state of the Meinel bands. Inde- 
pendently and at almost the same time as they reported the discovery of these 
bands we also observed a new system of emission bands of N¢ in the region 
2050-3070 A and temporarily attributed the system to the transition 
211,-A ?II,,* which turned out to be Janin’s assignment for the bands he 
observed. 

EXPERIMENTAL 

We used a z-shaped discharge tube, 2.5 cm in diameter, 30 cm in length, 
with two aluminum electrodes, which was attached to a 2-meter, normal 
incidence type vacuum spectrograph utilizing a grating of 15,000 lines per 
inch. The tube was separated from the spectrograph by a thin LiF window. 
Tank neon, purified by passage through a charcoal trap, was introduced into 
the discharge tube over a pressure range of 1-20 mm Hg. A small amount of 
pure nitrogen was added to the neon; the amount of Nz suitable for excitation 
of the present bands was about 1-5% of the neon pressure. If too much N2 
were added, strong bands, such as the Lyman-—Birge—Hopfield and the fourth 
positive bands, interfered with the identification of the present bands. 

1Manuscript received May 1, 1961. 

Contribution from the Photochemistry Laboratory, Geophysics Research Directorate, 
Air Force Cambridge Research Laboratories, L. G. Hanscom Field, Bedford, Massachusetts. 

*The results with an explanation were presented at the Symposium of Molecular Structure 
and Spectroscopy, Ohio State University, June, 1958. 
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The mixture of Ne+N: was excited by an a-c. condensed discharge produced 
by a 15-kv, 0.96-kva transformer, combined with a laboratory-made condenser 
of about 0.005 uf, and an external spark gap. In most cases we used a flow 
method in which Ne and N» were suitably premixed and then pumped through 
the discharge tube. If one uses a closed system, then one must add nitrogen to 
the tube periodically in order to maintain conditions favorable for the pro- 
duction of the bands; this is necessary because the nitrogen disappears rapidly 
in the discharge in a closed system. The bands were quite strong so that a 
10-minute exposure produced a reasonable density on SWR film in the first 
order with a spectrograph slit width of 20 uw. The use of Ne is the essential 
factor for the production of the bands. With helium the bands do not appear 
at all, but the strong C-X bands and some unidentified very weak bands in 
the same region as that of the present bands do appear. With argon neither 
the present bands nor the C—X bands are excited. The use of an a-c. con- 
densed discharge is the other factor important in the production of the bands. 
In a non-condensed discharge in a mixture of Ne and Ne, band systems of 
N: such as y"II,-a’ '27, y I,-w "Ay, x'2>-a' 'Z7,a M,-X 'ZF, and the Herman- 
Kaplan system were strong throughout this region, while the bands presently 
being considered were extremely weak. 


RESULTS 

The spectrum as observed in the region of 2110-2530 A is reproduced in 
Fig. 1(a@). As one can see, the individual bands are composed of two heads, the 
longer wavelength one being stronger than the other, and both are degraded 
toward the red with relatively coarse spacing. The measurements were carried 
out for the heads, and the results are listed in Table I. The values listed in 
Table I are the averaged ones measured from several plates which were taken 
under different conditions. Janin and d’Incan’s values for the band origins 
are also listed in the same table. Recently, Janin et al. (Grandmontagne et al. 
1959) revised their previous vibration assignment and raised the upper 
vibration number by 3 after calculating the relative intensity distribution 
by means of the modified ‘‘distortion’”’ method (Pillow 1951; Grandmontagne 
1957, 1958). We adopted their revised assignment in Table I and Fig. l(a), 
although our original assignment was exactly the same as their old one. The 
intensity distribution alone cannot be considered as a reliable means of 
determining the vibrational numbering. However, at present, their new 
assignment has a plausible basis for explaining the actual band appearance, so 
it appears to be proper to accept their assignment until further evidence, such 
as isotope shift, becomes available. 

For the revised numbering shown in Table I, the long wavelength heads 
(R:) are well represented by the following equation calculated by the least- 
squares method: 


(1) vy = 43128.44+905.50(0’ +4) —11.69(v’ +3)?+0.0114(0’ +-3)° 
—1902.80(0”’ +4) +14.92(0’ +4)?. 













TABLE I 
Observed wavenumbers of the *II,-A *II,, bands. Janin and d’Incan’s values 


are for the band origins 
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R;, R2 heads 








R> heads 
v’ vo” (cm) I calc. 
10 1 48598 .4 
538.5 3 48539 .9 
9 1 47921 .9 
863.8 4 47864.8 
8 1 47224.4 
166.9 5 47166 .9 
7 1 46504.5 
447.9 6 46446 .3 
9 2 46077 .9 
020.8 6 46021.7 
8 2 45381 .3 
325.4 7 45323 .8 
10 3 44943 .2 
885.1 4 44883 .5 
a 2 44661 .6 
605.0 10 44603 .2 
6 2 43915.9 
857.5 Z 43859 .6 
10 4 43161.9 
102.1 2 43100.1 
5 2 43151.2 
093 .0 4 43092 .2 
Z 3 42863 .8 
788 .6 5 42789 .9 
4 2 42363 .5 
303.3 1 42303 .7 
6 3 42105.8 
045.4 5 42046 .4 
11 5 42053 .7 
41997 .7 6 41998 .9 
10 5 41399.3 
343.3 1 41346.5 
5 3 41335.5 
278.0 5 41279.9 
9 5 40729 .7 
669.9 1 40671.3 
4 3 40545 .2 
489 .6 2 40490. 4 
8 5 40026 .9 
39972 .0 2 39973 .5 
5 4 39552 .4 
495.1 2 39496.5 
7 5 39305. 7 
252.2 Z 39252.8 
4 4 38768 .6 
707.2 3 38707 .0 
8 6 —- 
38247 .0 ] 38249 .7 
3 4 37947 ..5 
893.3 3 37894 .5 
5 5 37801 .6 
742.8 2 37742.9 
7 6 37586 .5 
§28.2 6 37529.1 
4 ‘fF z 
36953 .4 2 36953 .4 
6 6 36844 .9 
786.! 36785 .6 


Janin and 


d’Incan (1958) 


45350. 


44629. 
43886. 
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TABLE I (Concluded) 


R» heads 





Ri, Ro heads Janin and 
v’ vo” (cm?) I calc. d’Incan (1958) 
3 5 — 
36140.5 2 36140.9 36167 
5 6 36079 .2 
020.8 1 36019.1 
9 8 35649 .4 
591.0 15° 35589 .6 
4 6 — 
35232.5 0 35229 .6 
6 7 35158 .0 
094.6 3 35091 .6 35120 
8 8 34951 .5 
894.0 4 34891 .7 34921 .0 
5 7 34386 .0 
325.5 2 34325 .2 
7 8 34227 .9 
170.5 10 34171.1 
7 9 32591 .2 
5 1 32536.8 32567 


532. 








*Janin’s unpublished data. 
‘Superimposed with another band. 


For the old assignment we have the following equation: 
(1’)  vyree = 45739.94+835.66(0' +3) —11.60(v’ + 3)?+0.012(v’ +3)* 
— 1902.81 (v’’ +3) +14.92(v’’ +4)?- 


No attempt has been made to formulate the second heads (R;) because o 
inaccurate measurements resulting from their weak appearance. However, the 
wavenumber differences of the two heads are fairly constant over the entire 
system, the average being Av = 57.8 cm™. This indicates that the equation 
for the second heads will be practically identical with equation (1) or (1’) 
except for a small difference in the electronic term. 

As shown in Table II, the vibration constants for the present lower state 


TABLE II 
The vibrational constants of "II, and A ?II, states. The values of Janin et al. are all based 


on the band origin, while the present values are from the head analysis 


Janin e¢ al. 


Present ee 
SI vo = 0 vo = 3° Douglas 
ve = 0 vo = 3? (1958) (1959) vo = 3° (1953) 
State (cm™) (cm) (cm7?) (cm?) (em~) (cm) 
“II, We 835.66 905 . 50 839.5 919.2 914.8 
WeXe 11.60 11.69 12.75 13.80 13.13 
WeVe 0.012 0.0114 0.116 0.117 0.083 
All, we 1902.81 1902.80 1903.42 — 1903.6 1902.84 
WeXe 14.92 14.92 15.00 — 15.0 14.91 


“Old assignment. 
>Revised assignment. 
*Janin’s unpublished data. 
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agree very well with those of the A II, state of N¢. Therefore, it is certainly safe 
to assume that the lower state is the A *II,, which is the upper state of the Meinel 
bands. Judging from the band features the upper state is also a 7II, and the 
two heads are formed by the R branches of subbands. In this type of transition, 
in general, the Q branch is much weaker than other main branches. With that 
assumption and the fact that the A II, is inverted, the upper state should be 
inverted also, and its spin separation turns out to be approximately 24 cm, 
which is a little less than one third of the spin separation for the A-state. 
The w, of the upper state is much smaller than the w, of the A-state indicating 
that the former state is much less stable than the latter. The linear extra- 
polation of the AG,,3;: v curve of the upper state gives 2.12 ev as a value for 
Dy (1.81 ev for old assignment). Because this value results from a long extra- 
polation it is not necessarily an accurate one; however, interestingly enough, 
it gives 8.52 ev (8.54 ev for old assignment) as the dissociation limit of the 
upper state which is very close to the value of the dissociation energy, 
Dy = 8.66 ev, of the ground state, X *2t, of Nf. 


DISCUSSION 


As mentioned earlier, only a mixture of Nz with Ne is suitable for the pro- 
duction of the bands. Furthermore, a condensed discharge will produce these 
bands with much more intensity than an uncondensed discharge will. Also, 
when the bands are produced, strong C—X bands of N¢ always appear at the 
same time while the N2 bands, which usually appear in the same region, are 
almost gone. These observations establish quite firmly that the bands belong 
to Nj and not to Nz». 

The energy level diagram of N¢ is shown in Fig. 2 together with potential 
curves. Except for the state *II,, the rest of the curves have been duplicated 
from Gilmore (private communication). As stated previously, the upper 
state, *II,, of the present bands very probably dissociates at the same dissocia- 
tion limit with the ground state whose dissociation products are N(2p* 4S) 
+N+(2p?’P). In addition to the three known states, X 22+, A ?II,, Bf, 
nine more states are predicted from these products: *II,, *2%,,, *IIu.,, *Xt.4, 
and ‘II,,,,. However, according to Mulliken’s predictions (1957) these additional 
states are all repulsive, although he did not mention their degree of instability. 
In approximately the same energy region, he also predicted two stable states 
42+ (To = 175,000 cm~ = 21.7 ev)* and 4A, (To = 190,000 cm—! = 23.6 ev),* 
both of which arise from the dissociation products N(2p* 2D)-+N+(2p? ?P).t 
However, it appears obvious that neither of these two states could be the 
upper state of the present bands although one of them, ‘2{, has energy 
which is fairly close to that of the present upper state. The only possible 
explanation is that the upper state has to be the *II,, which arises from the 
products, N(4S)-++-N+(?P), mentioned above. This explanation disagrees with 
Mulliken’s prediction; however, on the other hand, the state is considerably 
less stable (Do ~~ 2.1 ev) than other known states. 

The ionization energy of Ne is 21.47 ev, which is slightly below the energy 


*The energy measured from the ground state of No. 
tGilmore’s state ‘2 (see Fig. 2) dissociates to N(4S) +N+*(8P). 
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Fic. 2. Energy levels and potential curves of N¢. 


of the present upper state. Apparently, the excitation of the bands occurs in a 
collision between Net+ and No». This is very similar to the process for the 
excitation of the C-X bands by Het as suggested by Watson and Koontz 
(1934) and Setlow (1948). The ionization energy of He, 24.580 ev, is greater 
than the dissociation limit, 24.24 ev, represented by the products N(‘4S) 
+Nt(8P) and this is the reason why the present bands were not excited in a 
He atmosphere, although the strong C-X bands appeared with an intensity 
anomaly in v’ = 3 (Douglas 1952; Tanaka 1953; Wilkinson 1956; Carroll 
1959). 

Janin et al. (Grandmontagne et al. 1959) stated that the absence of the 
v’ = 0, 1, and 2 progressions on their plates was partly due to a selective 
excitation and partly due to superposition of these bands with other bands in 
the same system. As mentioned before, the ionization potential of Ne is 
slightly below the v’ = 0 level and there is no state in Ne which may cause 
selective excitation of the present bands at v’ = 3. As for their second reason, 
an attempt was made to identify the bands belonging to the progression of 
v’ = 0, 1, and 2. Some of the bands seemed to be present (see Table III); 
however, they were so weak that we cannot be sure of their identification. For 
reference, the wavenumbers of these are listed in Table III together with 


some others. 
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TABLE III 
Observed wavenumbers of uncertain band heads 





R,, Re heads R2 heads 
v’ vo” (cm~) I calc. 
10 0 — 
50409 .3 0 50412.9 
12 2 48031 .9 
47981 .2 0 47979.1 
12 5 — 
42625 .9 52 42628 .8 
12 8 37605 .3 
37550.8 1 37547 .0 
2 5 _— 
35303 .9 0 35305. 2 
1 5 34498 .6 
34444.2 0 34446 .3 
3 6 —_ 
34414.6 0 34417 .6 
3 7 32770.1 32723 .2 
(or 1 6) 32718.1 0 (or 32722.6) 





“Superimposed with another band. 


It might be interesting to take note of the difference in the relative intensity 
distribution of the C—X bands resulting from the two methods of excitation: 
in a He atmosphere and in a Ne atmosphere. The two spectra are compared in 
Fig. 1(6), which shows the difference more clearly than is shown in the paper 
by Takamine and his co-workers (1939). According to the mechanism proposed 
by them, the bands which originate in the level v’ = 1 must be abnormally 
strong in the case of Ne atmosphere. Actually, there is no conspicuous anomaly 
in the spectrum, and the distribution of intensity appeared normal to us. 
The C-X bands can also be excited in a condensed discharge of pure No». 
Although the bands are very weak in this case, the distribution of intensity 
is almost exactly the same as that when No is excited in a Ne atmosphere. 
These observations lead us to believe that the excitation mechanism proposed 
by Takamine and his co-workers is actually not the case. Further work con- 
cerning this problem is now in progress and additional information will appear 
in the near future. 
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NOTE ADDED IN PROOF 


Our recent study of the isotope shift in the *II,— A *II, bands substantiated 
the fact that the vibrational numbering as revised by Janin ef a/. is correct. 
An analysis of the *II, — A *II, bands of (N}°)* will be presented in a separate 


paper. 





SELF-DIFFUSION IN BODY-CENTERED CUBIC ZIRCONIUM! 


G. Kipson AND J. McGurn 


ABSTRACT 


Self-diffusion coefficients of crystal-bar zirconium have been measured between 
1500° C and 1100°C, using radioactive Zr® as a tracer. The results may be 
represented by 


D = 2.41074 ex cm?/sec. 


= ar 
P RT 
The pre-exponential factor is about three orders of magnitude smaller than that 
measured in most close-packed systems and the activation energy about one-half 
that anticipated from an empirical correlation with the melting point. The results, 
however, are similar to those of a few other recently studied body-centered cubic 
(BCC) systems, and agree quantitatively with work in the Soviet Union on 
zirconium. There is considerable evidence that the diffusion process occurs via 
vacant lattice sites. 


1. INTRODUCTION 


It is well established that the self-diffusion coefficients of pure metals have 
a simple Arrhenius type temperature dependence of the form 


~Drexp-| | 
(1) D = Doexp | 2. ; 


where Dy and Q are temperature-independent parameters characteristic of the 
particular system studied. It has been noted frequently that a rather good 
correlation exists between the activation energy Q and the absolute melting 
temperature for most of the metals studied. (See, for example, Kidson and 
Ross 1958.) It has likewise been pointed out that the values of Do usually lie 
in the range of 0.1 to 10 cm?/sec (LeClaire 1953). These two observations have 
made it possible to predict values of the self-diffusion coefficients for systems 
in which no direct measurements are available. 

Of the various models proposed to account for the fundamental diffusion 
process, the two most discussed are the vacancy and the ring mechanisms. The 
latter, proposed by Zener (1950), is a generalization of the simple interchange 
of two neighboring atoms to a group or ring of m atoms which move syn- 
chronously to produce the effect of a jump for any particular member of the 
ring. While the resulting expressions for D for the two mechanisms are similar, 
and lead to equations of the form (1), a fundamental difference does exist. 
This became significant when Smigelskas and Kirkendall (1947) showed that 
the different rates of interdiffusion of two chemically distinct species can give 
rise to a measurable mass flow in the diffusion zone. The phenomenon is 
readily interpretable in terms of a vacancy mechanism, but it cannot be 
accounted for if a ring mechanism is operative. This conclusion, plus the very 

1Manuscript received February 2, 1961. 
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considerable amount of supplementary evidence accumulated during the past 
decade, has led to a general acceptance of the vacancy model for at least the 
close-packed metal systems. The situation is not quite so well defined for 
BCC systems, but here again the ring mechanism cannot account for the 
observed ‘Kirkendall shifts’’. 

Zener (1951) has suggested that, regardless of the particular model chosen, 
the restricted range of Do is in keeping with the thermodynamic requirement 
that the associated entropy of activation must be positive. During the past 
two years, however, results of diffusion measurements in a number of BCC 
systems have been reported that appear anomalous both in the fact that a 
negative entropy of activation is implied by the Do values, and that the activa- 
tion energies are about one-half those expected on the basis of the empirical 
correlation with the melting points. These results are summarized in Table I. 


TABLE I 


Self-diffusion in body-centered cubic metals 


Do Qobserved (Q calculated 





Authors Metal (cm?/sec) (kcal/mole) (kcal/mole) 

Adda and 

Kirianenko (1958) y Uranium 1.811073 27.5 49 
Bochvar et al. (1958) y Uranium 1.17 XIE? 26.6 49 
Rothman et al. (1959) =y Uranium 2.38xXie- 28.5 49 
Mortlock and Tracer Cr® 

Tomlin (1959)* in 6 titanium 5X1073 35.3 73 
Paxton and 

Gondolf (1959) Chromium 10-* 52.0 76 





*Qualitatively similar results have been obtained recently for a series of dilute solutes in 8 Ti (Lidiard 1960). 


An attempt to account for the low Do values has been made recently by 
Pound, Bitler, and Paxton (1961). In essence, they suggest that an additional 
factor must be incorporated into the previous expression for Do derived on the 
basis of absolute-reaction-rate theory. This factor takes into account the 
requirement that not only must a migrating atom possess sufficient energy in 
the direction of a jump, but also that its actual movement must be in that 
direction or very close to it in order that a successful jump can occur. That is to 
say, the components of its vibration in the plane perpendicular to the jump 
direction must be very small or zero. 

This restriction gives rise to a new factor [1/f,]*” where f, is the Einstein 
partition function for vibration, and n is the number of atoms involved in the 
jump. 

For a vacancy model, ” was assumed equal to unity, and the effect was not 
significant. For a four-atom-ring model, however, m is 4, and the additional 
factor is large enough to account for the observed values of Do in the BCC 
systems. On this basis, Pound et al. conclude that a low Dy in a BCC system is 
evidence for the ring mechanism being operative. They suggest that the large 
chemical concentration gradients present in the Kirkendall experiments could 
give rise to abnormally large vacancy concentrations, and as a result, the 
vacancy mechanism becomes dominant in these systems. The rather 
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extraordinarily low activation energies tabulated in Table I were not discussed. 

In view of the considerable intrinsic interest in the problem, and the tech- 
nological importance in reactor metallurgy, a study of self-diffusion in the 
BCC 8 phase of high-purity zirconium has been made. During the course of the 
present investigation, the authors became aware of similar studies reported by 
three Soviet groups, Borisov et al. (1958), Gruzin et al. (1958), and Lyashenko 
et al. (1959). Their results will be compared with the present work later. 


2. EXPERIMENTAL PROCEDURES AND RESULTS 

Three different techniques were used to measure the diffusion coefficient 
of zirconium as a function of temperature. The first of these was the basic 
sectioning technique described below. In view of the rather surprising results 
so obtained, two additional methods were employed as checks on the first. 
One of these avoided the possibility of so-called ‘‘short-circuit”’ effects, while 
the other was a sintering experiment designed to provide information concern- 
ing the mechanism of the diffusion process. 


Experimental Procedure. I 
Crystal-bar zirconium was obtained from the Foote Mineral Company. The 
analysis is given in Table II. Buttons weighing approximately 80 grams were 


TABLE II 


Analysis of crystal-bar zirconium* 





<0.005°% Ni <0.005% 
Cu Not detected 


Mg <0.005% 
Hf/Zr <0.05% 


Fe <0.005% 
Ti <0.005% | 


: | 
Al <0.005% || Ca 0.02% 





*Received from Foote Mineral Company, Lot No. 810-8. 


prepared by arc-melting sections of the crystal bar in an argon atmosphere, 
using a tungsten electrode.* Specimen disks ranging from ? in. to 3} in. in 
diameter and about 4 in. thick were machined from the buttons, cleaned in 
a chemical dip, and annealed at 1100° C for 4 hours in a dynamic vacuum 
of 3X10-* mm, to remove any absorbed hydrogen. One face of each disk was 
carefully polished mechanically and chemically, using a 50% H:20, 45% HWNOs, 
5% HF solution. 

The radioactive tracer isotope Zr® was prepared by irradiating high purity 
zirconium wire in the NRX reactor at 1.510" n/cm? sec for 14 days. The 
use of the Zr® isotope was complicated by the presence of the radioactive 
daughter decay product Nb", since y-ray energies of both species are so close 
as to make discrimination between the two impossible. The contribution of the 
Nb® to the total activity measured was calculated. The results are shown 
graphically in Fig. 1. It was concluded that chemical separation of the Nb® 
could be avoided justifiably provided the diffusion anneals were completed 
within 4 days of the removal of the irradiated zirconium from the reactor. 
This precaution was observed in all of the runs in this work. 


*Subsequent spectral analysis gave no indication of tungsten contamination. 
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Fic. 1. Calculated curves for the radioactivity of Zr® and Nb® after 14 days’ neutron 
irradiation at 1.510" n/cm? sec. 


A thin layer of the tracer metal was evaporated onto one face of each of the 
two disks, using the apparatus shown in Fig. 2. The diffusion specimens 
were annealed in pairs, in a ‘‘Globar’’ furnace and a dynamic vacuum of 2 X 10-* 
mm of Hg. In order to avoid both contamination of the specimens and loss of 
the radioactive deposit by evaporation, the disks were placed with active faces 
close together in a pure zirconium “gettering boat’’. The boat was held in the 
cold portion of the annealing tube until the furnace was up to temperature, 
and then pulled into the hot zone to reduce the heating time to a minimum. 
The reverse procedure was followed at the end of the anneal. A typical heating 
and cooling curve is shown in Fig. 3 from which appropriate corrections to the 
annealing times were made. 

Temperatures were controlled to +3°C and measured to +2°C, using 
a Pt — Pt Rh thermocouple. 

Following the diffusion anneal, thin sections accurately parallel to the face 
of the specimen were removed, using a special lathe device. Surface diffusion 
effects were avoided in the usual way by first reducing the diameter of the 
specimen by 0.030 inch. The turnings from each section were collected in 
bottles, weighed, and dissolved in a solution of 45 parts HNO;, 50 parts 
H,O, 5 parts HF. The y-ray activity was measured over the 0.75-Mev peak 
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Fic. 2. A schematic drawing of the apparatus used for the evaporation of radioactive 


zirconium onto the diffusion specimens. 
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Fic. 3. A typical heating and cooling curve of the diffusion specimens. 
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using a well-type crystal scintillation counter and a single-channel analyzer. 
Appropriate corrections for background and decay were applied, and sufficient 
counts taken to ensure a statistical error of less than 1%. 


3. RESULTS 


If the diffusion paths of the atoms were through the lattice proper, as opposed 
to the grain boundaries, the resulting distribution should be given by 


(2) OMe) = is ew [35 


where C is the concentration of the diffusing species (in this case, the y-ray 
activity of Zr® in counts/mg), So is the total Zr® deposited, D the diffusion 
coefficient, ¢ the annealing time, « the penetration distance, and A the specimen 
cross-sectional area. All of the runs carried out exhibited close linearity in a 
plot of log C vs. x*, a typical example being shown in Fig. 4 for a specimen 
annealed at 1210° C for 3 hours. 

The results are summarized in Table III, and a plot of log D vs. 1/T° K is 
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Fic. 4. A typical penetration vs. radioactivity plot. 
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TABLE III 


Summary of diffusion measurements 





Annealing 





Spec. No. temp. °K D in cm?/sec 
H 1503° C 4.2X10°* 
F 1410°C 2.5X10-° 
O 1362° C 2.2X10-* 
B 1362° C 2.1X10-° 
kK 1310° C 1.4X10°° 
G 1251 1.1X10* 
E 1251 1.2X10-° 
L 1210 7.3X10~° 
A 1210 7.3X10~ 
J 1168 6.1X10~ 
I 1168 5.7X10~ 


shown in Fig. 5 by the points in open circles. Dp and Q were determined by a 
least-squares fit to these data giving 


3 30,100 2 
: =? * exp— : : 
(3) D = 24X10 exp | RT | em /sec. 





It is apparent that here again D, lies well outside the usual range of values, and 
the activation energy Q is somewhat less than one-half the 70 kcal/mole 
anticipated on the basis of the melting-point correlation (pure zirconium melts 
at 2118° K). 

The results are compared with those of the Soviet workers by drawing the 
line of best fit to the present data through the assembled points of the former. 
See Fig. 6. In spite of some scatter, it would appear that the anomalous Dy 
and Q values are substantiated. 


4. EXPERIMENTAL PROCEDURE II 


Gruzin (1958) pointed out that all the zirconium specimens necessarily 
pass through the CPH a to BCC 8 phase change during the heating up to the 
annealing temperature. The resulting 6 grains could therefore contain either 
a high dislocation density or internal interfaces, either of which would provide 
short-circuit paths for the diffusion, and hence account for the observed Do 
and Q values. 

To check this possibility, the apparatus shown in Fig. 7 was constructed. 
An inactive zirconium specimen was placed in the small tungsten-wound 
furnace as indicated in the offset diagram. The furnace and specimen were 
brought up to the annealing temperature in the 6 phase and held there for 
3 to 1 hour, again in a vacuum of about 2X10-* mm Hg. 

This promoted the growth of large perfect 6 grains. (In some runs, as few 
as two or three grains were contained in disks of about 3 in. diameter.) The 
molybdenum radiation cover and zirconium specimen cover were then raised 
and swung out of the way, meanwhile maintaining the specimen and furnace 
at temperature. Finally, a layer of radioactive Zr® was evaporated from the 
tungsten filament onto the specimen, the cover replaced, and the diffusion 
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Fic. 7. The apparatus used for the measurements on samples in which no structural 
transformation took place prior to the diffusion anneal. 


anneal completed. Three runs at temperatures between 1280° C and 1100° C 
were thus made, and the diffusion coefficients determined as described in 
Experimental Procedure I. The results are plotted in Fig. 5 as triangles. 
Again the previous data are confirmed. 


5. DISCUSSION 


From the foregoing it seems safe to conclude that there exists a real dis- 
crepancy between the measured diffusion parameters in zirconium and similar 
BCC systems and those anticipated from empirical correlations or theoretical 
models (excluding that of Pound et al.). The question to be answered is 
whether this discrepancy provides a clue to a specific mechanism as Pound et al. 
suggest, or perhaps a modification of the vacancy model, or whether there are 
extraneous factors common to all of these measurements leading to consistent 
but nevertheless not truly representative results of volume diffusion. Such a 
common factor, for example, could be the high reactivity of Zr, U, Ti, and Cr, 
which makes it difficult to prepare specimens of a purity comparable to, say, 
the noble metals. In the present study, the use of a zirconium “gettering”’ 
boat and high-vacuum annealing was considered necessary since zirconium 
can absorb up to 25 atomic per cent oxygen substitutionally. The fact that 
the specimens were easily machined following an anneal, whereas the gettering 
boats were frequently quite hard, led us to believe the contamination of the 
former was kept very low. 
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As mentioned previously, Pound and his co-workers maintain that, while 
the Kirkendall experiments exclude the possibility of a ring mechanism operat- 
ing in BCC binary systems, this conclusion does not necessarily hold true for 
the pure metals. Fortunately, however, the sintering experiments described by 
Kuczynski (1950) provide a similar critical test of the ring mechanism for 
one-component systems. 

Kuczynski established that the transfer of mass during carefully controlled 
sintering experiments on silver wires can be accounted for in terms of volume 
diffusion, involving the flow of vacancies. The relationship between the 
diffusion coefficient D, the surface energy y, and the annealing time ¢ is given by 


5 
(4) D= 
wa é yt 
where x is the half width of the sintered neck between two wires of diameter a, 
6 is the lattice constant, 7 the absolute temperature, and k is Boltzmann’s 
constant. 

If the diffusion of the atoms takes place via a ring mechanism, no appreciable 
sintering should occur. Figure 8 shows a cross section of a bundle of 0.010 in. 
diameter zirconium wires wound on a zirconium thimble and sintered at 
1200° C for 22 hours. Assuming a value of 1000 dynes/cm? for the surface 
energy y, and averaging the neck widths over 100 measurements, one obtains 
a calculated value of D in order-of-magnitude agreement with the previous 





« - 





Fic. 8. <A cross section of a sintered bundle of zirconium wires, indicating the considerable 
degree of sintering occurring. 








1156 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


work, as seen in Fig. 5. While this is to some extent ambiguous in view of the 
assumed vy, the value for the latter is not likely to be more than +30% in 
error, as can be seen by comparing it with those of other similar metals for 
which y has been measured. 

This result, together with the observed Kirkendall shift in the BCC phase 
of U-Zr (Adda et al. 1959) lends strong support to the vacancy model. 

A further piece of evidence again pointing toward a vacancy mechanism is 
given by the observation of rapid pore formation in zirconium when bom- 
barded with monoenergetic a@ particles and annealed in the 8 phase. This 
effect, described in detail by Barnes (1959) and Ells and Evans (1959) for 
copper and aluminum, is attributed to the clustering of the injected helium 
atoms, with a subsequent migration of vacancies to the clusters, to form 
observable pores. 

Dr. Ells, who made the observations on the zirconium, pointed out that the 
pore formation was very rapid, consistent with the low activation energy of 
self-diffusion found in the present work. 


CONCLUSIONS 


On the basis of the foregoing, it appears that the ring mechanism must be 
excluded as an explanation of the observed results for the BCC phase of 
zirconium. While there exists the possibility of impurity effects lowering both 
Dy and Q, it is considered unlikely that they would be large enough to account 


for the observations. We conclude that the diffusion process involves vacant 
lattice sites as in the close-packed systems. If the parameters obtained in this 
study are representative of volume diffusion, it may be necessary to consider 
the detailed atomic configuration in the immediate vicinity of the vacant site 
in order to account for the departure of Do from the more usual range of values. 
In particular, the extent of relaxation of the neighboring atoms is generally 
thought to be considerably greater in BCC lattices than in, say, copper, and 
this could lead to a higher degree of co-operative motion being required to 
produce a jump. 

Until such time as a proper calculation of Q for these systems can be made, 
the ‘‘abnormality”’ of the activation energy cannot be assessed. It is felt that 
this latter aspect of the results is the more interesting of the two, and further 
work along these lines is looked for. 
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ANGULAR DISTRIBUTION OF FAST PHOTONEUTRONS! 


R. G. BAKER AND K. G. MCNEILL 











ABSTRACT 


The angular distributions and the yields of the high-energy neutrons emitted 
in photodisintegration have been studied by silicon detectors (25 elements) and 
aluminum detectors (6 elements). With the silicon detectors systematic variations 
are apparent in the coefficient a2 of W(@) = ao( Po+a2P2), and these variations 
are interpreted in terms of the Wilkinson shell model of photonuclear reactions. 









INTRODUCTION 


In recent years the Wilkinson picture of the nuclear photoprocess in the 
region of the ‘‘giant resonance”’ has achieved considerable success in explaining 
many of the phenomena connected with this process (Wilkinson 1956, 1959). 

In this picture, the absorption process is associated primarily with the 
excitation of levels which correspond to single-particle excitations. The dis- 
integration process results from the decay of the excited single-particle state 
which proceeds either by direct emission of the particle, or by reabsorption 
of the excited particle into the nucleus and eventual de-excitation of the system 
by allowable competing processes. In principle, the former process is much 
the less complicated of the two, and definite correlations may be expected, 
say, between the angular distribution of the emitted nucleon and the initial 
single-particle state of the nucleus—photon system. On the other hand, the 
great internucleon interactions of the latter process, described by a statistical 
or evaporation model, smears out any definite angular correlations and the 
disintegration products would be expected to have, and do have at least for 
neutrons, an isotropic distribution. One also expects a differentiation of the 
decay mode in terms of the energy distribution of the emitted nucleons. Directly 
emitted nucleons leave the residual nucleus at or near the ground state and 
carry all or most of the available energy, while those emitted by the statistical 























process have a typical evaporation spectrum. 

Focussing one’s attention on the photoneutron cross section in the region 
of the maximum of the ‘‘giant resonance” at a gamma-ray energy E,,, and 
assuming that most of the single-particle excitations occur for those neutrons 
with the least binding energy (£,), one expects the directly emitted neutrons 
to have an energy in the neighborhood of E,,—E, ~ 6-7 Mev, while those 
emitted by the evaporation process will have a Maxwellian distribution charac- 
terized by the nuclear temperature. For a typical nuclear temperature of 
0.75 Mev, characteristic of the heavy elements and independent of the excita- 
tion energy (Gugelot 1951; Graves and Rosen 1953), less than 1% of the 
evaporation neutrons would have energies greater than 5 Mev, and in the 
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neighborhood of 10% would have energies greater than 3 Mev. Thus, a 
simple threshold detector should serve to differentiate the direct interaction 
spectrum from the evaporation spectrum. 

While there has been a considerable amount of work already carried out in 
this field, the results of different workers are not in agreement. The disagree- 
ment may be attributed largely to an ill-defined detector threshold and 
variations in spectral sensitivity, both of which vary widely with the experi- 
mental technique. However, the experimental results obtained for the angular 
distribution can always be expressed in terms of an expansion in Legendre 
polynomials, 


W(0) = aod) (1+a,P,(cos 8)) 
n=1 
where only the coefficient a2 makes an appreciable contribution, indicating a 
process consistent with a mainly dipole type of interaction, although some- 
times the possible addition of a significant a; term arises from dipole—quadru- 
pole interference. Table I contains some of these results, with notes as to 











TABLE I 
(En)min 
Element —as Ey (Mev) (Mev) Method Reference 

Pb <0.07 22 ZnS Geller et al. (1954) 
0.40 65 5 ZnS Johansson (1955) 
0.05+0.05 70 ZnS Dixon (1955) 
0.17+0.06 17 ZnS Asada et al. (1958) 
0.36+0.04 22 3 Al(n,p) Price (1954) 

0.13 23 4 Emulsion Toms and Stephens (1957) 

Bi 0.40+0.04 22 3 Al(n,p) Price (1954) 
0.31+40.05 30 3 Al(n,p) Ferrero et al. (1956) 
0.38+0.05 30 5 Si(7,p) Ferrero et al. (1956) 
0.57 19 4 Emulsion Zatsepina et al. (1957) 

Ag <0.07 22 ZnS Geller et al. (1954) 
0.18+0.06 22 3 Al(n,p) Price (1954) 

Cu <0.07 pt ZnS Geller e¢ al. (1954) 
0.13+0.09 70 ZnS Dixon (1955) 
0.19+0.04 22 Al(n,p) Price (1954) 

Ta 0.37 65 5 ZnS Johansson (1955) 

Al 0.43 65 5 ZnS Johansson (1955) 
0.19+0.15 70 ZnS Dixon (1955) 

Fe 0.06+0.17 70 ZnS Dixon (1955) 


the method used in detection of the direct (fast) photoneutrons, the maximum 
photon energy of the bremsstrahlung spectrum, and an estimate of the detector 
threshold. 

The present paper reports results of measurements made using silicon and 
aluminum threshold detectors on a series of 25 elements in the medium- and 
heavy-weight parts of the periodic table. Some preliminary results of this 
work have been reported previously (Baker and McNeill 1960). 
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EXPERIMENTAL TECHNIQUE 


In the present work, use has been made of Si?* and Al?? threshold detectors. 
For both these nuclides, the (”,p) reaction is endoergic requiring neutron 
energies of greater than 3.9 Mev and 1.8 Mev respectively, and for both 
elements the end products of the reactions are beta radioactive with half-lives 
of 2.3 minutes and 9.5 minutes respectively. Owing to the small cross section 
near threshold, effectively the silicon detector responds only to neutrons of 
energy greater than 5 Mev, while the aluminum detector responds to those 
with energies greater than 3 Mev. Thus, the response of these detectors could 
be expected to show the characteristics of the ‘‘direct interaction’’ spectrum 
provided it constituted an appreciable fraction of the total spectrum, e.g. 
more than 1% of the total for the case of the silicon detector. 

The aluminum detector is also activated by the (n,y) reaction which is 
exoergic, and has a cross section with a well-defined £,~!”* dependence. The 
reaction product, Al*’, is beta active with a half-life of 2.3 minutes. This 
activity will reveal the characteristics of the low-energy or evaporation 
spectrum, while in the same detector the Al?’ (m,p) Mg?’ reaction, with a 
half-life of 9.5 minutes, shows the characteristics of the direct emission spec- 
trum. Other possible competing reactions, such as (y,m) or (n,a), which all 
have products with half-lives of the order of a few seconds or many hours, are 
easily eliminated or dealt with in the analysis of the activation data. 

In the silicon detector, possible interfering activities are the (7,p) reactions 
in Si?® and Si® giving rise to half-lives of 2.3 minutes and 6.6 minutes, and 
the (”,a) reaction in Si*® giving rise to a 9.5-minute activity. Compton scattered 
photons from the target would have energies below the (y,p) threshold for all 
practical angles of observation and the 9.5-minute Mg?’ produced by Si®°(,a) 
could in principle be distinguished by its half-life. However, in practice a pure 
2.3-minute activity attributed to Si?*(7,p)Al?* was found in all cases for the 
silicon detector. 

Each detector was in the form of a cylindrical annulus of inside diameter 
3 in. designed to slip over a Victoreen 1B85 Geiger counter in a lead castle. 
The aluminum detectors had 3-in. walls and were 6 cm in height. The silicon 
detectors were formed from powdered silicon packed in an annulus formed 
of 0.002-in. Mylar film to give a }-in. wall and 7-cm height. 

A diagram of the experimental arrangement is shown in Fig. 1. The brems- 
strahlung beam from the Allis-Chalmers Betatron at the Ontario Cancer 
Institute was collimated by lead and iron. Extraneous neutrons were moderated 
and absorbed by paraffin and cadmium. A jig held the target samples in the 
center of the beam at 120 centimeters from the bremsstrahlung source. Detec- 
tors could be placed in shielded positions at 10 centimeters from the samples 
and at 30°, 60°, 90°, 120°, and 150° to the direction of the incident beam. 

“No sample background”, due mainly to fast neutrons generated in the 
shielding, was always appreciable, and was sampled before, during, and after 
a series of runs. 

Throughout the experiments, the betatron energy was maintained within 
limits of +0.1 Mev of a calibrated value of 22 Mev. The absolute magnitude 
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SCHEMATIC EXPERIMENTAL ARRANGEMENT 
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Fic. 1. Experimental arrangement for the detection of the fast photoneutrons, showing the 
shielding and collimation. 


of the betatron energy was determined by measuring the peak voltage induced 
in a leg turn of the magnet at the thresholds for the (y,m) reactions in Cu® 
and C” taken as 10.73 and 18.73 Mev respectively (Spicer 1955; Tobin et 
al. 1958). Since the relation between peak voltage and peak bremsstrahlung 
energy is linear (McElhinney et al. 1949), the two points suffice for calibration. 

The y-ray flux from the betatron was monitored by a built-in ionization 
chamber situated in the direct beam immediately in front of the collimator. 
This chamber was calibrated against a standard Victoreen chamber with a 
5-cm lucite wall placed at the sample position. Corrections were made for 
the beam attenuation in the sample and for the angular variation of intensity 
of the beam. 

The center of gravity for the production of neutrons, taking into account 
the attenuation of the beam in the target samples, was calculated and the 
samples so placed that this point was at the center of the angular distribution 
jig. For cylindrical lead and bismuth samples of 3 cm diameter, the displace- 
ment from the center of mass was 3.5 mm. 

The finite sample size and consequent self-scattering of neutrons will tend 
to reduce the anisotropy produced by the nuclear reactions, but as this effect 
has been calculated to produce less than a 10% attenuation in the distribution 
coefficients for the worst case, no further account has been taken of it. 

During a run, the samples were irradiated for a given duration at a constant 
exposure dose rate, and then the detectors removed immediately to a counting 
room. The counting data together with the irradiation data give the relative 
activation due to Al(m,p) and Al(n,y) reactions or the Si(m,p) reaction as a 
function of target material and angle of observation. After normalization to a 
standard sample and exposure of 100 r/minute mole, the saturated yields for 
each element were fitted to the Legendre polynomial Y(@) = ao{1+a:P:(cos 6) 
+a2P.(cos 6)|. The average observed yield factor do is related to the average 
absolute neutron flux produced by the sample by an efficiency factor dependent 
on the neutron reaction cross section and on a geometry factor. The detection 
efficiency factors for the two types of detectors were found using a standard 





1162 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


Ra-Be neutron source obtained from Atomic Energy of Canada Limited, by 
considering the neutron spectrum from this source (Hill 1952; Teucher 1949) 
and the Si(u,p) and Al(n,p) cross sections (Neutron Cross-Sections, USAEC- 
BNL 325). The calibration factor so obtained allowed the data to be presented 
in terms of 
fo(n,p)d(En)dE = 6(n,p)® 

where o(,p) is the cross section of the detector nucleus in millibarns at neutron 
energy E,, 

¢(E,) is the number of neutrons in the energy interval dE, at E, throughout 
the total solid angle per 100 r per mole of sample, 

é(n,p) is the average cross section of the detector nucleus weighted in accor- 
dance with a given neutron energy spectrum, and 

® is the total number of neutrons per 100 r per mole of sample. 

The Al(n,y) detection efficiency could not be estimated to better than an 
order of magnitude since the combination of the short half-life of the reaction 
product and the absence of an appreciable low-energy neutron flux from the 
Ra-—Be source gave an activation which was not statistically significant in the 
presence of the Al(m,p) activation. 

The yields for six elements using the aluminum detectors are tabulated in 
Table II. 

The two sets of results shown for bismuth were from runs about two months 
apart and are indicative of the reliability and reproducibility of the results. 

While the statistics for the detection of the low-energy portion of the 
neutron spectrum are poor, and one would expect isotropy of emission, there 
does appear to be a suggestion of greater emission at the larger angles. How- 
ever, without further evidence, it would be unwarranted to conclude that 
this is a real effect. Accordingly only the value of ado is shown in the table. 

The angular distributions for the fast neutron component, detected by the 
Al(n,p) reaction, show considerable anisotropy for the neutron closed shell 
nuclei La, Pb, and Bi and considerably less for nuclei removed from closed 
shells. 

Published photoneutron spectra for Pb (Toms and Stephens 1957; Dixon 
1955), Ta (Cortini e¢ al. 1958), and Cu (Dixon 1955) were used to calculate 
& for the Al(n,p) reaction. Estimates of @ at 22 Mev were taken from pub- 
lished data (Montalbetti et al. 1953; Goldemberg and Katz 1954). Table III 
gives the comparison of ¢(Al(,p)) for this work and that derived from the 
photoneutron spectra. The agreement is good, considering that the calculations 
derive from material with large and uncertain correction factors and large 
statistical errors, and therefore gives confidence in the calculations based on 
the Ra-—Be calibration. 

The aluminum detectors have the advantage, in principle, of differentiating 
between the fast direct emission component and the compound nucleus 
emission component, and of giving information relative to the two processes 
simultaneously. Unfortunately, the counting statistics are subject to large 
standard deviations so that the angular distribution coefficients are not well 
defined. 
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TABLE III 


o Al(z,p), millibarns 
@ at 22 Mev ——— 





Element _ per mole/100 r This work From spectra Spectrum reference 
Pb 2.72X10° 6.87 6.47 Toms and Stephens (1957) 
5.64 Dixon (1955) 
Ta 2.50109 4.49 6.1 Cortini et al. (1958) 
9.50 9.02 Dixon (1955) 


Cu 0.30 10° 





With the silicon detectors, the statistics were somewhat better, and the 
analysis was not subject to multiple regression errors since the data could 
be treated as merely a 2.3-minute activity associated with a background. The 
results are given in Table IV and are also shown graphically in Figs. 2 and 3. 

Figure 2 portrays the data of Table IV, column II, showing the variation 
of ap with atomic number and also with neutron number. The values of a, 
as shown in column III, do not show any unambiguous trend with the excep- 
tion that the over-all trend is to a small positive coefficient. The next section 
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Fic. 2. The activity of the silicon detectors (proportional to the yield of fast photoneutrons ) 
as a function of atomic number Z and neutron number N. 
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COEFFICIENT (-a,) of W(O)=R ta, PR + ae Pp 


Si(n,p) Activation by Photoneutrons 


Eo = 22 Mev 


30 40 50 60 70 80 90 
Z 


Fic. 3. The variation of the coefficient a2 in the angular distribution W(@) = Po+a,;P,+a2P2 
as a function of atomic number and neutron number. 


contains the discussion of the trend and interpretation of the a2 coefficients 
of column IV which are portrayed in Fig. 3 as a function of atomic number 
and neutron number. The values of (¢ Si(m,p) ®) of column V are simply 
related to the a» of column II by a detector efficiency factor obtained from 
the Ra-Be calibration. For our particular experimental arrangement, 
(¢@) = 2.47X107 a> millibarn-neutron. Column VI gives the number of 
neutrons per mole per 100 r at 22 Mev as given by Montalbetti et al. (1953) 
and Goldemberg and Katz (1954). An estimate of the ratio of the fast, direct 
neutron flux to the total neutron flux, as shown in column VII, was made 
by assuming that the fast neutrons have a Gaussian energy distribution 
about E, = 7 Mev with a half-width of 4 Mev for which ¢ Si(n,p) = 245 
millibarns. It is noted that from 6% to 15% of the photoneutrons may be 
interpreted as fast, direct interaction products. These rough figures are in 
agreement with the results of Ferrero et al. (1956) for Bi, and of Poss (1950) 
for Bi, Pb, and W. 

The data of columns V and VI were used to calculate ¢, which is compared 
to the value of ¢ obtained by use of published spectra. The comparison is 
made in Table V. The agreement is satisfactory considering the statistical 


nature of the data. 
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TABLE V 








¢ Si(n,p), millibarns 
Element Thiswork From spectrum Spectrum reference 


Pb 20.8 16.25 Toms and Stephens (1957) 
17.55 Dixon (1955) 
15.4 23.2 Cortini et al. (1958) 
28.6 33.5 Dixon (1955) 
23.5 17.0 Cortini et al, (1958) 








DISCUSSION OF THE ANISOTROPY COEFFICIENT, a2 


It will be seen from Fig. 3 that although there are obvious magic number 
effects, they are not obviously straightforward effects; thus, whereas |a2| is 
a maximum for NV = 82 and 126, it is a minimum for NV = 28 and 50. 

In the Wilkinson picture, a neutron in a shell model level may be raised 
by an £1 transition to a higher unoccupied level from which it has probability 
of escaping directly into the laboratory before colliding and sharing its energy 
with other nucleons. The relative probability of a particular / — /’ transition 
for an individual neutron is given by (Courant 1951) 


tet for l’ = 1+1 


2/+1 


l a 
341? for l’ = 1-1 


where 


2 


D= J Ruy Ra r'dr 
0 


The radial overlap integrals are expected to be larger for transitions in- 
volving no change in the number of radial nodes (i.e. 1 > /+1) than in those 
involving such a change (i.e. /—+/—1) by up to an order of magnitude. The 
probability of direct escape is dependent on /’ and on the energy of the photon 
responsible for the transition minus the binding energy of the affected neutron. 
Values of the direct escape probability, 7), were taken as the transmission 
factor for neutrons in a square well potential (Blatt and Weisskopf 1952, 
Ch. VIII, Eq. (5.7)). While this is a crude approximation, it is adequate for 
a first-order qualitative interpretation. Other estimates, based, for example, on 
the empirical optical model transmission factors (Beyster et al. 1957) or on the 
overlap integral evaluated outside the nuclear radius (De Sabbata and Toma- 
sini 1959), exhibit the same trend as the transmission probabilities for the 
square well for which approximate values for neutron energies of the order 
of 7 Mev are given below in Table VI. 

Table VI assumes most of the transitions occur from the same oscillator 
level. Thus, 1g, 2d, and 3s neutrons occupy oscillator level number V and £1 
transitions 1g-2f and 2d-2f have the same final state and the same escape 
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TABLE VI 
Approximate values of Ty 





Transition Tr 
1] > 1(1+1) 0.0-0.: 


21 > 2(1+1) , 
11 + 21-1) oe: 


31 — 3(1+1) 
21 — 3(1—1) 





probability. The main point to note is that for the same initial orbital angular 
momentum, the escape probability favors the 1] + /—1 transition. 
The angular distribution terms for the direct interaction are (Courant 


1951) 


}-+-2 2 
W(@) = ad Po(cos 6) aes P»2(cos 0) | for /—1+1, 
l-1 
2(2/+1) 


The strength functions and angular distribution functions with and without 
spin-orbit coupling considerations are also given in terms of the Racah algebra 
coefficients by Fujii (1959). 

It is necessary to have a knowledge of the ordering of the nuclear shells 
and subshells before an interpretation of the experimental results is possible. 
In the following discussion, an ordering based on experimental evidence is 
used as displayed in Fig. 4. The references to the work, pertinent to a par- 
ticular region of the periodic table, are given in the caption of the figure. 

Using the above ideas, and the fact that the silicon detectors will not 
observe neutrons of less than 5 Mev, it is possible to explain the variation 
of the angular distribution coefficient a2 with neutron occupation number. The 
main things to note are that (1) /—>/+1 transitions give —a2 > 0.25, while 
1— 1-1 transitions give —a, < 0.25; (2) / -+/—1 transitions are less prob- 
able than the / +/+1 type; and (3) on the other hand, / > /+1 transitions 
resulting in high /’ will be inhibited by the angular momentum barrier from 
appearing as direct interactions. 

The increasing anisotropy from N ~ 28 to N ~ 36, in which region the 
2p3/2 and fs;2 subshells are filling, may be attributed to the 2p — 2d tran- 
sitions with a small admixture of 2p — 3s. Considering the yields (Fig. 2) 
from the Z odd and Z even nuclei separately, they are noted to increase 
as the 2p shell fills. The alternation of the Z odd and Z even yields may be 
associated with the (y,”) threshold energy which is higher in the Z even 
nuclei and results in less energy being available to the directly emitted neutron 
and thus, a lower barrier penetration probability and detection probability. 
From N ~ 35 to N = 50, decreased anisotropy reflects the angular momentum 
barrier effect in suppressing the g9/2 > Aiij2 and fs;2 — gz/2 direct emission 
products and favoring the 1/ — 2(/—1) transitions which have —a2 < 0.25. 
The interval from N = 50 to N = 76 marks the filling of the 2d5/2, gz)2, and 


P.»(cos 0) | for 1 — 1-1. 


a,| Po(cos 6) 
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COMPOSITE NEUTRON LEVEL DIAGRAM 


(based on experimentol evidence) 


Neutron Reference 
Energy| Arrangement Number 
(Level) (Total) 
4 





Fic. 4. The assumed ordering of nuclear levels used in the discussion. The references cited 
in the figure are: (a) Holm e¢ al. (1960); (b) Blomquist and Wahlborn (1960); (c) Kurath 
(1960); (d) Zeldes (1958); (e) Lawson (1960); (f) Cohen and Price (1960). 
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hii /2 subshells, with the 2d5 2 — 2f7/2 transitions producing greater anisotropy. 
Other transitions from these subshells are rather ineffective due to either 
small transition strengths or suppression by the angular momentum barrier. 
The marked increase in anisotropy at N = 82 is presumably due to the 
3s — 3p and 23/2 — 2f5,2 transitions. Beyond N = 82, subshells of higher 
angular momentum are filling, and only / — (J—1) transitions give an appre- 
ciable contribution, supplemented by transitions from the N = 82 core. The 
yield from the core is reduced as the subshells fill, and the yield from the 
| — (J—1) transitions no more than compensate for this reduction yieldwise, 
but in so doing gives a decreased anisotropy in the angular distribution. At 
N = 126, the 3p — 3d transitions could account for the greater anisotropy 
if one accepts that the 3p — 4s transition strength is very weak, or that the 
4s level is depressed so that direct emission from it is at too low an energy 
to be detected efficiently. Thorium (25% photofission) and uranium (50% 
photofission) (Price and Kerst 1950) show the effect of isotropic fast photo- 
fission neutrons in their more isotropic distribution. 


CONCLUSION 


Comparing the results obtained here with those of other workers, it is 
noted that the only previous study of a variety of elements with a similar 
technique was by Price (1954), some of whose results were quoted earlier 
(Table I). The elements chosen by him happen to be such that they indicated 


an increase in anisotropy with increasing neutron number. While in essential 
agreement with the work reported herein, the significant variations of the 
angular distribution coefficient exhibited in Fig. 3 could not be anticipated. 

The only report appearing in the literature attempting a theoretical corre- 
lation of experimental yield and angular distribution of fast photoneutrons 
from an element in the range of those studied here is that of Ferrero e al. 
(1956). Their results, quoted earlier, give —a2 = 0.38+0.05 for bismuth for 
30 Mev bremsstrahlung, which is only slightly different from the present 
result of 0.45+0.04 for 22 Mev bremsstrahlung. They consider the total 
transition strength to be from allowed transitions of the shell model type 
nl — n(l+1), and ignore ni — (n+1)(/—1) transitions. Assuming an average 
energy of 7 Mev for the directly emitted neutrons and applying an appropriate 
transmission factor, they obtained good agreement with experiment. Although 
ignoring the /— (/—1) transitions does work in the region of N = 126, it 
fails for many other elements; indeed one could never get —a2 < 0.25 if the 
/— (/—1) transition strengths are neglected. 

From the correlation found here between the predictions of the Wilkinson 
model of photodisintegration and the experimental results, it is clear that the 
present results give strong evidence in favor of the independent particle model 
for the direct effect, at least as applied to the angular distribution phenomena. 
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CROSS SECTIONS OF (p,pxn) REACTIONS IN Au! 


T. M. KAVANAGH AND R. E. BELL 


ABSTRACT 


Cross sections of (p,pn), (p,p2n), and (p,p3n) reactions in Au!” have been 
measured by the activation method for incident proton energies up to 86 Mev. 
The curves of cross section as a function of energy have similar shapes for the 
three reactions. They rise from apparent thresholds at about 16, 21, and 30 Mev, 
respectively, to peak values of 180, 145, and 150 mb at proton energies about 
30 Mev higher than the apparent threshold energies. The cross sections are much 
larger than those predicted from the statistical and cascade-evaporation theories, 
and they are interpreted in terms of two-body collisions in the diffuse surface of 
the target nucleus. A combination of these results with measured (p,xm) cross 
sections yields an approximation to the total reaction cross section of a heavy 
nucleus. 


INTRODUCTION 


In an earlier paper from this laboratory Bell and Skarsgard (1956) reported 
the results of measurements of (p,xm) cross sections of heavy elements. This 
work has now been extended to (p,pxn) reactions. Cross sections of (p,pn) 
(p,p2n), and (p,p3n) reactions in Au’ have been measured in the proton 
energy range 18 to 86 Mev; the experimental method was similar to that 
described in the earlier paper. 

Reactions of the type (p,pxn), which involve inelastic scattering of the 
incident proton, are inhibited by the Coulomb barrier in both the entrance 
and exit channels. On the basis of statistical theory such reactions in a heavy 
target element are expected to have negligibly small cross sections compared 
with (p,xn) reactions. Bell and Skarsgard found that near the peak of a 
particular (p,xm) curve, the (p,p(x—1)m) cross section was indeed small by 
comparison, and they were led to expect that the sum of the (p,xm) cross 
sections would approximate the total reaction cross section. This sum curve, 
however, showed structure that could not be explained on the basis of 
statistical theory (Jackson 1956). 

It is now well known that inelastic scattering of medium-energy protons 
cannot be described in terms of compound nucleus processes (Eisberg and Igo 
1954; Cohen 1957). For example, Eisberg and Igo found that for 31-Mev 
protons incident on heavy nuclei, inelastic scattering accounts for about 15% 
of the total reaction cross section; furthermore the scattered protons have 
distributions in energy and angle that are characteristic of direct interactions. 
While the interpretation of these medium-energy direct interactions has been 
the subject of some controversy, they probably result from collisions of the 
incident protons with nucleons in the surface of the target nucleus (Griffiths 
and Eisberg 1959; Eisberg 1960). 

The present experiments constitute a study of cross sections of some in- 
elastic scattering reactions. Cross sections have been measured as a function 
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of proton bombarding energy over a range in which the surface mechanism is 
expected to predominate. The use of a heavy target nucleus makes compound 
nucleus contributions unimportant. At the higher energies internal direct 
interactions can lead to the emission of protons, but our experimental cross 
sections exceed those predicted by cascade-evaporation theory. The results 
indicate that inelastic scattering can account for the structure in the Bell and 
Skarsgard (p,xn) sum curve, and an addition of (p,pxn) cross sections to this 
curve gives an approximation to the total reaction cross section of a heavy 
nucleus. 

Previous measurements of (p,pxm) cross sections of heavy elements are not 
very comprehensive. Cohen et al. (1954) have measured the (p,pm) cross 
sections of silver, palladium, and tantalum for proton energies up to 23.5 Mev. 
The (p,pm) cross section of gold has been measured by Sosniak (1958) for 
energies up to 32 Mev, and by Yule and Turkevich (1960) at 82 Mev and at 
several higher energies. Gusakow et al. (1960) have published (p,p”) and 
(p,p2n) cross sections of gold for proton energies in the range 40 to 155 Mev, 
with results in strong disagreement with our own. The (p,pxn) cross sections 
of cesium for x = 1, 2, 3, 4, and 6 have been reported by Fink and Wiig (1954) 
for proton energies of 60, 80, 100, 150, and 240 Mev. 

Measurements made on medium- and light-weight elements are more 
extensive, but of less interest here. The C"(p,pn)C"™ reaction was used, 
however, for monitoring the proton beam. 

It should be pointed out that in the present work cross sections have been 
measured by the activation method. This method does not distinguish whether 
or not the outgoing particles are bound, and thus, for example, the measured 
(p,pn) cross section includes the (p,d) cross section. 


EXPERIMENTAL METHOD 


The experimental method was essentially that described by Bell and 
Skarsgard (1956). In each bombardment a gold target foil, } in. wide and 
either 20 or 40 mg/cm? thick, was mounted together with a graphite beam 
monitoring foil in the circulating beam of the McGill synchrocyclotron. The 
product nuclides formed by the various (p,pxn) reactions decay primarily by 
orbital electron capture, and cross sections were determined by absolute 
counting of K X rays and gamma rays that follow decay events. 

The monitor foils were sheets of pure carbon 20 mg/cm? thick, prepared 
with a surface grinder from fine-grain compressed graphite. These foils had an 
advantage over the teflon monitors used by Bell and Skarsgard in that they 
could withstand the heating caused by long bombardments. Beam monitoring 
was based on the reaction C"(p,pm)C"™, and the monitor cross section used 
was that given by Aamodt et al. (1952) modified according to the recommenda- 
tions of Whitehead and Foster (1958). Uncertainties in the energy of the 
proton beam limit the accuracy of carbon monitoring below about 30 Mev, 
since in this region the monitor cross section is a strong function of energy. 
Only the first few points on the (p,pm) curve are seriously affected, however, 
and these points are also in a region of rapidly varying cross section; thus errors 
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in monitoring have little effect on the appearance of the curve. Sosniak (1958) 
gives (p,pn) points at 18 and 21 Mev, monitored by the F"(p,pm)F8 reaction 
which is less energy dependent in this region, and these data define the (p,p7) 
cross section at energies where the present results are uncertain. 

The large amounts of radioactive mercury formed in the gold foils by 
(p,xn) reactions were removed by vacuum distillation. Two different techniques 
were used. A slow method involved heating the gold foils after bombardment 
in an evacuated quartz tube; the mercury was driven into a liquid air trap, 
and at temperatures of about 800° C, separation was complete in 15 to 20 
minutes. This procedure was used in initial studies of the (p,pm) and (p,p3n) 
reactions for bombarding energies up to 60 Mev. 

A more elaborate technique was required for the (p,p2m) measurements 
since the product nuclide Au™ is also formed in the target foil during a slow 
separation by the decay of relatively short-lived Hg’. Corrections for this 
mercury decay are difficult since the cross section for formation of Hg™ is 
not known, and more important, since an unknown amount of mercury is lost 
into the cyclotron vacuum during the bombardment due to heating by the 
proton beam. A fast separation technique was developed which involved 
heating the targets electrically in the cyclotron during bombardment. Gold 
targets 20 mg/cm? thick (half the thickness used for the slow method) were 
heated to 800° to 900° C by passing about 5 amperes of 100-kilocycle alter- 
nating current through them. The high frequency was necessary in order to 
prevent destruction of the foils in the 16-kilogauss field of the cyclotron. The 
target strips had a 0.01 in. wide slit cut longitudinally along the center of the 
strip from the base to within § in. of the tip, so that both electrical terminals 
were at the base of the target, and the target tip could be cut to fit the monitor 
foil. The target’s temperature was indicated by its resistance, and was 
controlled remotely during bombardment. 

A tracer experiment was performed as a check on the above technique. A 
gold target of the same form as those used in the cross-section measurements 
was bombarded at 28 Mev, where most of the resultant radioactivity was due 
to mercury. Low beam current was used to prevent mercury boil-off. This 
sample was heated electrically in an evacuated glass vessel provided with a 
liquid air trap, and K X rays originating in the sample were counted with a 
sodium iodide counter. The counting rates before and after the heating current 
was turned on were recorded on a chart recorder, and the results are shown in 
Fig. 1. The curve indicates a ‘“‘mean life’’ against boil-off for a mercury atom 
in the foil of about 40 seconds. This time is very short compared with the 
radioactive half-life of Hg’, 9.5 hours. 

A further check on the effectiveness of the fast separation was offered by a 
sample bombarded at about 23 Mev. The measured (p,p2m) cross section 
(8.1 mb) is comparable with the value of 9.5 mb reported for the (p,t) reaction 
at this energy (Cohen and Rubin 1959). The (p,3”) reaction leading to Hg™ 
probably has a cross section of about 400 millibarns at 23 Mev (Bell and 
Skarsgard 1956), and a decay of only 1% of the Hg™ before its removal from 
the foil would result in an apparent (p,p2n) cross section at least 30% greater 
than the (p,t) value reported by Cohen and Rubin. It is thus concluded that 
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Fic. 1. Results of a radioactive tracer experiment performed as a check on the fast mercury 
removal technique (see text). The curve indicates the rate of vacuum distillation of mercury 
from a 0.0005 in. thick gold foil heated to about 800° C. 


more than 99% of the mercury formed in (p,3”) reactions was removed before 
it decayed into Au™, Finally, several targets were maintained at an appreciably 
lower temperature (~600 degrees) during bombardment, with no significant 
effect on the measured (p,p2m) cross section. 

The fast mercury removal also facilitated the measurement of (p,pm) and 
(p,p3n) cross sections above 60 Mev, since there was significantly less inter- 
fering activity produced in the samples by decay of the mercury isotopes 
formed at high energies. 

X Rays and gamma rays from the active samples were counted with sodium 
iodide scintillation counters in conjunction with a multichannel pulse analyzer. 
The (p,pn) and (p,p3n) reaction products (Au™ and Au™) were counted with 
a single calibrated crystal; the method of determining reaction yields from 
counting rates for particular total energy events has been described by Bell 
and Skarsgard. Au™ activity from (p,p2m”) reactions was counted by sand- 
wiching the sample foils tightly between two sodium iodide counters connected 
in parallel. The counter gains were made equal by adjusting the high voltages, 
and the assembly functioned as a single counter with about 95% geometrical 
efficiency. Disintegration rates were calculated from total counting rates of 
all detected events. This method had the advantage of being extremely 
insensitive to uncertainties in the Au™ decay scheme, and also enhanced the 
observed counting rates (Au has a half-life of about 6 months and thus the 
disintegration rates were low). 


EXPERIMENTAL MEASUREMENTS AND RESULTS 


(p,pn) Cross Section 

Yields of the product nucleus Au were determined from the total photopeak 
counting rates of the 331-kev and 354-kev gamma rays that follow electron 
capture, and the 426-kev gamma rays that follow B- decay. The branching 
ratios and conversion coefficients used were those given by Thieme and 
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Bleuler (1956a). Cross sections were calculated from counting rates after the 
decay of 9.5-hour Au”, and no correction was made for reactions that 
proceed via the isomeric state; the error introduced is negligibly small. 
Decay curves of Au’ activity were followed for about 25 days. The measured 
half-lives were in good agreement with the value of 6.1 days reported by Yule 
and Turkevich (1960), and this value was used in the cross-section calculations. 
The (p,pn) results are summarized in Table I, and are shown graphically 


TABLE I 


Experimental Au’ (p,pxn) cross-section results 


Proton (p,pxn) cross sections (mb) 
energy 


(Mev) (bpm) —_—_(P,P2n)—_(,p3n)_ 


18.0 23.6 
22.3 83.3 
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by the open circles in Fig. 2. The solid curve was simply drawn freehand 
through the experimental points, and does not represent any theory. The 
standard error in the points, compounded from uncertainties in the individual 
quantities involved in the cross-section calculations, is about 13%; this is 
represented by an error bar in the figure. The crosses show (p,pm) cross sections 
measured by Sosniak (1958), and the triangle is the lowest-energy point from 
high-energy results of Yule and Turkevich (1960). The dashed curve, which 
disagrees strongly with the present work, is due to Gusakow et a/. (1960). 


(p,p2n) Cross Section 

Cross sections were calculated from the total counting rates recorded by the 
high-efficiency counter described above, that is, from integrated pulse height 
spectra. The pulse spectrum showed a strong K X-ray peak, a 99-kev gamma- 
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Fic. 2. Cross section of the reaction Au'?(p,pn)Au® as a function of incident proton 
energy. The circles are from the present work, and the error bar represents the estimated 
standard errors in the absolute cross-section values. The crosses are from Sosniak (1958), the 
“oe from Yule and Turkevich (1960), and the dashed curve was published by Gusakow 
et al. (1960). 


ray peak, a strong peak at about 130 kev due primarily to K+K sum events, 
and weaker contributions due to summing of K X-ray pulses with 99- and 
130-kev gamma-ray pulses. The decay scheme parameters used were not 


critical; calculations were based on values published by Bisi et al. (1959). 

The long-lived Au™ activity was counted for a period of from 7 to 9 months 
for the different samples, beginning 3 months after bombardment. The decay 
curves all showed a half-life of 1885 days, and this value was used in the 
cross-section calculations. It is in good agreement with the values of 185+3 
days reported by Wilkinson (1949), and 192+5 days given by Bisi et al. 

The (p,p2n) cross-section results are summarized in Table I and are shown 
graphically in Fig. 3. The standard error in the values reported is about 13%. 
The dashed curve in the figure is due to Gusakow et al. (1960). 


(p,p3n) Cross Section 

This reaction leads to the formation of 39-hour Au, and cross-section 
determinations were based on counting rates of K X rays, and 291- and 327-kev 
gamma rays. The gamma rays were unresolved from Au™ gamma rays, and 
their intensity, as well as that of K X rays from Au™, was determined from 
complex decay curves. Decay scheme parameters used were due to Thieme 
and Bleuler (19560). 

At lower energies (below 45 Mev) absolute cross sections were calculated 
from K X-ray counting rates. At higher energies, however, there is inter- 
ference from other activities, particularly 17.5-hour Au; the gamma rays were 
used as an indicator of relative cross sections over the whole energy range, and 
these relative cross sections were normalized to the absolute values obtained 
by K X-ray counting at low energies. 

The (p,p3m) results are shown in Table I and Fig. 4. The standard error in 
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Fic. 3. Cross section of the reaction Au'7(p,p2n)Au® as a function of incident proton 
energy. The error bar indicates the estimated standard error in the absolute cross-section 
values. The dashed curve is from Gusakow ef al. (1960). 

Fic. 4. Cross section of the reaction Au?(p,p3”)Au™ as a function of incident proton 
energy. The error bar indicates the estimated standard error in the absolute cross-section 


values. 


the cross-section values is somewhat larger than for the (p,pn) and (p,p2n) 
reactions since yields were calculated from complex decay curves; the error 
is estimated to be about 16%. 


DISCUSSION 
(p,pxn) Cross Sections 
The (p,pn) and (p,p2n) cross-section results from the present work are in 
strong disagreement with the values reported by Gusakow et al. (1960) (see 
Figs. 2 and 3). Their values are, for the most part, lower than the present ones 
in the energy range where the two investigations overlap. It should be noted 
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- that Gusakow et al. also show a (p,3m) cross-section curve which does not 
resemble (p,3m) curves for other heavy elements (Bell and Skarsgard 1956). 
The (p,p7) results from the present work are in very good agreement with the 
low-energy measurements of Sosniak (1958), and at 82 Mev they agree with 
the results of Yule and Turkevich (1960) to within the combined experimental 
errors of the two sets of measurements. 

The smooth curves drawn through the experimental points in Figs. 2, 3, 
and 4 are presented together in Fig. 5. The two triangles in Fig. 5 represent 
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Fic. 5. The smooth curves drawn through the experimental points in Figs. 2, 3, and 4 are 
shown together. The triangles at 23 and 31 Mev represent (p,p’) cross sections (Cohen 1959). 
The dashed curve is a theoretical one for the (p,pm) reaction based on cascade-evaporation 
calculations of Jackson (1956). 


(p,p’) cross sections at 23 and 31 Mev as given by Cohen (1959). (The (,9’) 
reaction may be regarded as a (p,pxm) reaction for x = 0.) Except for dif- 
ferences in threshold, the (p,pm), (p,p2n), and (p,p3m) curves are similar. The 
(p,pm) cross-section values are somewhat higher than the others, but they 
include contributions from (p,d) pickup reactions. Pickup effects are probably 
less important so far as the (p,p2m) and (p,p3n) reactions are concerned. The 
(p,t) reaction, which is energetically favored over true (p,p2n) events by 
about 8.5 Mev, does, however, lead to a relative lowering of the apparent 
(p,p2n) threshold; the effect is visible at the beginning of the (p,p2”) curve in 
Fig. 5. Finally, it is interesting to note that a smooth curve drawn from a 
reasonable threshold through the two (p,p’) points is similar in appearance to 
the other (p,pxn) curves. 

The experimental (p,pxn) cross sections are an order of magnitude larger 
than those expected on the basis of statistical theory. These large cross sections 
are characteristic of proton inelastic scattering reactions in general, and as 
mentioned earlier, such reactions should probably be interpreted in terms of 
nucleon-nucleon collisions in the diffuse surface region of the target nucleus 
(Griffiths and Eisberg 1959), at least for energies below about 50 Mev. 
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In terms of the surface mechanism, (p,pxn) reactions would take place as 
follows. The incident proton suffers a single collision, or a small number of 
multiple collisions (Kikuchi 1960), near the surface of the target nucleus and 
the incident proton or a struck proton emerges. Prompt neutrons may accom- 
pany the prompt proton, but more likely the other nucleons involved in the 
initial interaction enter the nuclear core to form a compound nucleus. Finally, 
the compound nucleus decays by neutron evaporation. 

At higher bombarding energies the nucleus as a whole becomes more trans- 
parent and direct interactions in the nuclear core can contribute to (p,pxn) 
events. This ‘‘cascade-evaporation” process has been treated by Monte 
Carlo methods (Jackson 1956; Metropolis et a/. 1958); in these calculations the 
nucleus was assumed to have a sharp boundary, so the cross sections obtained 
indicate the importance of internal direct processes. The dashed curve in 
Fig. 5 represents such a cross-section curve for the (p,pm) reaction, calculated 
according to the results of Jackson. Direct interactions in the nuclear interior 
thus account for about a third of the observed cross section at 90 Mev, and at 
lower energies are less important. Curves for the (p,p2”) and (p,p3m) reactions 
are not shown; they are similar to the (p,pm) one but lie 10 and 20 Mev higher 
on the energy scale, respectively. 

In a process involving prompt proton emission followed by compound 
nucleus formation, the residual energy of excitation following the prompt 
events has a distribution that is a mirror image of the energy distribution of 
the scattered protons. In order for a (p,pxm) reaction to occur, the residual 
energy of excitation must be appropriate to cause the evaporation of x neutrons; 
thus the relative magnitudes of the various (p,pxn) cross sections would 
depend in a simple way on the energy spectrum of the scattered protons. On 
this basis, and using neutron evaporation probabilities published by Jackson 
(1956), the experimental cross sections shown in Fig. 5 (including the two 
(p,p’) points) are consistent over the whole range of bombarding energies 
with energy spectra that are approximately flat above the energy of the 
Coulomb barrier. 


Comparison of (p,pxn) and (p,xn) Cross Sections 

Bell and Skarsgard (1956) have published (p,x«) cross sections of several 
heavy nuclides, for a range of values of x. The individual curves have narrow 
peaks due to compound nucleus processes, followed by low tails extending to 
higher energies; these tails are due to direct interactions. A surprising feature 
of the (p,xm) results is the relatively large peak cross section for the (p,2n) 
reaction; this has the effect of causing a minimum in the (p,xm) cross-section 
sum (i.e. sum for all values of x) at about 25 Mev. There are smaller ripples in 
this curve at higher energies, but they should probably not be taken seriously. 

A comparison of the Bell and Skarsgard results for Bi?® with the results of 
the present work shows that the minimum in the (p,xm) sum curve occurs at 
proton energies for which (p,pxn) cross sections are becoming significant. It will 
be indicated in the next section (see Fig. 6) that an addition of (p,pxn) cross 
sections for Au' (including those for x = 0) to the (p,m) sum curve for 
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Bi? yields a roughly monotonic curve, and one is led to speculate that the 
marked structure in the (p,«m) sum curve is due to competition from (p,pxn) 
reactions. 

Up to about 50 Mev, incident protons that enter the nuclear core will 
probably lead to compound nucleus formation with subsequent neutron 
evaporation, or to prompt neutron emission. Competition between (p,pxn) and 
(p,xn) processes would be unimportant, so far as this mechanism is concerned. 
It is interesting, however, to consider competition between these processes for 
interactions in the nuclear surface. One might suppose that at low energies, 
say below 20 Mev, interactions in the diffuse surface also lead primarily to 
compound nucleus formation, or to prompt neutron emission, and hence to 
(p,xn) reactions. That is, the surface of the target nucleus appears relatively 
more dense at low energies, making direct interactions less probable, and direct 
interactions would be less likely to lead to proton emission due to internal 
reflections at the edge of the potential well. Thus at low energies surface 
interactions would lead to (p,x«n) events, and the dip in the (p,«n) sum curve 
at higher energies could be attributed to an enhanced cross section for direct 
interactions leading to charged particle emission. 

The (p,pxn) cross sections are about the same as (p,(x+1)m) cross sections 
in the energy region of the direct interaction tail of the (p,(x-+1)m) cross- 
section curves. For example, the (p,pm) cross section of gold is approximately 
the same as the (p,2”) cross sections of Bi? and Pb? at 35 Mev, and at 45 
Mev the (p,p2”) and (p,3”) cross sections are similar. The observed (p,2”) and 
(p,3n) cross sections can be accounted for largely in terms of internal cascades 
(Jackson 1956), and thus, as far as the surface mechanism is concerned, 
(p,pxn) reactions predominate over (p,(x-++1)m) reactions. On the basis of an 
elementary surface collision model (p,pxn) cross sections are expected to be 
about twice as large as (p,(x+1)m) cross sections (from simple statistical 
arguments), but the actual disparity is probably greater. Cohen (1959) has 
shown that at 23 and 31 Mev, (p,p’) cross sections are typically 5 to 10 times 
larger than (p,m) cross sections. A discussion of this point in relation to the 
nucleon-nucleon collision mechanism is contained in a recent paper by Eisberg 
(1960). 

The Total Reaction Cross Section 

The sum of the (p,pxm) cross sections from the present work and the pre- 
viously published (p,x) cross sections gives an approximation to the total 
reaction cross section of a heavy nucleus. The experimental (p,pxn) cross 
sections include (p,d) and (p,t) cross sections, and other reactions are probably 
relatively unimportant. (Bell and Skarsgard have shown, for example, that at 
32 Mev (p,a), (p,an), and (p,2pn) cross sections of lead are each less than 
2 millibarns.) 

The sum of the individual (,x7) cross sections of Bi?® is represented by the 
solid curve in Fig. 6. This curve was taken from the paper of Bell and Skarsgard, 
and was adjusted in accordance with a more recent beam monitoring cross 
section (Whitehead and Foster 1958). The small circles show the smoothed 
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Fic. 6. The solid curve represents the sum of (p,xn) cross sections of Bi? (Bell and 
Skarsgard 1956). The circles represent this sum plus the (p,pxm) cross sections from the present 
work. The triangles at 23 and 31 Mev include the (,p’) cross section (Cohen 1959), and the 
dashed curve drawn through these points approximates the total reaction cross section of a 
heavy nucleus. At higher energies the circles give a rough lower limit to the total reaction cross 
section. The crosses represent directly measured total reaction cross sections of lead (from 
Gooding 1959; Meyer et al. 1960; and Goloskie 1961). 


(p,pxn) results from the present work added to the (p,«m) curve. The triangles 


include the (p,p’) cross sections reported by Cohen (1959), and the dashed 
curve results from the addition of a reasonable (p,p’) curve based on these 
two values. The dashed curve approximates the total reaction cross section, 
with a standard error of about 15%. The small circles represent a rough lower 
limit to the total reaction cross section at higher energies, since the (p,p’) and 
(p,p4n) cross sections are not included (the (p,p5n) and (p,p6m) cross sections 
of bismuth are contained in the sum curve presented by Bell and Skarsgard). 

Three points representing directly measured total reaction cross sections of 
lead are shown in Fig. 6 as crossed error bars. The point at 34 Mev is due to 
Gooding (1959), the one at 62 Mev to Meyer et a/. (1960), and the point at 
77 Mev is from Goloskie (1961). The values at 34 and 77 Mev are in good 
agreement with the present work, but that at 62 Mev appears to be somewhat 
low. 

The total reaction cross section of a heavy nucleus is, to within experimental 
accuracy, a monotonic function of incident proton energy. Bell and Skarsgard 
were led toward a different conclusion because they did not know that (p,pxn) 
cross sections were so large; the present work does, however, support their 
conclusion that near the peak of a (p,xm) cross section, the (p,p(x—1)n) 
cross section is small by comparison. 

In this section and the preceding one we have compared and combined 
(p,pxn) cross sections of gold with (p,m) cross sections of bismuth and lead. 
This procedure is of course not strictly correct, but cross sections are not 
expected to vary strongly among heavy elements. In the addition of gold and 
bismuth cross sections, the gold results were shifted slightly toward lower 
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energies to account for the difference in neutron binding energies. Bismuth 
(p,xn) cross sections were chosen since Bi?®, like Au™, is an odd A, odd Z 
nucleus. Since the total (p,px”) contribution to the cross-section sum shown 
in Fig. 6 is only about 20%, the addition procedure is justified to well within 
the errors in the experimental measurements. 


One of us (T. M. K.) wishes to thank the National Research Council for 
scholarships received during the course of this work. 
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THE RATIO OF THE RESONANCE INTEGRAL TO THE 
THERMAL NEUTRON CROSS SECTION FOR Sm?! 


W. H. WALKER AND R. E. GREEN 


ABSTRACT 


Cadmium ratio measurements have been made in similar lattice positions in 
ZEEP with thin foils of Sm and gold. From a comparison of these cadmium 


ratios we find 
(I'/goo) smise = 14.65+0.41. 


If goo is assumed to be 212+12 barns, then J’ = 3100+200 barns. On the 
assumption that only one resonance in Sm’®, at 8 ev, contributes appreciably to 
both the resonance integral and the thermal cross section, it follows that g = 1 
and that T, the width at half-maximum of the total cross section resonance, is 
193-£5 milli-ev. Because of the appreciable disagreement between this value and 
one reported earlier, a new time-of-flight measurement of the resonance 
parameters has been made, and it is reported in the next paper (Chrien 1961). 


INTRODUCTION 


Although the presence of a large resonance at 8 ev in Sm!” was first noted in 
1954 (Sailor, Landon, and Foote 1954) only recently have the parameters 
become generally available (Hughes et al. 1960). Because the neutron scattering 
width is unusually large for a resonance energy of the order of a few electron 
volts, it has not been possible to determine either width to the degree of 


accuracy usually associated with time-of-flight measurements in this energy 
region. 

Even before the resonance was located, activation measurements were 
made of the epicadmium resonance integral (Sunyar and Goldhaber 1949) 
and of the ratio of the resonance integral to the thermal cross section (Harris 
et al. 1950). More recently the resonance integral and thermal cross section 
have been determined using the Harwell pile oscillator (Tattersall et al. 
1960), and a crystal spectrometer has been used to measure the 2200 m/sec 
cross section (Pattenden 1958). The measured cross sections and values 
calculated from resonance parameters cover such a wide range that a further 
measurement of the cadmium ratio using thin foils seemed justified. A tabula- 
tion of the earlier results is included in Table III. 


THEORY 


The theoretical basis for relating the cadmium ratio and I’/gao has been 
discussed in detail elsewhere (Walker e¢ al. 1960). Only an outline sufficient 
for the present experiment will be given here. We assume that the neutron 
flux can be resolved into a Maxwellian component and an epithermal com- 
ponent having a 1/E distribution, that the radiative capture cross sections of 
the target nuclei follow closely the 1/v law in the subcadmium region, and that 

‘Manuscript received May 18, 1961. 
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the cadmium cover is sufficiently thick to eliminate all thermal neutrons. Then, 
for two different nuclei, A and B, irradiated in the same spectrum, 






a) (FRea=1)a _ (Gel’/go0)n+/4/2K 
(FRea—1)p  (Grl /goo)at+ V0 /2K ’ 


where Rog is the measured cadmium ratio, in the present case the ratio of the 
activity of a detector wrapped in aluminum to that of a detector wrapped in 
cadmium, 

ay is the cross section for neutrons with a velocity of 2200 m/sec, 

goo is the effective absorption cross section for a Maxwellian neutron 
distribution with characteristic temperature 7, 

I’ is the “reduced” resonance integral, i.e. the usual absorption integral 
fo,(E) (dE/E) over the entire epithermal spectrum minus the contribution 
of the 1/v extrapolation of the ‘‘Maxwellian”’ cross section goo, 

F is a factor accounting for the reduction of resonance neutrons due to the 
high energy wing of the cadmium resonance at 0.2 ev, 

G, is the resonance self-shielding factor for detectors that are not infinitely 
thin, and 

1/K is proportional to the 1/v activation above the cadmium cutoff. 

Since o¢ and J’ are known for gold (a summary of experimental data is 
given by Walker et a/. 1960), a comparison of the cadmium ratios of gold and 
Sm!” will yield a value of J’/goo for Sm’ provided the other parameters can 
be estimated. 

























EXPERIMENTAL PROCEDURE 


The sample material used in these measurements contained >98% 
samarium. The main contaminants were Cb (<0.1%), Pb (<0.1%), Zr 
(<0.1%), Ce (<0.5%), Pr (<0.2%), and Eu (~0.1%). The samarium was 
highly enriched in Sm'; the relative abundances of the various samarium 
isotopes are shown in Table 1. 










TABLE I 


Relative abundances of samarium isotopes 
















Isotope Abundance (wt.%) 
Sm!# 0.03 
Sm? 0.20 
Sms 0.20 
Sm49 0.27 
Sm}50 0.30 
Sm! 97.20 










Sm 1.80 






A solution of samarium nitrate in lacquer was used in the foil preparation. 
The samarium was deposited as samarium oxide (Sm,Q;) on a 0.0025 cm thick 
sheet of spectroscopically pure aluminum by painting on a thin coat of lacquer 
containing about 8 wg/cm? of samarium and then burning off the lacquer by 
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‘firing’ in an oven at ~500° C. This process was repeated until a layer 
containing ~140 ywg/cm? of Sm’ was built up. From this foil six 1.31 cm 
diameter disks were cut. 

The average deposit thickness was estimated from the total amount of 
material used, assuming that the deposit was evenly distributed over the 
painted area. In a target of Pu™° prepared in a similar way (Walker et al. 1960) 
scanning with an a-probe showed that no local deviation exceeded 15%. The 
relative amounts of Sm’ on each of the six foils were determined, by activation, 
to have a standard deviation from the mean of 9%, and a maximum spread of 
26% on the mean. On this basis we conclude that the mean deposit thickness, 
140 ug/cm?, is correct to +13%. 

The gold target material was prepared by vacuum evaporation of pure 
gold onto spectroscopically pure aluminum foils. Average deposit thicknesses 
were determined by weighing the backing foil before and after deposition. Six 
gold targets were cut; four of these were 1 cm by 1 cm square foils, ~450 ug /cm? 
gold thickness, and the other two were disks, ~210 yug/cm? gold thickness 
and the same diameter as the samarium foils. The maximum spread in the 
relative amounts of gold in these foils was found to be 10% by activation. The 
mean deposit thicknesses are therefore assumed to be correct to +5%. 

For the cadmium ratio measurements the foils were placed in aluminum and 
cadmium capsules having 0.076 cm thick walls. These capsules were mounted 
on thin aluminum channels which were, in turn, inserted into aluminum 
thimbles attached to the ZEEP reactor lid (see Fig. 1). The thimbles used in 
these irradiations were D and G, located at cell boundary positions in the 
lattice on opposite sides of the pile and equidistant from the center. They 
were not filled with moderator. 

In addition to the central assembly, which was changed from time to time, 
the lattice in which these measurements were made consisted of the following: 
a core of eighteen 19-element U—UQO, clusters (each cluster contained fifteen 
1.32 cm diameter UO, elements and four 1.31 cm diameter U elements), an 
intermediate region of eighteen 19-element UO: clusters, and an outer ring 
of 42 ZEEP rods (3.26 cm diameter U rods) which acted as a driver region for 
the lattice. The cluster-type rods formed a hexagonal array with a center-to- 
center spacing of 24.6 cm. (For further details of these fuel elements see Hone 
et al. 1958.) The central site contained either a loop assembly (Bigham e¢ al. 
1961) or a normal 19-element UO, cluster, or was empty at various times 
throughout these experiments. Our measurements were always made in 
thimbles D and G in the U—UO: cluster region (see Fig. 1). 

The usual procedure was to put the cadmium-covered foils in one thimble 
and the aluminum-covered ones in the other with the foils of a given type at 
the same elevation. At other times both foils of one type were mounted in the 
same thimble at different elevations spaced so that the cadmium capsule did 
not depress the thermal flux at the position of the aluminum-covered detector. 
Copper foils were also placed at various elevations in both thimbles to check 
any radial asymmetry (always <1%) and to map the axial flux distribution 
for the case when the foils were at different elevations. 











1188 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 





Since the fuel rods in the U-UO, cluster region of the pile had 51 cm long 
mock-ups of pressure and calandria tubes (see Fig. 1) the epithermal fraction 
of the neutron density was not constant all along the axis of the pile but was 
higher in the region of the double-walled tubes. For most of the measurements 
reported here the foils were placed well away from these tubes but for the 
early runs some of the foils were in or near this region. The latter measurements 
have been corrected for the effect of the double-walled tubes using epithermal 
neutron densities determined from the cadmium ratios of 0.025 cm diameter 
gold wires placed along the axis of the pile. 

The irradiated detectors were counted with a system which included two 
Nal(T1) scintillation counters and an automatic sample changer and readout 
unit. The counters were biased to accept all pulses corresponding to y-ray 
energies >30 kev. The data were recorded on typewritten sheets and on 
punched tape and were reduced by feeding the latter into the Burroughs-205 
computer at Chalk River, which made corrections for decay, counter deadtime, 
and background. 

The resulting cadmium ratios are listed in Table II. In order to obtain 


TABLE II 
Measured values of Rca. Derived values of (G,I'/goo) gm! 














Irradiation (Rea) sm'5? (Rea) au (G,I'/gao) Sm!52 
1 3.099 2.942 13.82 
2 3.051 3.030 13.902 
3 3.124 2.958 13.88 
4 3.162 3.061 14.322 
5 3.118 3.097 14.38 
6 3.135 3.024 13.892 
7 3.096 3.041 14.12 
8 3.120 3.018 13.90 
9 3.130 3.023 14.07 

10 3.124 3.056 14.05 
11 3.150 3.009 13.84 
12 3.141 3.045 13.85 


Mean value of (G,I’/gao) gm82 = 14.00+0.06 


@These values were corrected for differences in the epithermal fraction of 
the neutron density at various gold and samarium detector sites. 


G,I'/goy for Sm'” from the measured ratios, it is necessary to know the value 
of K appropriate to the irradiation conditions, values of F for both gold and 
Sm! for the thickness of cadmium used, and the value of G,I’/gao for gold. 
Values of K have been calculated for a 1/v detector in a 1/E spectrum using 
the known variation of the cadmium cross section with energy (Westcott et al. 
1958). For an infinite slab (or thin foil) in an isotropic flux, and 0.076 cm 
(0.030 in.) of cadmium, K = 2.20. Values of g have been calculated for Max- 
wellian distributions having temperatures from 20°C up (Westcott 1960). 
The effective neutron temperature in the moderator of this ZEEP lattice is 
known to be about 10° above moderator temperature, i.e. ~30° C (Bigham 
and Tunnicliffe 1961), so that for these irradiations gay = 1.006. 

The capture by cadmium of neutrons at resonance energies can be computed 
from the parameters of the 0.2-ev resonance, assuming that contributions 
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from other resonances are negligible. The cadmium capture cross section at 
the gold resonance energy, 4.9 ev, is 0.2 barn, and at 8.0 ev, the Sm! resonance 
energy, it is about 0.06 barn. The corresponding F-values for 0.076-cm 
cadmium are 0.994 and 0.998. 

For resonances in which capture predominates, as in the case of gold, 
resonance self-shielding can be calculated using the method described by Roe 
(1954). For the gold foils used here we estimate a value for (G,)4y of 0.951 
for the 450 ug/cm? foils and 0.971 for the 210 ug/cm? foils. Using these factors 
and taking (¢0)a4u = 98.6 barns and (J’),4y = 1525 barns (Walker et al. 1960), 
the measured cadmium ratios can be used to derive the values of G,/’/goo for 
Sm!» listed in the last column of Table II. Those results which required correc- 
tion for differences in the value of the epithermal neutron fraction have been 
noted; in every case the correction was less than 3% in the ratio (G,1'/goo) gms2. 
The mean value of (G,/’/goo)smis: for 12 irradiations is 14.00+0.06; the error 
is based on the reproducibility of the individual measurements. 

In order to derive I’/gao for Sm’ we require (G,) sms:. A straightforward 
calculation using Roe’s method is not possible since in Sm'® the probability 
of scattering is of the same order of .magnitude as that of radiative capture, 
and the effect of resonance scattering is not well understood. In the estimates 
which follow we use the recently measured values of the resonance parameters 
(Chrien 1961). An upper limit of 0.955 for (G,) sms: is obtained if resonance 
scattering is assumed to have no effect; on the other hand, if scattering is 
considered to be just as effective as absorption then the value is 0.94. However, 
recent work with manganese detectors indicates that scattering is about 
one-quarter as effective as absorption, whence (G,)smis: = 0.95; this is the 
value we will adopt here. 

A further correction is required to allow for the fact that the epithermal 
spectrum is not exactly 1/E in this energy region. Calculations of the slowing- 
down spectrum have been made for a fuel arrangement very similar to the one 
used in these experiments, namely, 19-element uranium oxide clusters in a 
24.6-cm hexagonal array (Bigham and Pearce 1961). These calculations 
indicate an average deviation of <0.2% per ev in the region 3-8 ev, in a 
direction which makes the flux smaller as E decreases. Since the resonance 
activation is proportional to the resonance flux times J’, a correction factor of 
0.994 must be applied to our result. 

The average value of the ratio of the reduced resonance integral to the cross 
section in a Maxwellian spectrum, (1’/goo) sm1s2, will also have the following 
errors in addition to the one assigned in Table II to the average ratio 
(G,I’/goo) gmis?: 

(i) The error in the ratio (J’/a0)au, 2.6%. This is due almost entirely to 
uncertainties in the resonance parameters, which give a calculated reduced 
resonance integral J’ of 1525+40 barns (Walker et al. 1960). 

(ii) The error in the self-shielding factor (G,) gms, 0.7%. For foils of this 
thickness G, is relatively insensitive to errors in the resonance parameters or 
estimates of foil thickness. The 13% uncertainty in foil thickness gives rise 
to an error of 0.13 X0.04, or 0.5% in (G,) gmisz. The error due to uncertainty in 
the effect of resonance scattering is taken to be about 0.5%. 
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(iii) The error in the self-shielding factor (G,) au, 0.8%, due to the uncertain- 


ties in the deposit thicknesses. 

(iv) The uncertainty in the correction for the fact that the flux is not 
exactly 1/E, 0.5%. 

(v) Errors due to uncertainties in the quantities gay, Fau, Fsms:, and K 


are all negligible. 
The derived value for the ratio of the reduced resonance integral to the 


cross section in a Maxwellian spectrum is therefore 
(I'/goo) sms: = 14.65+0.41. 


This result is compared in Table III with values derived from earlier cross- 
section measurements. 















TABLE III 
Comparison of results of cross-section measurements for Sm® 





Resonance integral Thermal cross section 
, 


















&o0 
Source (barns) (barns) I'/ga0 

Sunyar and Goldhaber (1949) 2400 — ee 
Harris et al. (1950) —_— — 14 
Pattenden (1958) —_— 200+6 — 
From resonance parameters 

(Hughes et al. 1960) 2100 100 21° 
Tattersall et al. (1960) 2850+300 22447 12.7+1.4 
This work —_ = 14.65+0.41 













“The original data have been corrected for resonance self-shielding in the thick samarium detector, estimated 


to be a factor of 0.8. 
’Contributions from the only known resonance, at 8 ev. 








DISCUSSION 


It is unlikely that any other level, virtual or bound, is contributing signi- 
ficantly to the cross section in the thermal region, since Sm", an even-even 
nucleus, is expected to have an average level spacing of several hundred 
electron volts. This assumption is supported by the crystal spectrometer 
measurements of Pattenden (1958), who found no significant deviation from 
a 1/v variation in this region. We therefore assume that ggmis: = 1.0 to better 
than 1%, and (I’/o0) sms? = (I’/gae) sms? = 14.65+0.41. 

Two recent measurements of the thermal cross section give oo = 200+6 
barns (Pattenden 1958) and goo = 224+7 barns (Tattersall et al. 1960). In 
view of the discussion above, it is unlikely that any significant part of the 
discrepancy can be assigned to g. Assuming a mean value op = 212+12 then, 
we obtain from our measured ratio a value for the reduced resonance integral of 


I’ = 3100+200 barns, 


















where the error is predominantly due to the uncertainty in the 2200 m/sec 
cross section. 

Because the level spacing in Sm!® is expected to be several hundred electron 
volts, it is probable that the 8-ev resonance dominates the resonance integral 
as well as the thermal cross section. In this case it can be shown that 
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,_ oly 
[= OR. +small terms, 


Gt, 


4 VE, /Eo+small terms, 


go = 


I'/oo = 20 VE OE mall terms, 


where I is the width at half-maximum of the total cross-section peak, 

I’, is the partial width for radiative capture, 

E, = 0.02530 ev, the energy of 2200 m/sec neutrons, 

E, is the energy at the resonance peak, and 

o, is the total cross section at the peak of the resonance, excluding 
potential scattering. 
T, l',, and E£, are in electron volts. For the 8-ev resonance in Sm! the small 
terms in J’ and ao are about 0.2% and 0.1%, respectively, of the first term, so 
that for the ratio J’/ao the error in neglecting the small terms is about 0.1%. 
Thus, 


T = 24 \/EcE; + = 19345 milli-ev. 


I 


This value for [ disagrees appreciably with the value inferred from the 
compilation by Hughes et al. (1960). They quote T,, = 79+8 and l', = 57+30, 
so that Tl = 136+31 milli-ev. On the basis of this disagreement a time-of-flight 
measurement of the 8-ev resonance was proposed, using the BNL - Chalk 
River chopper facility. The results of this investigation are reported in the 
succeeding paper. 
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RESONANCE PARAMETERS FOR THE 8-ev LEVEL OF Sm? 


R. E. CHRIEN 


ABSTRACT 


The joint BNL-AECL Fast Chopper Facility operating at the NRU reactor 
at Chalk River has been used to determine the resonance parameters of the 
8-ev level of Sm. Using a resonance shape analysis method, the following 
parameters were obtained: Ey = 8.02+0.02 ev, Tl = 205+15 milli-ev, and T, 
= 79+3 milli-ev. With these values, a reduced resonance capture integral of 
3090+220 barns is calculated, in agreement with recent cadmium ratio measure- 
ments. 


INTRODUCTION 

There is considerable interest in the resonance capture properties of the 
nuclide Sm? because of its importance as a fission product, particularly for 
Pu*® fission. Measurements of the resonance capture integral for this nuclide 
are in disagreement with a value calculated from previously reported resonance 
parameters for the 8-ev level (Hughes et al. 1960), a large level which dominates 
the resonance integral of Sm!®. Recent measurements at Harwell (Tattersall 
et al. 1960) and at Chalk River (Walker and Green 1961) lead to reduced 
resonance capture integrals of 2850+300 barns and 3100+200 barns, re- 
spectively, compared with 2100 barns as computed from resonance pa- 
rameters. This discrepancy suggested a remeasurement of the parameters of 
the 8-ev level with the BNL—AECL Fast Chopper Facility at Chalk River. 


EXPERIMENT 

The BNL-AECL Fast Chopper Facility consists of the BNL-type fast 
chopper operating in conjunction with an 88.2-meter flight path at the NRU 
reactor. The detector is made up of an array of 384 counter elements in a 
common atmosphere of BF; gas. The time analysis of the neutron counts is 
done with the Brookhaven 1024-channel time analyzer. A detailed description 
of the Chalk River Facility will soon appear elsewhere (Zimmerman et al. 
1961). Except for minor instrumental changes the chopper is identical with 
the one in operation at the Brookhaven research reactor (Seidl et al. 1954). 

A sample of samarium enriched to 97.2% Sm! was made by preparing a 
solution of samarium nitrate in a lacquer base and applying it to an aluminum 
backing according to the prescription given by Walker and Green (1961). In 
the present experiment 15 successive layers of Sm! were built up using various 
brush-stroke directions. Strips of approximately 5X2cm were cut from the 
center of the foil and four such strips were stacked to give a sample thickness 
of 4.60 X10'8 atoms/cm? of Sm?, 

The chopper was run at 2000 r.p.m. giving a repetition rate of 267 bursts/ 
sec with the eight-armed rotor. 

1Manuscript received May 18, 1961. 
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The use of an 88-meter flight path with a rotor designed for 20-meter 
operation introduces the problem of ‘‘overlapping”’ of neutrons from successive 
bursts. The low-energy cutoff of the rotor is such that fast neutrons from the 
previous burst overtake slower neutrons before they reach the detection 
station. In the present experiment the use of a boron filter reduced the overlap 
background component to an amount much smaller than ‘‘room’’ back- 
ground. The amount of overlap was determined by running a thick (j in.) 
sample of U?* and observing the counting rate in the region of resonances, 
where, for a sample of this thickness, the transmission is very nearly zero. 

The untreated data of this experiment are shown in Fig. 1. About 4000 
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Fic. 1. The counting rate curve in the neighborhood of the 8.02-ev level of Sm, The 
total background indicated is the sum of the room background and the overlap components. 
This curve is the summation of several runs totalling 7639 minutes’ duration. 


counts per 2-ysec channel were recorded in the neighborhood of the resonance 
in several runs of 7639 minutes’ total duration. The room background, that 
is, the beam-independent background, and overlap backgrounds are indicated 
in the figure. The counting rate in the ‘‘open beam’”’ or sample out position 
as shown is not measured, but is computed by fitting a least-squares straight 
line to the counts in the neighborhood of the resonance. This line was then 
used for conversion of the data to transmission form. 


ANALYSIS 
Because of the high resoiution available, a ‘‘shape’”’ analysis of the 8-ev 
level was performed. The resolution function of the chopper was calculated 
by compounding the following contributions: a triangular chopper burst of 
5.0 usec full width at half-maximum, a time uncertainty of 4.8 usec resulting 
from the flight time of an 8-ev neutron in the detector, and an analyzer 
channel width of 2.0 usec. The compounding of these contributions results 
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in a function which is very nearly Gaussian in form with a standard deviation 
of 2.51 usec. This corresponds to 52.5 milli-ev and is to be compared with a 
Doppler width of 74 milli-ev and a level width of about 200 milli-ev. Thus the 
resolution corrections on the observed shape are quite small. 

The Doppler-broadened single level formula tables of Sailor (1953) were 
used for the shape analysis. The observed transmission curve was compared 
with a series of resolution-broadened, Doppler-broadened curves for various 
values of I’, with the total area of the resonance held fixed to the experi- 
mentally observed value. That value of I was selected which corresponded 
to a minimum in the sum of the squares of the deviations between the observed 
and expected data points. The best fit was obtained for the following set of 
parameters: 

T = 205+15 milli-ev, 
T,, = 79+3 milli-ev, 
Eo = 8.02+0.02 ev. 


Figure 2 shows the fit of the curve based on these parameters to the experi- 
mental points. 
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Fic. 2. The result of the shape analysis for the 8.02-ev level. The solid line is the resolution- 
broadened, Doppler-broadened Breit-Wigner single level curve for the set of parameters 
shown. 
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CONCLUSION 


With the reasonable assumption that only the 8-ev level is important in 
determining the capture properties of Sm!, the following quantities may be 
computed from the presently reported resonance parameters for comparison 
with experiments done with pile neutrons (for precise definitions of the 
following quantities refer to Walker and Green (1961)): 

(1) The reduced resonance capture integral, 

I’ = 3090+220 barns. 
(2) The thermal capture cross section, 
oo = 224+18 barns. 
(3) The ratio of the resonance integral to the thermal neutron cross section, 
I’/oo = 13.8+1.0. 

The values of parameters (1), (2), and (3) are in good agreement with 
recent pile neutron measurements. A detailed listing of these previously 
obtained parameters is given in the accompanying paper of Walker and 
Green. 

The value of I, inferred from the present experiments is 126+15 milli-ev, 
in disagreement with the value of 57+30 quoted in the recent supplement 
to BNL 325 (Hughes et al. 1960). The value of 126415 milli-ev is signi- 
ficantly greater than that measured for Sm47, Sm™9, and Sm"™!, which all have 
radiation widths of about 60 milli-ev. 
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ALFVEN WAVES IN PARTIALLY IONIZED GASES! 


ToMiyA WATANABE? 


ABSTRACT 


The conditions for a wave, propagated in a partially ionized gas along an 
external magnetic field, to be of Alfvén type have been obtained. 


1. INTRODUCTION 

It is well known that hydromagnetic waves can be propagated in a good 
conducting medium permeated by an external magnetic field. For instance, a 
sinusoidal transverse wave can be propagated along the external uniform 
magnetic field, which is characterized by the following two properties (Alfvén 
1950): 

(i) The velocity of propagation of the wave is given by 
(1.1) V = Ay//4rp, 


where Hy is the intensity of the external magnetic field and pp is the density 
of the fluid. 

(ii) The density of the kinetic energy of the fluid is equal to that of the 
magnetic field energy. 
These particular waves, called Alfvén waves, are the simplest kind of hydro- 
magnetic waves. 

Mathematically speaking, the preceding are inevitable results from the 
equation of motion of the fluid, 


du 


(1.2) Poa 


j A Ho, 


1 
c 


and the law of electric conduction in an infinitely good electric conductor, 


(1.3) E+" uA Ho aa 0, 


as well as Maxwell’s equations of the electromagnetic field. Equation (1.3) 
is the extreme case of the following ohmic law as the conductivity o tends 
to infinity: 


(1.4) £ogt+iu nM 
o c 


It assures the condition of “freezing in’’ of magnetic lines of force (Cowling 
1957). Any volume element of fluid which is extremely conducting cannot 
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move across magnetic lines oi force, because no infinitely intense electro- 
motive force may be induced. Due to this “freezing in’’ condition, a tube of 
magnetic lines of force can be simulated by a tensile string, since the Maxwell 
stress can be resolved into the isotropic pressure H?/8x, and the tension 
H?/4x along a magnetic line of force. This analogy leads to the velocity of 
propagation given by (1.1) (Cowling 1953). 

Certainly, eq. (1.3) is mathematically idealized. Actually, there is no 
material of infinitely good conductivity. One of the conducting fluids of 
particular interest is an ionized gas. For ionized gases, the equation of motion 
(1.2) and the law of electric conduction, eq. (1.3), do not always hold true. 
The condition under which Alfvén waves can exist in a completely ionized 
and very rarefied gas have been sought by Astrém (1951), Schliiter (1952), 
and Ferraro (1955). Waves in slightly ionized gases have been studied by 
Hines (1953), and in relation to the ionosphere, by Dungey (1954a), Watanabe 
(1957), and Fejer (1960). 

One of the aims of this paper is to derive the conditions for the occurrence 
of Alfvén waves. The approach is not to manipulate the dispersion relation 
as done by most of the authors, but to discover the conditions under which 
the equation of motion (1.2) and the ohmic law of perfect conductivity (1.3) 
hold true. Another aim is to find out the exact expression for the rate of 
dissipation of energy of Alfvén waves. 

In Sections 2 and 3, the author will derive fundamental equations to describe 
the electrodynamical behavior of the partially ionized gas and several other 
useful relations. The conditions for the occurrence of Alfvén waves will be 
sought in Section 4, and the expression of dissipation rate of the waves in 
Section 5. Application to the ionosphere will be discussed in Section 6. 

2. FUNDAMENTAL EQUATIONS OF WAVES IN PARTIALLY IONIZED GASES 


Let us consider an ionized gas which consists of electrons, ions, and neutral 
particles. Ions are assumed to have been singly ionized. The number density 
of the electrons is assumed to be equal to that of the ions, and it is represented 
by N,. Then, the equations of motion of the electrons, ions, and neutral 
particles are given by the following, neglecting the pressure terms (Schliiter 
1951): 


7 1 
(2.1) Pe Ue + Npotie(Ue— U1) +NpNaten (We Un) = — EN, (E+— ue A H)), 


ou 2 1 
(2.2) Pi “ap + Noette(Ui— We) +NpNnetin (Us — Un) eN, (E+ u,A Ho), 


dUn 
(2.3) Pa 5p + NpNaetin(Un—W1) +N pNacen(Un—Ue) = 0. 


The quantities with subscripts e, 7, and m refer to the electrons, ions, and 
neutral particles. pe, p;, and p, represent the mean densities of electrons, ions, 
and neutral particles respectively. N, is the number density of neutral par- 
ticles. U., U;, and u, represent the mean velocities of the electrons, ions, and 
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neutral particles. Hy is the intensity of the external magnetic field. E is the 
induced electric field. The coefficients aje, dn, and ay are defined by 


— Vie_MiMe_ _ Viele 
(2.4) ai N, m,+m,. N, ’ 
Ven ™MeMn_ __ VenMe 
(2.5) Cn = Rica <a e 
— Yin MM _ Vin 
(2.6) Ca Namitie Ny’ 
(2.7) Vin = 3¥ in- 


vie is the collision frequency of an ion with electrons; ve, and vy», have similar 
meanings. me, m,, and m, are the masses of an electron, ion, and a neutral 
particle. The ion mass is assumed to be equal to that of a neutral particle, 
i.e. m;, = m,. The mass of an electron is assumed much smaller than the ion 
mass. The second and third terms of each equation from (2.1) to (2.3) represent 
the forces due to collisions between the different kinds of particles. Since we 
are considering only the case of transverse waves, we may obviously omit 
the pressure terms. The amplitudes of waves are assumed to be so small that 
all the quantities may be linearized. 

From eqs. (2.1), (2.2), and (2.3), we can derive the following equations 
(cf. Schliiter 1951): 


ou, ~ is Vie 


(2.8) Pn "Ot = a ae? a" o Pp¥2(U,—Up), 
du, __ Hon; a Ao 
(2.9) Pp _— ¢ We j- Pp¥2(Up un)+ C jAl, 


where 1 is the unit vector parallel to the external magnetic field. p, and u, 
are the density and the mean velocity of the ionized part of the gas, which 
we shall call here the plasma gas: 


(2.10) Pp = PitPe, 
(2.11) u, = Pili t pelle . 
Pp 


The electric current density is given by 
(2.12) j = eN,(ui—u,). 


u, and u, are related to j and u, as follows: 


ce 1 Pes 
(2.13) u,;= Ut oN, bp) 

_ 1 Bs 
(2.14) u. = Up— oN st, 
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The third term on the right-hand side (abbreviated herein as r.h.s.) of eq. 
(2.9) represents the ampére force, which of course acts on the plasma gas 
only. The second term of each equation (2.8) and (2.9) may be considered 
as an internal friction force between the plasma gas and the neutral particle 
gas. Writing the equations of motion in the forms of (2.8) and (2.9) implies 
that the behavior of a partially ionized gas can be described by a model of 
two kinds of fluids, the plasma gas and the neutral particle gas, with the 
mutual interactions between them due to both collisions and electromagnetic 
field (Schliiter 1951). 

By defining the density and the mean velocity of the whole gas, p and u, 


as follows, 


(2.15) P = PptPn; 
(2.16) u = Peller etn 
p 


the equation of the whole gas becomes 


ou ] 
aod.) — =-—jA Ho. 
ain re ee 

One can also derive the following equation, which may be called the general- 
ized ohmic law (cf. Schliiter 1951): 


aj : . 0 Hy vy 
(2.18) 2d (Fe rbted A 1) = E+‘ u, A 1+ S = (u,—wu,). 


Wp 


The coefficients involved are defined bv 


9 2 4 rNy,e™ . 
(2.19) y= @p = plasma frequency, 
Me : 
eH ‘ 
(2.20) oe =—, we = electron gyvro-frequency, 
MC , : 
9) « eH * 7 
(2.21) o,=—, w; = 10n gvro-frequency, 
mC : 3 
(2.22) Mina Yen —Viny 
oe , , Me 
(2.23) Vo = Vin +— Ven, 
my; 
9° Me 
2.24) V3 = VietVen to Vin- 
i 


The terms of eq. (2.18) are defined as follows. The first term in the left- 
hand side (abbreviated herein as |.h.s.) may be called the term of plasma 
oscillation, i.e., the plasma oscillations come from this term. The second and 
third terms in I.h.s. are the terms of ohmic and Hall electric fields. The second 
term in r.h.s. is called the dynamo field. The last term in r.h.s. is called here 
the term of friction electromotive force (abbreviated as the friction e.m.f.). 











WATANABE: ALFVEN WAVES 1201 


3. SOME USEFUL RELATIONS 


Let us consider sinusoidal transverse plane waves propagated along the 
external and uniform magnetic field. Wave quantities of the first degree may 
be assumed to depend on time and space by only one factor 


(3 1) etwttk.n) ; 


Therefore, the operator 0/dt on the wave quantities of the first degree may 
be replaced by iw. Hence, the following result: 


T 1 v 
39 3, 470Pn ot E : \ th pts) 1], 
(3.2) —P aoe a 1-i® rn jAl 
3 a [swt ; intl 
(3.3) Gm ee j+i2 jal}, 


by making use of eqs. (2.8) and (2.9). Then, by use of eq. (2.18), the following 
are derived: 


(3.4) u, = —-2— ows = u-i®2 uA »|, 


pvo+ paw Pp V2 Pp We V2 


(3.5) u, = —#—|u a> ,, (a A |. 


Using eqs. (3.2) and (3.3) or eqs. (3.4) and (3.5), we have 


ee Be E + Pn ¢s | 
(3.6) 8 ee — (GjA1)], 
or 
(3.7) u,—u, = i* —©2 — #| ay 2 ~ (uA » |. 
Pp pv2+ipaw voLp We 


With eqs. (3.6), (2.18), and Maxwell’s fundamental equations, we can show 
that j must satisfy the relation 


(3.8) aj = Bj A 1, 
where 
(3.9) a = —(iwtv;) +i" = teas $8) ae Pe ; Po ¥2\ 
a V~  pvottppw (\we a: a 
and 

= 2 pnw 1 
(3.10) B= meee 


In order that j may not vanish identically, the following relations should be 
satisfied : 


(3.11) 8 = +ia 
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and accordingly we have 
(3.12) j= Fij A 1. 


This shows us that there are two kinds of circularly polarized waves. In one 
of them, eq. (3.12) with + sign, all the field vectors rotate in the same sense 
that the electrons gyrate. This kind of wave is called the extraordinary wave; 
the other kind is the ordinary wave. 


4. CONDITIONS FOR ALFVEN WAVES 


Equation (3.11) is the dispersion relation. One could examine the charac- 
teristics of the waves with this relation as done in Hines’ work (1953), and 
further seek the condition under which waves may become Alfvén waves. In 
this paper, however, the author will seek the conditions in a way so as to 
reduce the equation of motion and the equation of the generalized ohmic 
law to eqs. (1.2) and (1.3) respectively. 


A. SLIGHTLY IONIZED GASES; kp = p,p/p K1, pa = p 


A(a) Strong Magnetic Field, w% > 
In this case, eqs. (3.2) and (3.3) or eqs. (8.4) and (3.5) tell us that 


ju,| > |u,l, if w > Kkpr2, 


and 


lu,| < |u,|, if w K Kkpyr. 


A(a)i. w > Kkprr 

In the frequency range w > kpv2, the importance of the friction e.m.f. relative 
to the dynamo field is of the order of v1/a., which is always smaller than 1 in 
this case. In the equation of motion of the plasma gas (2.9), the first term 
of |.h.s. is always negligible compared with the third term. The friction force 
is negligible compared with the ampére force only if 


(4.1) w > v2. 
Under these conditions, the equation of motion of the plasma gas is reduced 
to 


du, _ Ho 


2 —- = 
(4.2) Poa) 7 IAI, 


and with this relation we easily find that the relative importance of the 
plasma oscillation term, the ohmic electric field, and the Hall electric field 
term to the term of the dynamo electric field becomes w?/w we, wv3/w we, and 
w/w, respectively. Therefore, if the following conditions as well as condition 
(4.1) are satisfied, 

(4.3) w K vV/wiwe , 

(4.4) w K wywe/V3, 


(4.5) w Koi, 
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the dynamo field predominates over the other fields, and thus, we have the 
relation 





(4.6) E+" u, A Hp = 0. 





Since w; < we, condition (4.3) is automatically satisfied if condition (4.5) is 
satisfied. Conclusively, we can say that both eqs. (4.2) and (4.6) hold true 
for the frequency range 


(4.7) ve Kw K min [4 4, wwe/vs], 








and, therefore, Alfvén waves will be propagated along the external magnetic 
field with the velocity (cf. Fejer 1960) 






Ho 


4p 


(4.8) Vi 





A(a)ii. w K kpve 
In this case, we have 










up = Uz = U. 






Replacing u, in the dynamo field term of the generalized ohmic law by u 
given by (2.17), we find that the relative importance of the friction e.m.f. 
to the dynamo field is given by »1/we.w/kpv2, which is smaller than 1 in our 
case. With the following conditions satisfied, the dynamo field term will pre- 
dominate over the other terms (plasma oscillation, ohmic electric, and Hall 
electric fields): 












(4.9) wK WV kpwiwe, wKky a , OK Kw. 


3 





In this case, the generalized ohmic law is reduced to 


(4.10) E+ uy A1=0, 







where u is the mean velocity of the whole gas. Therefore, in the frequency 
range 





(4.11) w <min E KpW1, Kp = , Vina | ; 
3 





we have Alfvén waves propagated with the velocity 







(4.12) 






which is smaller than V,; given by eq. (4.8) (cf. Fejer 1960). 





A(b) Weak Magnetic Field, we. K vy 
In the frequency range where 





wo > Kya" (<K ve) 
1 





1204 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


we have |u,| > |u,|, and thus, the term representing the friction e.m.f. pre- 


dominates over the dynamo field. 
When 


(4.13) w K kya. =, 
1 


we have u, = U, ~ U, and the dynamo field term predominates over the 
terms representing the plasma oscillation field, the ohmic electric field, the 
Hall electric field, and the friction e.m.f., if condition (4.9) and the following 


are satisfied: 


(4.14) w< on(®) 
Vi 


Then, we have Alfvén waves propagated with the velocity V2, in the frequency 


range 


; m ; We \- RaaeTe W We 
(4.15) wo <K min] «r(£) ; V kp@ We, is wal ; 
Vy) V3 


B. HIGHLY IONIZED GASES; Kk, = pr/p K1, pp = p 


B(a) Strong Magnetic Field, v; K we 

In this case, the friction force term and the first term in r.h.s. of the equation 
of motion of the plasma gas is always negligible compared with the ampére 
force term. Therefore, the equation of motion is reduced to eq. (4.2) and the 
conditions that the terms of the plasma oscillation, the ohmic electric field, 
and the Hall electric field are negligible compared with the dynamo electric 
field term are again given by (4.4) and (4.5). 

The importance of the friction e.m.f. relative to the dynamo electric field 
is conveniently discussed case by case. 

B(a)i. Very Strong Magnetic Field, v\/we K ty K 1 

In this case, the friction e.m.f. is always negligible compared with the 
dynamo electric field. In addition, it can be seen that u, = Uy, if w XK vo/kp, 
and u, > |uU,| if w > v2/K,. Therefore, we have the following conclusions: If 


Vo 6 WiWe 
<min| wy ‘| : 
Kn V3 


the Alfvén waves are propagated with the velocity V, in the frequency range 


. V2 , W {0 
(4.16) Kw «K min] oy, ‘| f 
Kn V3 


and with the velocity V»2 in the frequency range 


(4.17) w@ K vo/ Ky. 


If 
V2 , WiWe 
77 MIN] @4, » 
Kn V3 
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they are propagated with the velocity V2, in the frequency range 


(4.18) w < min| 0, ent | ; 

B(a)ii. ky K v1/ae K 1 

In this case, it can be seen that u, = U, if w XK v2.a,/v; and that |u,| > |u,| 
if w > vo(w,/r1). 

The friction e.m.f. is always negligible compared with the dynamo electric 
field if (v1/we)? K kp. When (1/we)? > ka, it is negligible only if w K ve(we/v1)?. 
Therefore, we have the following results. 

(a) Moderately Strong Magnetic Field, (v1/ we)? K ky K (vi/@e) K 1 

If 


® . @) 4), 
Vo—~ «min wo, = ‘| : 
V1 V3 


the Alfvén waves are propagated with the velocity V, in the frequency range 


(4.19) ve Kw min) ws, vee | 


Vi V3 
and with the velocity V, in the frequency range 


(4.20) w Ky_.—2 


Vi 


If 


Ge ° WiWe 
v2— > min) ws, ~ at 
Vv; V3 


the waves are propagated with the velocity V2, in the frequency range, 


(4.21) w <min| wo, vee | ; 


(8) Fairly Sirong Magnetic Field, ky & (v1/a.)?K «1 


We 


If 


. WiWe We 
min) w, | Kr—, 
V3 Vv) 


the Alfvén waves are propagated with the velocity lV’, in the frequency range 


aa W We 
@® QM @;, > : 
Vs 
wW ; ww 
"we < min] ws, ‘| 2 
Vi V3 


the Alfvén waves are propagated with the velocity V, in the frequency range 
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2 
(4.23) nMt<omin | 0, 2, (%) | 
Vi V3 


V1 


and with the velocity in the frequency range 


(4.24) ain. 
V) 


B(b) Weak Magnetic Field, wu K 

In the frequency range where w X v2(w,/v1), we have Uy = U, = U and we 
find that the relative importance of the friction e.m.f. to the dynamo electric 
field is of the order of (w/v2) . (v1/we)?. It is much smaller than 1 if w K vo(we/v)?. 
The conditions under which the terms representing the plasma oscillation 
field, the ohmic field, and the Hall field are negligible compared with the 
dynamo field are given by (4.9). Conditions in (4.9), however, are practically 
the same as in (4.4) and (4.5), since x, is almost equal to 1 in this case. There- 
fore, we have Alfvén waves propagated with the velocity V2, which is almost 
equal to V; in highly ionized gases, for the frequency range 


2 
or ° ® {W. ® 
(4.25) wo< min| os, ——, (ts) | ; 
V3 Vi 


If w > vo(w/v1), we have |u,| > |u,|, and we find that the friction e.m-f. 
term predominates over the dynamo field term. 
The frequency ranges where Alfvén waves can exist are shown schematically 


in Table I. 
5. RATE OF DISSIPATION OF WAVE ENERGY 
It is important to deduce the exact formula of the rate of dissipation of 
wave energy in problems of damping of waves and heating of ionized gases. 
Scalar multiplying both sides of equations (2.1), (2.2), and (2.3) by u,, uy, 
and u,, we can derive the following equation: 


0 ‘ 2 : 
~ Ot (Fpeuet 410i + dpnun) = —G-E)+Q, 
(5.2) Q = viem NV, (Uy —-U,) +venmN, (Ue —u,)’+v,m N,(u, - u,). 


On the other hand, one can also derive the following equation from Maxwell's 


equations: 


e ~ d BE’ -) ° 
(5.2 (i.E) = — es S 
).3 ) j E) aft + div 


where Poynting’s vector § is defined by 


(5.4 § = —E Ah. 


4n 


Then, we have 
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eee ee ee 5 Wadi 
at (2Pelet 301M i+ Fpnun) at \Se' 8a = divS+@Q. 


Equation (5.5) is the extended Poynting’s theorem. Q gives the accurate 
rate of dissipation of wave energy. Every term in Q involves a collision fre- 
quency in its coefficient, suggesting that the energy of the systematic motion 
of the gas particles is changed into the thermal energy through collisions. It 
is to be noticed that the quantity (j. E) does not give the rate of dissipation 
in our case. 

By making use of eqs. (2.13), (2.14), and (3.6), we haveT 


1/1 


: s DiI ie 
5.6 = =| —4.— jt t+— VS (j-7* 
(5.6) Q ee x) j VF (jzjt), 


where a1, o2, o2, x, and y are given as follows: 


] MeV ] M ie ] MeV 3 


o1 ea N,e* F i Npe*’ i Nye” F 

a Pnwi—p vi(me/m;) _ 2 Pn—P (Vi /we) 

(5.8) {= ee ae es 
p v2+ paw p v2+ pw 


and 
(5.9) 


j* represents the complex conjugate of j; hence 

(5.10) j.j* = frit t+iviv- 

In circularly polarized waves, we have 

(5.11) I (jejy) = tiejs = tiviy = +H-J*, 


where the + sign corresponds to the extraordinary waves, and thus, Q becomes 


— aft 3 
(5.12) Q = Dg 
-_ 3 2 


The energy of the magnetic field averaged over one period is given by 
(x /w*)( V?/c?)j.j*, with the condition that the propagation velocity of the wave 
is much smaller than the velocity of light. Under this condition, the amount 
of energy of the electric field becomes negligible compared with that of the 
magnetic field. In Alfvén waves, this condition is tacitly assumed to be satis- 
fied. Besides, in Alfvén waves the amount of the kinetic energy of the gas 
particles is equal to that of the magnetic field as mentioned in Section 1. 
Therefore, the density of the total wave energy averaged over one cycle is 
(24/w*)(V?/c*)j.j*. The equation of the extended Poynting’s theorem shows 
that dissipation of the wave energy becomes negligible if the amount of the 
energy dissipated in a time interval 1/w = (1/2) X(one period) is negligible 


tO is the dissipation rate averaged over one period. 
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compared with the amount of the total wave energy, which we have calcu- 
lated above. One can show that the condition is satisfied in the cases men- 
tioned in Section 4. 


6. APPLICATION TO THE IONOSPHERE 
The ionosphere is slightly ionized; thus the discussions in case A in Section 3 
are appropriate. 
Table II shows the values of w;, we, and vin, and vy, in the ionosphere 


TABLE II 


Ionospheric data due to Francis and Karplus (1960). #; = ion gyro-frequency, and we = elec- 

tron gyro-frequency at the geomagnetic latitude 45°. vin = the collision frequency of an ion 

with neutral particles. ven = the collision frequency of an electron with neutral particles. The 
meanings of 7, and 72 are given in Section 6 


Altitude, 


km Qi We Vin Ven T; T: 
80 1.60? 8.47° 2.055 3.298 3.9272 6.14-5 
90 1.60? 8.438 3.244 5.425 3.927? 3.8874 
100 1.62? 8.398 5.373 1.005 3.887? 2.34-3 
120 1.71? 8.32 3.12? 7.938 3.687? 4.037? 
140 1.76? 8.246 7.44! 2.66% 3.57-* 1.697% 
160 1.81? 8. 16° 2.83! 1.37 3.477 4.44-1 
180 1.87? 8.098 1.33? 8.86? 3.367 9.467! 
200 1.91? 8.026 6.94° 6.46? 3.20°* 1.81° 
250 2.047 7.84 2.01° 3.97? 3.087? 6.25° 
300 2.16? 7. 66° 8.257! 3.83? 2.917? 1.53! 
350 2.23? 7.498 3.917 4.20? 2.827? 3.21! 
400 2.31? 7.33° 1.94-! 4.43? 2.727? 6.48! 
450 2.34? 7. ke* 1.007! 4.32? 2.687? 1.267 
500 2.37? 7.01° 5.267? 4.14 2.657? 2.48? 
2.34? 6.86° 2.867? 3.73? 2.687? 4.407 


550 


(Francis and Karplus 1960). One can see that vy = yen and wm = vig. Vie iS at 
most of the order of 10?, and thus we also have, approximately, v3 ~ ven. 
Therefore, we have wya,/v3 = @@./Yen and the lower limit of this is of the 
order of 4X10? sec, which is larger than w, throughout the ionosphere. And 
thus, Alfvén waves propagated with the velocity V; can occur in the frequency 


range 
(6.1) vin Kw K wh, 


or in the period range 


(6.2) (Tye) cere tnt. 
@| Vin 

As shown in the table, vi, > w; (72 « 71), below an altitude of 100 km. 
Therefore, no Alfvén wave can appear below this altitude. The lower and 
upper limits in the period range are also shown in the table. Waves whose 
periods are shorter than about 3X10-? second cannot be of Alfvén (V}) 
type. The upper limit increases with altitude. It has been suggested that 
Alfvén waves can occur in the exosphere (Dungey 19545) or in the lower 
ionosphere (Watanabe 1959). For instance, a wave of period 1 second can 
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penetrate to an altitude of 180 km as an Alfvén wave. Below that it is trans- 
formed into an electromagnetic wave as in metallic substances (Watanabe 
1957; Piddington 1959). On the other hand, a wave of period 15 seconds 
can occur only above 300 km. 
The rate of dissipation of wave energy is simplified as follows in the iono- 
sphere: (m,./m;)v'n K veq and vp = v2, so that we have 
2 1 


w4 oe 29 
Vin tw 


WW 
ae 2% 
Vin tw 


a et) j -j* +40 en wry Iii)’ . 
in 


5 
Me Vin tw 


7. CONCLUSIONS 


In a slightly ionized gas, there appear two types of Alfvén waves if the 
external magnetic field is strong. One type, in which only charged particles 
partake of the wave propagation (V,), appears in a higher-frequency range. 
The other, in which neutral particles also partake of the wave propagation 
(V2), appears in a low-frequency limit. If the magnetic field is weak, only the 
latter type of low-frequency wave appears. 

In a highly ionized gas, the situation is more complicated, as shown in 
Table I. Practically speaking, however, there is little difference between the 
two kinds of propagation velocities, V; and V2. And, if we do not distinguish 
between them, the Alfvén waves can be always propagated in a frequency 


range limited by wo, 
(7.3 0< w<XK wo. 


The upper limit we is given by 


me . W 5H) 
(7.2) wy = min| oy | , 
< V3 


if the conditions «o vw, and (v;/w,.)? <«, are satisfied. In the other cases 


(fairly strong or weak magnetic fields), wo is given by 


\2 
. WjWe We 
oq = DUN @,, » V2 > 
. V3 Vy 
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ANGULAR CORRELATION OF BETA-GAMMA COINCIDENCES IN THE 
COMPTON EFFECT! 


F. W. VAN NAME, JR. AND JOHN W. Kocu 


INTRODUCTION 


In the theory of Compton scattering of photons it is well known that the 
relation between the angle @ of the scattered photon and the angle ® of the 
recoil electron is as follows: 


(1) cot 6 = (1+a) tan(6/2), 


where @ represents the energy of the incident photon in Mev divided by 0.511. 
Shankland (1937) verified equation (1) to within +20°. Cross and Ramsey 
(1950) verified equation (1) to within +1° using a single target and a single 
pair of angles. In view of the improved detection and coincidence circuits now 
available, it was decided to verify equation (1) for a variety of angles and with 


two different scattering targets. 


APPARATUS AND EXPERIMENTAL PROCEDURE 

The source was a 10-millicurie interstitial needle of Co® which was placed 
in a brass rod for ease in centering and for elimination of unwanted beta-rays. 
A hole } in. in diameter was drilled through three lead bricks, thus providing 
a collimating tube 5.5 in. long. The targets were beryllium, 0.028 in. thick, 
and carbon, 0.050 in. thick. Each target, in turn, was mounted at one end of a 
steel pipe 13.5 in. long, which had three Mylar windows each 0.001 in. thick. 
A vacuum of 15 microns was maintained in this pipe, so that recoil electrons 
would not be appreciably scattered in their passage from the target to the 
beta-counter. The large Mylar window behind the target allowed the escape of 
photons traveling at small angles to the initial direction. The arrangement of 
the apparatus is shown in Fig. 1. 

The recoil electrons were detected by an anthracene crystal and the scattered 
photons by a Nal(TI) crystal. The pulses from the two photomultipliers were 
amplified and sent to a coincidence circuit similar to that described by Garvin 
(1953), which has a resolving time of 4X10~-® second. The coincidences were 
recorded after passing through a scaler. 

The position of the gamma-counter was set at some particular angle and 
coincidences were measured for 1-hour periods for each position of the beta- 
counter. The beta-counter was moved in 2° steps up to 6° on each side of the 
expected angle of maximum coincidences. The data are shown in Table I. 

‘Abstracted from a thesis submitted by John W. Koch in partial fulfillment of the require- 
ments for the M.S. degree in physics at the University of Delaware. 
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Fic. 1. Diagram of the apparatus. 






TABLE I 












0 ® | 6 ® 













angle of angle of {| angle of angle of 
scattered recoil || scattered recoil 
photon electron photon electron 
Carbon | Beryllium 
100° 14°+2° 90° 16°+1° 
80° 18°+2° 70° 22°+1° 
65° 24°+2° 60° 27°+2° 
57° 28°+2° 53° 30°+2° 
46.5° 34°+2° 






















DISCUSSION OF RESULTS 
The two runs yielding the narrowest and widest peaks about the expected 
angle are shown in Fig. 2. The estimated uncertainties in locating the centers 
of these peaks are +1° and +2° respectively. A plot of the angle of the 
scattered photon, 6, versus the angle of maximum coincidences of the recoil 
electron, ®, is shown in Fig. 3. The estimated uncertainties in the values of 
are shown by vertical lines. Since the beryllium target was about half as thick 
as the carbon target, the beryllium target provided more accurate points. 
Although Co® emits gamma rays of energy 1.17 and 1.33 Mev, an average 
value of 1.25 Mev was used in computing the theoretical curve shown in Fig. 3. 
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Fic. 2. Coincidences plotted as a function of © for the narrowest and widest peaks. 
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Fic. 3. Plot of &, the recoil angle of the electron, versus 0, the angle of the scattered photon. 


Since the difference in energy of the two gamma rays should produce a change 
in ® slightly less than 2° on the average, the use of the average energy appears 
to be consistent with the accuracy of the measurements. All of the configura- 
tions used agreed with equation (1) within +2° or better for both targets. 





NOTES 1215 


Our thanks are due to Mr. Melvin Weaver for his help in the mechanical 
construction and to Mrs. John W. Koch for her help in taking measurements. 


Cross, W. G. and Ramsey, N. F. 1950. Phys. Rev. 80, 929. 
Garvin, R. L. 1953. Rev. Sci. Instr. 24, 618. 
SHANKLAND, R.S. 1937. Phys. Rev. 52, 414. 


RECEIVED May 15, 1961. 
DEPARTMENT OF PHyYSsICcs, 
UNIVERSITY OF DELAWARE, 
NEWARK, DELAWARE, U.S.A. 


FRANCK-CONDON FACTORS AND r-CENTROIDS FOR SOME BANDS OF 
THE SiO A'lI-X'=+ BAND SYSTEM 


A. T. McGrecor, R. W. NICHOLLS, AND W. R. JARMAIN 


A partial (v’+v"’ < 8) array of Franck—Condon factors gy for the SiO 
(A1II-X'!+) band system was computed to assist the interpretation of spectra 
of shock-excited powdered SiO». Private requests for Franck—Condon factor 
arrays of this and other band systems of diatomic inorganic oxides which we 
have studied prompts us to report Franck—Condon factors, r-centroids, and 
wavelengths for SiO (A‘II-X'Z*+) in Table I. 


TABLE I 


Franck—Condon factors qv'»", r-centroids 7,"», and wavelengths A,» for the SiO A'NI—X1zt 
band system 











a 


3 


0.17; 
1.66; 
2563.8 
0.09 
1.639 


LEGEND: Qp'v’’s Tv'v' (A), Av'n’(A). 
*Band wavelength calculated. 
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The wavelengths are mainly those of Jevons (1924); wavelengths of weak 
bands which he does not report have been calculated. The r-centroids were 
evaluated by methods discussed by Nicholls and Jarmain (1956) and the 
Franck—Condon factors were calculated by the method of Fraser and Jarmain 
(Fraser and Jarmain 1953; Jarmain and Fraser 1953; Fraser 1954) for Morse 
molecules using the basic data for the transition tabulated by Herzberg (1950). 

The band system involves a moderately large change in internuclear separa- 
tion (Ar, = 0.11 A) and exhibits primary and secondary Condon ‘parabolae’ 
in the Franck—Condon factor array. The entries for bands through which 
they pass are italicized. The eye estimate intensity distribution reported in 
Pearse and Gaydon’s (1950) compilation agrees well with the Franck—Condon 
factor array. 


The work has been supported in part by research grants from the National 
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of Scientific Research, and the Office of Naval Research. 
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THE HALF-LIFE OF Sm"3* 
R. E. GREEN AND W. H. WALKER 


During the course of some measurements to determine the resonance integral 
of Sm! (Walker and Green 1961) it was observed, when correcting the product 
activity, Sm'’, for decay, that best results were obtained using a half-life for 
Sm!3 which was smaller by 2% than the best published value of 47.1+0.1 
hours (Cork et al. 1958). Since we were concerned about the presence of 
impurities which would affect our resonance integral results we decided to 
investigate more fully the decay of Sm", 

The target material used consisted of a layer of Sm,03 (~140 ug/cm? of 
Sm!) baked onto 0.0025 cm thick superpure aluminum foil (see Walker 


*Issued as A.E.C.L. No. 1280. 
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and Green (1961) for details of foil preparation). The original sample material 
contained >98% Sm, the main impurities being Ce (<0.5%), Pr (<0.2%), 
and Eu (~0.1%). The samarium was highly enriched in Sm! (97.2%), so 
that radiations from the other samarium isotopes could be neglected. 

The foils were irradiated either in the ZEEP reactor or in the thermal 
column of the NRU reactor for an integrated flux of ~10° neutrons/cm?. After 
activation the targets were counted in a system which included two Nal(TI) 
scintillation counters and an automatic sample changer and readout unit. The 
counters were biased to accept all pulses corresponding to y-ray energies 
>30 kev. To check on counter stability, a source of long-lived Am*! 
(T1;2 = 470 years) was also counted throughout each experiment. The read- 
out unit recorded on punched tape the sample number, the counts in the 
two detector scalers, the counting time interval, and the time of day. The 
timing unit was a tuning fork oscillator with a frequency of 1000 c.p.s. This 
information was fed to a Burroughs-205 “Datatron’’ computer where each 
sample count was corrected for background and counter dead time (~3.0 usec). 

The counting was continued until the Sm'* activity had decayed to a few 
times background. The results of one run are shown in Fig. 1. Here the cor- 
rected activities of the Sm! and Am*! are plotted as a function of time. The 
constancy of the Am*! results indicates the stability of the system. In all, 
four separate measurements of the Sm'* decay were made. In each case more 
than 200 individual counts were obtained and in the least-squares fits 200-215 
points were used. In fitting the data it is important to correct for the presence 
of any contaminant activity. The main interfering activity is expected to 
come from the ~0.1% europium impurity since Eu! has a large thermal 
neutron capture cross section (~1400 barns) leading to a 9.35-hour isomeric 
state in Eu". In order to correct for this impurity the counting results were 
fitted by the method of least squares to a double exponential function of the 


form 
(1) A(t) = A,(O)es'+ A o(O)e2!, 


where A(t) is the corrected activity at time ¢, 
A,(0), 42(0) are the Sm! and Eu!” activities at zero time, respectively, 
andj, Az: are the corresponding decay constants. 

For the fitting A, was fixed at the value corresponding to the Eu'® half- 
life of 9.35++-0.02 hours obtained by Guében and Govaerts (1958), and 4,(0), 
A,(0), and A; were allowed to vary. The results of the least-squares fits are 
given in Table I for the four experiments. The table lists the Sm! half-lives 
obtained using both double (eq. (1)) and single exponential decay (.4.(0) = Q). 
The run labelled Sm'**/Am**! is a fit of run 4 with each Sm" activity point 
normalized to the nearest Am*! activity. This procedure should correct at 
least partially for any long-term drift in the counting system. 

The errors shown for the individual values of the half-lives are based on 
the ‘‘goodness-of-fit’”’ of the results to the function assumed, while the errors 
given for the mean values of 7,2 are based on the reproducibility of the 


individual results. 
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TABLE I 
Sm'3 half-life results from least-squares fitting 
















Number of Half-life, hours 
points Ba 
Run in fit Double exponential Single exponential 
















1 215 .14+0.09 45.97+0.06 
2 212 46.22+0.18 46 .33+0.13 
3 212 46.27+0.03 46.19+0.02 
4 200 46 .34+0.08 46.20+0.06 










Mean 46.24+0.04 46.17+0.08 


46.16+0.09 





200 










Sm53/Am*4! 








The value obtained for A2(0)/A (0) varied somewhat because of the different 
irradiation and counting conditions but was always <0.02, in agreement with 
that expected on the basis of the amount of europium impurity present. This 
fact was verified by a y-ray spectrum study made with a 100-channel kick- 
sorter; using y-ray spectra taken at intervals for several days it was found 
that all resolvable lines decayed with the Sm" half-life. It was concluded 
therefore that any contaminant activity present was initially <2% of the 
initial Sm" activity. The close agreement between the half-lives obtained 
using either single or double exponential decay functions supports this con- 
clusion. 












TABLE II 
Summary of Sm" half-life measurements 















Reference T1;2, hours 











Kurbatov et al. (1942) 47+1 
Miller and Curtiss (1946) 46 (no error quoted) 






Rutledge et al. (1952) 46.5+1 
Cork et al. (1958) 47.1+0.1 
Wiile and Fink (1958) 4543 






Guében and Govaerts (1958) 46.7+1.6 
Present work 46.2+0.1 

















In addition to counting statistics, the main sources of error in these measure- 
ments are: (1) the error in the absolute frequency of the tuning fork, (2) the 
variation of tuning fork frequency with temperature, (3) the uncertainty in 
the counter dead time, (4) non-statistical variations in the background, and 
(5) a long-term variation in the efficiency of the counting system. Errors (2), 
(4), and (5) have already been dealt with satisfactorily since they affect the 
reproducibility of the individual results. Also, error (2) is small since the 
frequency temperature coefficient is —0.015%/°C and the maximum tempera- 
ture variation in the air-conditioned counting room was 4° C. We have cali- 
brated the fork against a quartz crystal controlled oscillator and found it to 
be within 0.02% of the rated frequency of 1000 c.p.s. This will produce an 
uncertainty in 71/2 of +0.02%. The maximum dead-time correction applied 
in these experiments was ~0.1%; if we assume that the dead time is known 
only to +30% this will lead to a maximum error in 7,2 of +0.03%. 
















1220 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


As a result we take the total error in 71,;2 to be ~0.2% so that the value 
we obtain for the half-life of Sm'* is 46.2+0.1 hours. This result is ~2% 
less than that found by Cork et al. (1958), who also quoted about the same 
error. These and other values are listed in Table II. Except for the disagree- 
ment noted above, for which we have no explanation, the measurements all 
agree within experimental error. 
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CALCULATION OF THE vy, AND INVERSION ENERGY LEVELS FOR AMMONIA 
AND SOME OF ITS ISOTOPICALLY SUBSTITUTED SPECIES 


T. P. Norris AND J. M. DOWLING 


In 1935, Manning used a potential function with a double minimum pro- 
posed by Morse and Stueckelberg (1931) to calculate energy levels for 
ammonia (NH3;) and deuterated ammonia (ND;). The results were in good 
agreement with the available experimental data. 

Since then more and better data have been obtained from high resolution 
studies of the infrared (Plyler and Benedict 1957) and microwave (summarized 
by Townes and Schawlow 1955) spectra of these molecules. The purpose of 
this investigation was to refine Manning’s calculations in order to obtain 
better parameters for the potential function, and to calculate the energy 
levels for some isotopically substituted species, namely NH2D, ND:;H, and 
NT;3. 


The potential function used by Manning was 


(1) V(r) = 3 i 2 (241) 40] sech” iS —D sech* ey 


where p, 8, and D are parameters to be determined, r is the co-ordinate 


describing the motion, and 


(2) 
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where c is the velocity of light, 4 is Planck’s constant, and u is the reduced mass. 
The potential energy function has a double minimum at r = +ro, where 


(3) sech” &) = [2 (241)40]/ 2D. 


The Schrodinger wave equation has been solved (Morse and Stueckelberg 
1931) and the energy levels are given in terms of the roots of a continued 
fraction, which is a function of 8, D, and kp?. Since the potential energy function 
must remain invariant under isotopic substitution, it is seen that the pa- 
rameters 38(38+3), and D must increase or decrease in the same ratio the 
reduced mass increases or decreases. 

In the first set of calculations 79 was taken as the equilibrium value of the 
height of the NH; pyramid (0.3810 A, Plyler and Benedict 1957) and the angle 
between a N—H bond and the plane formed by the three hydrogen nuclei was 
taken as 22.1°. Essentially this reduced the number of parameters in the 
potential function from three to two by means of eq. (3). Computations 
were carried out to fit the inversion transition (0, <— 0,) and the higher transi- 
tions (3, <— 0,, and 3, —0,) for ammonia and deuterated ammonia. A wide 
range of parameters were investigated, but the calculated values were always 
in poor agreement with the experimental values. In particular it was found 
that when the inversion transitions were fitted closely, the calculated values 
for the higher transitions were considerably lower than the experimental 
values. 

Because of this difficulty calculations were carried out to find the three 
parameters 8, D, and kp? without using eq. (3). The best set of results is 
summarized in Table I and the parameters obtained are given in Table II. 
The over-all agreement is somewhat better than that obtained by Manning. 
In order to require eq. (3) to hold for 7) = 0.3810 A it was necessary to treat 
the reduced mass as unknown. The reduced mass is a function of the geometry 
and under these conditions it was found that the angle between a N—H 
bond and the plane of the three hydrogen nuclei is reduced by 8° from its 
equilibrium value. This is not unusual in view of the large amplitude of this 
vibration. The value of the reduced mass obtained in this manner should be 
considered as some sort of average reduced mass for this kind of vibration. 

The parameters from Table II give the following potential energy expression 


r 


y = 4 —111,482 sech’ { —— }+67,587 sech* = \ cm", 
2p 2p f 


which for ro = 0.3810 A yields p = 0.4262 A, anda barrier-hindering inversion 
of 2085 cm=. 

With all the parameters determined, the energy levels of NH2D, ND:H, and 
NT; were calculated. The results are summarized in Tables I and II. 

From this investigation it has become apparent that if better results are to 
be obtained, with this type of potential function, at least two modifications will 
have to be made. They are: (1) the reduced mass should be considered as a 
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NOTES 


TABLE II 


Parameters of the double minimum potential function for ammonia and some of its 
isotopically substituted species* 








Param- 
eter NH; NH.D ND-.H ND; NT3 


kp? =. 2.83554 X10 fF 3.63657 X10? = 4.37355 X10 ~=5.06937X10-? 6.93332 x10-? 
B 70.021 79.364 87 .082 93.792 109.772 
D 1916.5 2457.8 2955.9 3426.2 4686 .0 


*The parameters are expressed in units such that the potential ener5y function has the unit of cm. 
i number of significant figures is not justified but they were carried to secure internal consistency in the 
calculatious. 





variable, because of the large amplitude of vibration, and (2) the interaction 
with the symmetric N—H stretching mode »; should be included, since this 
vibration is mixed to a fair degree with the v2 mode. This may alleviate to a 
large extent the difficulty mentioned previously, namely the fitting of the 
upper- and lower-energy levels, since a term describing the N—H stretching 
vibration would increase the slope of the potential energy functions for 
r > ro and thus raise the upper levels. 
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EFFECT OF SOLUTE ATOMS ON GRAIN BOUNDARY MIGRATION 
IN PURE METALS 


E. L. HoLMEs AND W. C. WINEGARD 


INTRODUCTION 


It has been suggested by Liicke and Detert (1957) in a theory of grain 
boundary migration in the presence of solute atoms that, if the interaction 
between the impurity atoms and the boundary is sufficiently large, the rate 


Can. J. Phys. Vol. 39 (1961) 
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of boundary migration will be controlled by the speed at which solute atoms 
can diffuse along with the boundary. The measured activation energy in this 
case would be approximately equal to Qp, the activation energy for diffusion 
of solute atoms in the solvent matrix. Liicke and Detert also suggested that at 
high temperature or low concentrations, the boundary might “break away”’ 
from the solute atoms and the measured apparent activation energy would be 
greater than the value Qp. 

This theory did not appear to be supported by the investigations on the 
effect of solute atoms on the rate of migration of a single boundary by Aust 
and Rutter (1959, 1960), or by those on grain growth by Bolling and Winegard 
(1958) and Holmes and Winegard (1960). Recent work, however, using 
concentrations of solute higher than those used previously, may be interpreted, 
at least in part, by the theory of Liicke and Detert. 


THE EFFECT OF LEAD ON GRAIN GROWTH IN ZONE-REFINED TIN 
Grain growth experiments on alloys of lead in zone-refined tin have been 
performed using the technique of Bolling and Winegard (1958). Figure 1 
shows the results of these experiments and the activation energies obtained. 
The activation energy for grain growth increases, with increase in lead concen- 
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Fic. 1. Logarithmic rate of grain growth at 100 seconds plotted against 1000/7 °K for 
zone-refined tin with various concentrations of lead added. 
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tration, to a value of 18+2 kcal/g-atom before the effect of solute on the rate 
of growth close to the melting temperature becomes markedly evident. Further 
additions of solute do not affect the measured activation energy but do 
decrease the rate of grain growth at the melting temperature. The activation 
energies measured for the various concentrations are plotted in Fig. 2. 


Measured Activation Energy 
k cal /g-atom 


Po aan acta ca 
: 0.10 


Concentration, Atomic % Lead 


Fic. 2. Graph of activation energy plotted against concentration for grain growth in 
zone-refined tin with lead added. 
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Fic. 3. Schematic representation of graph of activation energy plotted against concentra- 
tion for (a) single boundary migration, (b) Liicke and Detert theory, and (c) a possible 
combination. 
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DISCUSSION 


Shown schematically in Fig. 3 are: (a) the results of Aust and Rutter on the 
effect of solute atoms on single boundary migration in zone-refined lead; 
(6) what might be expected from the theory of Liicke and Detert; and (c) a 
possible plot of activation energy versus concentration if both (a) and (6) 
are correct. At high concentrations of solute, the experimental results shown in 
Fig. 2 are in agreement with the theory of Liicke and Detert. At low concentra- 
tions, the value of the activation energy varies with concentration in a similar 
manner to that in the Aust and Rutter work, up to the value Qp. The dif- 
ference between Fig. 2 and Fig. 3c may be due to the fact that the grain growth 
experiments may not be sufficiently sensitive to determine the ‘‘breakaway”’. 


The authors gratefully acknowledge the support of the National Research 
Council of Canada (grant No. BT-48) and the University of Toronto. 
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A NOTE ON SOME PROPERTIES OF HAMILTONIANS INVOLVING 
SPIN MATRICES FOR HALF-INTEGRAL VALUES OF I 


G. M. VOLKOFF 


In discussions of problems of nuclear magnetic and quadrupole resonance 
a prominent role is played by the expressions for the magnetic dipole (Zeeman) 
interaction of the nuclear magnetic moment uw = yAI with a given uniform 
external magnetic field H = HP 


(1) ©, = —yH = —yhHP-l 


and for the electric quadrupole interaction of the nuclear electric quadrupole 
moment tensor eQ with a given external electric field gradient VE = —VV@ 


(2) Sq = —QO-VE = Ci[¢.2(312—1’) + (¢22— by) (12— T})] 
Cl(3—-1)+9(12-I5)] 


with 


— eQdz. _ 


- we *< Gage 2 —, 
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Here J, y, eQ are the nuclear spin, gyromagnetic ratio, and electric quadrupole 
moment, H is the magnitude, and P the polarization vector of the external 
magnetic field (constant or time dependent, and representing a linearly, 
circularly, or elliptically polarized field), x, y, z are the principal axes of the 
electric field gradient tensor $;; = 0°¢/dx,0x,; (defined by 12, = dy: = dz 
= 0), 7 is its asymmetry parameter, and I is the usual spin operator. 

Expressions (1) and (2) are simple examples of spin-Hamiltonians involving 
odd or even powers of the components J; of the spin operator I. Hamiltonians 
involving homogeneous polynomials in the J; of degree higher than the 
second can also be written down to express the interaction of higher nuclear 
moments with higher order inhomogeneities of given external electric or 
magnetic fields. However, in this note we shall be particularly concerned with 
(1) and (2). 

Since I? commutes with any polynomial in the J; the eigenvectors y of 
(1), (2), etc., can always be chosen as eigenvectors of I? describing nuclei of 
definite spin J. If T is the time-reversal operator (cf., for example, Heine 
1960, or Wigner 1959) it may be shown that for any such eigenvector 
(4) Ty = (—-1)"y. 

Any spin-Hamiltonian which is a polynomial with real coefficients containing 
only terms of even degree in the J;, in particular Hg defined by (2), may be 
shown to be invariant under time reversal. From this it follows that if wz, is 
one of its eigenvectors belonging to the eigenvalue E£,, then yy, = Tz, 
is also an eigenvector belonging to the same eigenvalue. It follows from (4) 
and from the invariance of the Hamiltonian under time reversal that when J 
is an odd half-integer (giving a minus sign in (4)) then py, = Tz, must be 
linearly independent of y,,, and that 7~y, = T°~z, = —Wz,- In other words, 
for half-integral J the eigenstates of (2) must necessarily be degenerate in 
pairs. This is an example of Kramers’ theorem (cf. Heine 1960; Wigner 1959). 

Therefore, it is known that for odd half-integral J the 2J+1 eigenvectors 
of (2) split up into two sets wz, and Yy,, of (1+43) eigenvectors each with one 
set being obtainable from the other set by time reversal. To demonstrate this 
explicitly, and to prepare the way for later theorems we have found it con- 
venient to introduce a special representation of the operators J; for half- 
integral values of J. 

The (27+1)-rowed Hermitian matrices representing the components of 
spin angular momentum I have their familiar form in the representation in 
which J, is diagonal and the order of the rows and columns of the matrices is 
chosen so as to exhibit the eigenvalues m of J, in the standard order 
m = I,I—1,..., —(IJ-—1), —I down the principal diagonal. In this repre- 
sentation non-zero elements of J, and J, appear only in positions one row or 
one column away from the principal diagonal. 

In discussing half-integral values of J we introduce a different representation 
which leaves I? and J, in diagonal form, and which is obtained by separating 
the 27+1 eigenvectors x4m of J, (and also of T°) into two equal sets L and M 
of I+4 elements each (which go over into each other by time reversal) by 
means of the definition 
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Mr, = —My, = [—2n, n=0,1,...(J-—}4). 


After rearranging the order of rows and columns so as to obtain a representa- 
tion in which the set of m, values is followed by the set of my values down the 
principal diagonal of J,, we can write 


(5) inte hate Bate, 


Here 3,, =,, 2, are (27+1)-rowed generalized Pauli matrices 


» nal), say es) 
1/0 . 0; -1/’ 


in which 1 stands for the (7+3)-dimensional identity matrix and the phase 
of =, has been chosen to agree with the ordinary Pauli matrices for J = 4 
and to make the traces of X, Y, Z all positive for all J. The 2; satisfy the 
generalized Pauli relations: 


‘7 2 2 9 i 0 
= e’%,, etc. 2, =>, = 2,= (; > ; 


X, Y, Z are ‘‘block-diagonal”’ matrices defined by 


p | e | e | 
| jell 2) Fe (-. 2) ” (419) 
8) - < ( I, ; ~ 0 ' i ; . ” 0 ¥ . 


in terms of the “‘half-matrices”’ /,, Fo I, which are (J+34)-rowed real symmetric 
matrices (Jf, diagonal), all with the same eigenvalues m, = J— 2n, 
n=0,1,... (J — 3), and with positive traces }(1+4), whose matrix elements 
(for m, m' both belonging to the my, set) are defined in terms of the matrix 
elements of the usual spin operators [; by: 


(9) Chileans = (fo ae ¢ ae = Cer. nuit lia incon! = Ora, 
Jn virtue of the properties of the J; matrices for half-integral values of J the 
matrix elements of the X, Y, Z matrices in the (lower) M subspace are identical 
with those in the (upper) Z subspace. As a result of this X¥, Y, Z commute 
with the =;. Examples of [,, [,, 7, are given for I = 3/2 and J = 5/2: 


je v3) at 0 i) = " 2% 
10) tn ( 9 ’ I ae —/3 2 ’ oe 0 of ’ Cm as 


0 V5 
3 —272 ff, 
; 0 


The commutation relations for the =; and the J; lead to: 


(1) 10,400, = (-1)7, ete, 24742 = 1(1+1)1. 


y 
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Expressed in terms of these matrices the nuclear quadrupole coupling 
Hamiltonian is given by 





(12) Sq = C[3Z?—I(I-+1)1+9(X?— ¥°)] = c(% *.) 


where the (J+3)-rowed real matrix 


(13) Ho = 3f2—1(1+1)i+n(f2- £2) 

has exactly the same form both in the L and the M subspaces. This is an explicit 
demonstration of the result embodied in Kramers’ theorem that for half- 
integral J the eigenvalues of Sg are all doubly degenerate, and that the pairs 
of degenerate eigenfunctions y, and yy in this representation consist of pairs 
of linear combinations (connected by time reversal) of class Z or class M 
eigenfunctions x,, of J, with exactly the same real coefficients. In other words 
the real unitary transformation D which diagonalizes Sg can be written in 
‘‘block-diagonal” form with the elements in the M subspace being identical 
with those in the Z subspace: 


A 


tai c(? | 9) (Bee (ee ) 4 (FS ) Z - D) 
) 0! D/\OIA/\ 0 D7 0 DHD-"/ \OE 


where £ is a (J+3)-rowed diagonal matrix with real diagonal elements 
E,(i = 1,2,..., 7+4) which in the cases considered so far have always been 
found to be distinct. 

The magnetic dipole coupling Hamiltonian (1) is given in terms of the 





representation (5) by: 

(15) Gz = —yhHPel = —yhH(P,X 2,+P,Y¥ 2,+P,Z =,). 

The transformation D which diagonalizes Hg, being in ‘‘block-diagonal”’ form 
with the same ‘‘half-matrix’’ D in both the Z and M subspaces, commutes 
with the =;, and therefore puts z in the form 


(16) D9;zD™ = —yhH(P,DXD"2,4+P,DYD" 3,+PDZD"  .). 
Here 
tip « 
™ xD = (PLE 0) te. 
(17) DXL 0 bib, etc., 


are still in “block-diagonal” form, and D/;D~' are still real symmetric matrices 
but with many more non-zero elements than the 7;. Whereas (12) does not 
involve the S;, and is therefore “‘block-diagonal”, both (15) and (16) do 
involve the Y; and are therefore, in general, neither “‘block-diagonal” nor 
“block-antidiagonal” in this representation. 

For perturbation theory calculations in which z given by (1) is treated as 
a small time-dependent perturbation on Sg defined by (2) it is useful to 
classify the eigenstates of Og into those which are connected by zero and 
by non-zero elements of DS,D~' as this determines the selection rules. Also 
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in connection with an investigation by Volkoff and Jog (to be published) of 
the generalized Bloch equations we were interested in the conditions under 
which D§ zD~ could be put into “block-diagonal’”’ or ‘‘block-antidiagonal”’ 
form (with zero matrix elements in two quadrants of the matrix, as in (6) 
or (8)) without upsetting the diagonal form of DO gD -'. Since the only trans- 
formations which will not upset the diagonal form of DSgD- are unitary 
transformations in the two-dimensional subspaces defined by each doubly 
degenerate eigenvalue E; we find it convenient to rearrange once again the 
rows and columns by putting them in the order L;, Mi, Lo, Mo,... Li, Mi,..., 
each pair L;, M; belonging to the appropriate eigenvalue E;. In this representa- 
tion the (27+1)-rowed matrices DOgD- and DSzD~ now take on the 
‘‘mosaic”’ form of (J+4)-rowed square arrays (1,7 = 1, 2,...I+4) of two- 
rowed matrices defined by 


(DSeD™");; 6,E,l, 
(18) (D&zD™)i, = —vhH[P.(BL,D7) 02+ (—1) 7 *P, (BLD) 0 
+P,(D1,D™) ,;0;] 


= —yhH (a jj0,+5ij0,+Ci;0;). 


Here (Df,D-),,; are one-dimensional matrix elements of f; between either 
the L,, L; or the M,;, M, states, 1, oz, o,, ¢, are two-rowed unit and ordinary 
Pauli matrices for J = 3, and a;;, bij, ¢4j are constant or time-dependent 
scalars. 

The transformations which leave DOgD-! diagonal are those unitary 
transformations in 27+1 dimensions, which in the “‘mosaic’’ representation 
consist of two-dimensional unitary matrices U; along the principal diagonal 
of the “‘mosaic’’ array, and zero elsewhere. Application of such a transforma- 
tion to DHzD~" yields for the two-dimensional 7,7 element of the ‘“‘mosaic”’: 


(19) U(D&2D™") yU5* = —yhH (aijU U5 +0 UU; +0 4jU U5") 
= —yhH[P,(D1,D™) U o2U7'+(—1)74P, (DL,D~ U 0,U7? 
+P,(D1f,D™") J e.U;" . 


In order for the whole D§zD~ to be either “‘block-diagonal” or ‘‘block- 
antidiagonal’’ each two-dimensional element of the “‘mosaic”’ representation 
must be either diagonal or antidiagonal. We consider the case when 
U; = U; = Us independent of 7. It is known that of the three Pauli matrices 
Or, Ty, 7, Only one at a time can be made diagonal, and only two at a time 
antidiagonal by a unitary transformation. This means that we can bring 
(DH zD-') ;; into diagonal form simultaneously for all values of 7,7 by the same 
transformation U if, and only if, any two of P,, P,, P, vanish (magnetic field 
along one of the principal axes of the electric field gradient tensor), and we 
can bring it into antidiagonal form if, and only if, at least one of P,, P,, P, 
vanishes (magnetic field in one of the principal planes, but not necessarily 
along a principal axis of the electric field gradient tensor). 
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If none of the three P,, P,, P, vanishes then, in general, we can diagonalize 
or antidiagonalize (DzD~");,; only in a particular 7,j position (and its trans- 
pose) in the ‘‘mosaic”’ array. In this two-dimensional subspace it has the form: 

~ e c a—1b 

(20) (D$zD day — a re ) 
with real values of 
(21) a@=P,(D1.D"),,, b= (-1)P,(DL,D—)y,  ¢ = P.(DI.D-),;. 
Diagonalization of (20) yields: 
c—r a-—1b 
a+ib —c—X 
The square of such a diagonal matrix element \? appears in the calculations 
by Volkoff and Lamarche (1954) of the intensity of nuclear magnetic and 
quadrupole resonance absorption lines in spodumene. The linear combinations 


¥* of the ¥z,, Yar; (andof the ~z,;, ~;) states corresponding to the eigen- 
values +d are: 


=NN--a?-0 =0 A= +(a?+b?+c’)}. 


(22) 








y? a * (KP yp teKO wy) 


(23) 
Ei 1 = bene 
yp‘ ) a A (e teK( Wi —-K vy) 
with 
(24) e=tan(b/a), K® = [l+c/(a?+b?+c?)4}}. 
yo = Ty, and y = —Ty™ are connected by time reversal. 


Antidiagonalization of (20) yields 
aio 


= 2 re ee 0 
(25) U(D$zD ) UV ‘= ~ tt eta 0 ) 
with 6, w arbitrary and 
y —_ es [yor+ i(d+e) = aloe ¢* Kw 4 “_) 
(26) —et( ge 8) gal, 
a —i(d+e) _» (+) (-) e™ -(+) 18, —i(8+6) 
Va = 7 le y-yv J= yl —e je Vi 
+-(K +e *K) vy). 


wv, and wz are connected by time reversal. 

(DS zD-");; can be diagonalized in all the 1,7 subspaces by the same trans- 
formation U defined by (23) only if e, K™ are all independent of 7,j, i.e., 
only if the ratios a/c, b/c are independent of i,j. As noted earlier, this happens 
if any two of P,, P,, P, are zero. This cannot happen if any two of the P; 
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(say, P;, P,) are different from zero, for this would imply that the correspond- 
ing I; (say, [,, f,) are multiples of each other, and this is not the case. 

(DS zD~")i; can be antidiagonalized in all the 7,7 subspaces by the same 
transformation U defined by (26) only if either ¢ is independent of 7,7 (which 
implies that either P, or P, is zero, since [, is not a multiple of f,) and we 
take 6 = w/2 and w = tan” (K©/K™), or P, = 0 (so that K* = 1) and 
we take 6 = 0,y = 0. 

Thus, for a nucleus of half-integral spin J in the general case a magnetic 
dipole interaction (15) can interconnect by non-zero matrix elements all 
eigenfunctions of the nuclear quadrupole interaction (12). But for special 
orientations of the magnetic field (confined to one of the principal planes of 
the electric field gradient tensor) linear combinations (26) of the doubly degen- 
erate eigenfunctions ~,;, Yu; of the quadrupole Hamiltonian (12) can always 
be found such that all the states can be divided into two equal classes 
Wva,, We, defined by (26) and having the property that the magnetic dipole 
interaction (15) has only zero matrix elements between states of the same 
(A or B) class, i.e. that §z is “‘block-antidiagonal” in this representation. If 
the magnetic field not only lies in a principal plane but is also parallel to one 
of the principal axes, a recombination of all the states into two other equal 
classes y{*? and ¥~ defined by (23) is always possible, which makes the dipole 
interaction (15) ‘‘block-diagonal’’, i.e. having zero elements between states 
of opposite class (+ and —). 


In the general case of a magnetic field oriented in an arbitrary manner with 
respect to the principal axes of the electric field gradient tensor, z given by 
(15) can always be diagonalized or antidiagonalized with respect to the eigen- 
states of Sg only in one particular two-dimensional subspace at a time defined 
by any two arbitrarily selected doubly degenerate levels E; and E, of the 
quadrupole Hamiltonian 5g given by (12). If 7 = 0 in (12) several, but not 
all, such subspaces of z can be diagonalized or antidiagonalized simul- 


taneously. 
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SOME REMARKS ON THE SPECTROSCOPIC STATE OF PARAMAGNETIC IONS 
IN DILUTE ALLOYS EXHIBITING RESISTIVE ANOMALIES 


E. W. Co.iinGs,* F. T. HEpGcock,t AND T. Sakupot 


Recent experimental investigations of the spectroscopic state of manganese 
in a weakly magnetic solvent reported by Collings and Hedgcock (1961) have 
led to the suggestion that the manganese ion is not in an S state as previously 
suggested by Owen et al. (1956, 1957). In fact appreciable crystal field splitting 
of the paramagnetic ions may be present. The immediate interest in the spectro- 
scopic state of paramagnetic ions in dilute metallic solution is the fact that 
these alloys exhibit anomalous electron transport properties at low tempera- 
tures (see Gerritsen 1959). The purpose of the present note is to report recent 
electron spin resonance experiments carried out on Au-Fe, Cu—Fe, and Cu-Co 
alloys, all of which exhibit so-called ‘resistive minima’’ at low temperatures. 
Table I contains the results of the present investigation where 4 /B is a measure 
of the symmetry of the resonance absorption curve and S/N is the signal-to- 
noise ratio of the absorption. The resistance ratios of the samples given in 
Table I indicate that the alloys contain appreciable paramagnetic ions in 
solution. The sensitivity and resolution of the electron spin resonance spectro- 
meter used is such that 5X10" spins could be detected in a resonance line of 
width over 300 oersted. Although both specimens of Cu-—Fe exhibited a 
resonance at room temperature the symmetry of the resonance absorption, 
the temperature dependence of the absorption, and the g value observed all 
indicate that the resonance was not due to paramagnetic iron in solution in 
copper but to precipitated agglomerates of ferromagnetic iron. It can be 
concluded therefore that no detectable resonance could be found for para- 
magnetic ions in solution in the dilute alloys studied. 

Included in Table I are values of the effective magneton number for the 
paramagnetic ions in the appropriate solvent. It can be seen that the effective 
magneton number is a marked function of concentration in most of these 
systems and seldom takes the limiting value corresponding to the magneton 
number calculated from the spin value of the isolated ion. In the previously 
mentioned publication of Collings and Hedgcock (1961), it was suggested that 
mixtures of spectroscopic states of the manganese ion would be consistent 
with the experimental results in the Mg—Mn system, where a mixture of 
states could result if either: 

(a) crystalline field splitting of the d state is sufficiently large to put the 
doubly degenerate state close to the conduction band Fermi level, or if the 
ionization energy for a 3d electron is too great to permit this, then 

(b) the exchange energy could place the ground level of the next filled state 

*Department of Physics, Victoria University of Wellington, Wellington, New Zealand. 

+Present address: Franklin Institute, Philadelphia 3, Pennsylvania. 


{National Research Council Pestdoctoral Fellow. On leave from Electrotechnical Laboratory, 
Nagatacho, Tokyo, Japan. 


Can. J. Phys. Vol. 39 (1961) 





TABLE I 


CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


Magnetic and electrical properties of some dilute copper and gold alloys. The data on Cu-Mn 


has been included for comparison purposes. 


The effective magneton number is calculated 


from magnetic susceptibility data and comes from the following sources. Au-Fe, Krongvist 
(1952); Cu-Fe, Bitter et al. (1950), Hedgcock (1956); Cu—Co, Jacobs and Schmitt (1959); 
Cu-Mn, Owen et al. (1956, 1957) 





Alloy 
Atomic per cent 


P (Qcm) 


Resistance ratio 
Ry 2/Re.73 

No. of paramag- 
netic ions in 


r-f. field 


Electron spin 
resonance 


Px (units of Bohr 
magneton) 


No. of d electrons 


Cu-Mn 
1.14 Mn 


3.4X 107° 


2.6X 10" 


A/B ~2.6 
~ 2.0 
H ~ 330 oe 
S/N ~ 20 
at 77° K 


Au-Fe 
7.06 Fe 


1.76 X10-5 


0.847 


1.7X 10" 


No resonance 
at 293° K 

at 77°K 

at 4.2°K 


~ 3.9 
at low conc. 


at high conc. 

~ 3/2 at low 
conc. 

~ 2 at high 
conc. 

3 or 7 at low 


conc. 
4 or 6 at high 


Cu-Fe 
4.55 and 
0.55 Fe 


4.03 X 10~* 


0.59 


7.7X10" 


A/B~1.0 

g ~ 2.16-2.07 

AH ~ 250 to 
170 oe 

S/N ~ 6-5 

at 293° K 

no signal at 
lower temp. 


~1.9 

at low conc. 
~4.9 

at high conc. 


3-1 at low conc. 


~ 2at high 
conc. 


2 or 8 at low 


conc. 
4 or 6 at high 


Cu-Co 
0.57 Co 


4.38 X10 


0.575 


1.0X10" 


No resonance 
at 293° K 
at 77° K 


~2 

at low conc. 
~ 4,2 

at high conc. 


-1 at low conc. 


~ 2 at high 
conc. 


2 or 8 at low 
conc. 
4 or 6 at high 


conc. 


conc. 


conc. 








near the Fermi level of the conduction band. If a similar mechanism is suggested 
for the dilute alloys considered in Table I, then the following mixtures of 
states could be suggested: 
Fe in Cu or Au, (3d)° and (8d)§ (a), 
or (3d)§ and (3d)? (d), 
(3d)? and (3d)° (a), 
or (3d)? and (3d)8 (Dd), 
(3d)* and (3d)5 (a), 
(3d)5 and (3d)® (0). 


It should be noted that none of the suggested mixtures of states results in the 
presence of an S state for Fe in Cu or Au or for Co in Cu. However, either 
mechanism (a) or (b) results in the presence of some manganese ions in a 
(3d)> state when dissolved in either Cu or Mg. The fact that there are no ions 
in a (3d)* state could result in the non observance of a resonance in the alloys 
Cu-Fe, Au—Fe, and Cu-Co. 


Co in Cu, 


Mn in Cu, 
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An interesting speculation exists that if process (6) is chosen then the 
mixture of ions will always result in the presence of some ions in the (3d)® or 
(3d)? configuration. Due to the Jahn-Teller effect both of these ions show zero- 
field splitting when present in insulating solids. If a splitting of a similar nature 
persists in a conducting solid, then the ground state degeneracy is removed in 
zero applied field and the resulting levels would be available for the Schmitt 
(1956) scattering mechanism to operate at low temperatures. 
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DIRECTIONAL CORRELATION OF THE 769-609 KEV 
GAMMA-RAY CASCADE IN Po?!4 


H. W. TAyLor AND R. McPHERSON 


The level structure of ssPoi3) (RaC’) excited by the beta-decay of s3Bij31 
(RaC) has been the subject of many experimental studies (Mayer-Béricke 
1959). As many as 58 gamma rays have been reported in the literature. In 
addition the conversion coefficients and K/L ratios have been measured for 
a number of the most intense transitions, and spins and parities assigned to 
many of the levels involved. A consequence of this work has been the identifi- 
cation of some levels as collective excited states and others as single particle 
excitations. Since Po‘ is an even-even nucleus only slightly displaced in Z 
and N from doubly magic Pb*®®, the collective levels are expected to arise 
from simple phonon excitations of a spherical nuclear shape. 

The first excited state of Po?“ at 609 kev hasa spin 2+, and the multipolarity 
of the associated ground state gamma transition is E2 (Mladjenovic and Slatis 
1954). A0-+ level at 1416 kev has been identified with certainty from a-particle 
measurements, internal conversion work, and directional correlation measure- 
ments with the 807-609 kev cascade which de-excites this level. If this level 
is part of the expected collective triplet state with spins 0+, 2+, 4+ the other 
members of the set are less firmly established. Levels at 1281 (Van Patter 
1959) and 1378 kev (Alder et al. 1956; Bishop 1958) have been identified as 
the second 2+ state in the decay scheme. The second of these appears the more 
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likely since the existence of a 769-609 kev cascade and a 1378-kev crossover 
transition have been firmly established. Measurement of the K/L ratio by 
Mladjenovic and Slatis (1954) (K/Z = 3.9) indicates a pure E2 multipolarity 
for the 769-kev gamma ray. The measurement of Nielson, Nielson, and Wag- 
goner (1956) (K/L = 4.6+1.2) would permit a considerable M1 mixture in 
the transition. The present directional correlation study was undertaken to 
check the 2+ spin assignment for the 1378-kev level and to determine the 
mixing ratio (if any) for the 769-kev gamma ray. This information would also 
extend our knowledge of 2-2-0 spin sequences near magic nuclei. 

The coincidence spectrometer used for the experiment has been described 
in some detail elsewhere (Taylor and McPherson 1960a). About 30 microcuries 
of RaCl. in HCI in equilibrium with its decay products was used as source 
material. All the gamma rays with energies greater than 500 kev emitted by 
such a source belong to the Bi*'* decay. The material was confined to a 
cylindrical cavity (2mm diameterX2 mm high) in a cylindrical piece of 
plastic with 5 mm thick walls. 

The pulse-height spectrum of the gamma radiations in coincidence with the 
609-kev ground state transition from the first excited state of Po?!4 is shown 
in Fig. 1. The data were obtained with two 3.75 cm diameter X2.5 cm Nal 


;—— 769 
° 


COUNTS PER MINUTE 








PULSE HEIGHT 


Fic. 1. Gamma-ray spectrum seen in coincidence with the 609-kev ground state transition. 
Crystal size: 3.75 cm diameter X2.5 cm long; source-to-crystal distance 5 cm; angle between 
counter axes 90°. 
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crystals located 5 centimeters from the source with an angle of 90° between 
their axes. The two gate positions A and B were used for the moveable counter 
in the directional correlation measurements. Gate A was chosen to include 
most of the 769-kev photopeak and B to include most of the 807-kev photo- 
peak. By analyzing the spectrum using standard line shapes and the energies 
of the more intense gamma rays as reported in the literature, the contributions 
to the counting rates of gates A and B due to the gamma rays involved were 
determined. The Compton continuum of gamma rays with energies above 
935 kev remains relatively flat through these two gate positions. 

For the directional correlation work, the two counters were placed at the 
afore-mentioned source-to-crystal distance in a horizontal plane with the 
counter axes passing through the source center. The dependence of the singles 
rate of the moveable counter on the angle @ between the counter axes was 
constant to better than 1%. Two experiments were performed: (1) the gate 
of the fixed counter was positioned on the central portion of the 609-kev 
photopeak, the gate of the moveable counter at position A; the coincidence 
rate was measured every 10° over the interval 90° < @ < 180°; (2) the gate 
of the moveable counter was located at position B and the experiment repeated. 
More than 10° coincidences, of which about 2% were accidental, were recorded 
in each of the experiments. The coincidence rates for both experiments were 
least squares fitted to the correlation function W(@) = 1+A2P2(cos @)+A4P, 
(cos 6). To obtain the directional correlation of the 769-609 kev cascade, 
the correlation of the Compton continuum was subtracted from the result 
obtained with gate position A. The continuum correlation was calculated 
by correcting the results obtained with gate position B for the presence of 
the correlations of the 935-609 kev (3-2-0) and 807-609 kev (0-2-0) 
cascades. The results for gate A were also corrected for small contributions 
from these cascades. After a further correction for the finite angular resolution 
of the counters, the correlation coefficients for the 769-609 kev cascade 
were found to be 42 = —0.220+0.065, A, = +0.378+0.125. Most of the 
error in each coefficient arises from the uncertainty of the spectral analysis. 
Additional contributions come from the counting rates and the experimental 
errors of the directional correlation results of Bishop (1958). 

A plot of Aq vs. Az using the mixing ratio 6 of the 769-kev transition as a 
parameter is shown in Fig. 2 for the spin sequences 4(2,3)2(2)0, 3(1,2)2(2)0, 
2(1,2)2(2)0, 1(1,2)2(2)0, and 0(2)2(2)0 (Taylor and McPherson 19605). The 
experimental result is shown as a point with the appropriate error bars. The 
6 = © point gives the coefficients appropriate to a pure quadrupole 769-kev 
transition. Assuming that the pure £2 character of the 609-kev transition 
is correct, a spin assignment of 2 for the 1378-kev level would appear to be 
unambiguous. The value of 6 determined from the value of Ae only is 
—4.3t1 8, corresponding to a mixture of 95% quadrupole and 5% dipole 
radiation for the 769-kev gamma ray. 

The information concerning the 1378-kev level may now be summarized: 
the energy ratio E2/E, for the second and first 2+ levels is 2.26, the intensity 
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0(2)2(2)0 


2(1,2)2(2)0 


1(1,292(20 


Fic. 2. Plot of A, against A» using the mixing ratio 6 as a parameter for the spin sequences 
shown. 


of the 1378-kev crossover transition relative to the 769-kev stopover tran- 
sition is about unity, the spin of the level is certainly 2+ with the multipole 
order of the 769-kev gamma ray being predominantly E2. These data would 
appear to confirm the suggestion that this level is the second 2+ level of 
Po?!4, 

The nature of the level at 1281 kev is left somewhat uncertain by this 
result. It could be another component of the predicted triplet state in which 
case its spin would be 4+, or a single particle state. Bishop (1958) has sug- 
gested 0+ and Mayer-Béricke (1959) 2+ as the spin of the level but these 
assignments are quite uncertain as yet. 


The work has been supported by the Atomic Energy Control Board, the 
National Research Council of Canada, and Queen’s University. 
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- LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


The Identification of Vibrational Levels in H; 


Using an electrostatic electron selector (Clarke 1954; Marmet and Kerwin 1960) we have 
measured the energy intervals between several of the vibrational levels in Hz (Marmet and 
Kerwin 1960b). It has recently been pointed out by Stevenson (1960) that difficulties in the 
interpretation of appearance potential curves are implied by our results. To be complete, 
these required the establishment of the ionization potential of Hz as an additive constant to 
the vibrational intervals. It was not possible to obtain this in our case without mass analysis. 

A selector has now been installed in the mass spectrometer at Saint Francis Xavier Uni- 
versity, and further data have been obtained from the Laval instrument. The ionization 
potential of Hz has been obtained by comparison with He, and has been measured as 15.37 +0.05 
ev. This is the additive constant required to establish that the first vibrational interval measured 
corresponded indeed to the difference between the v’ = 0 and v’ = 1 levels of the 22, state of 
H;. Two other intervals have been measured, and our results are now as follows: 





Vibrational interval (vj —v}), ev +0.01 





Experiment - 4 2-% 3-4 

No mass analysis . 2 .230 

Mass analysis f . 233 0.245 

Average 22. ; . 233 . 237 

Reference 2: 223 
(Richardson 1934) 


The relative probabilities of ionization from the molecular ground state to the various 
vibrational levels of the ionic state may be calculated from the slopes of the various parts of 
the ionization curve. The probabilities are, for succeeding states: 0.45, 0.77, 1.00, 0.55, 0.35. 
Krauss and Kropf (1957) have given theoretical values for these probabilities as: 0.44, 0.87, 
1.00, 0.91, 0.72. Our low values for the higher states may be due to the unfavorable geometric 
conditions under which they are measured (Marmet 1960) and we hope to conduct experiments 
to investigate this point. 

All of these data appear to constitute evidence favoring the Franck—Condon principle 
as applied to the hydrogen molecule. The experimental evidence offered by Stevenson to the 
contrary seems to suffer from two ills: first, as Stevenson points out, the interpretation of 
ionization curves is open to question, the more so because of the recent observance of structure 
in the curves that have been traditionally used as references (Fineman and Bouffard 1960; 
Kerwin, Marmet, and Clarke 1961); secondly, there is some question in our minds as to 
whether the potential curve for the ?Z, state of Hz used by him may be in slight error, and in 
fact should not be a trifle lower, as required by his data. 
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LETTERS TO THE EDITOR 


Lifetimes of the Low-lying Levels in Net 


We have measured the lifetimes of the 1.63-, 4.25-, and 4.97-Mev levels in Ne” by the 
Doppler shift attenuation method (Devons et al. 1955; Litherland et al. 1961). 

Triply charged carbon ions of energies near 18 Mev produced by the Chalk River tandem 
accelerator were directed upon unbacked carbon targets and carbon targets evaporated on 
thin magnesium, aluminum, and copper foils. The energies of gamma rays emitted from 
Ne** following the C(C!,.a)Ne®* reaction were observed in gamma-alpha coincidences 
between a movable 5 in. X6 in. Nal crystal gamma-ray detector and a p-n silicon junction 
alpha-particle detector placed at zero degrees to the ion beam (Gove ef al. 1961a). 

Since the neon ions associated with the coincidences were moving forward with velocities 
near 4X 108 cm/sec, Doppler shifts in the gamma-ray energies would be expected with maximum 
values for the unbacked target of about v/c ~ 1.3%. The gamma-ray energy as a function 
of angle is given by E = Eo (1+ F(v/c)cos 9) where F is an attenuation factor which takes 
into account the slowing down of the neon ions before gamma-ray emission. 

Measurements of dE/dx for heavy ions in various materials made at Manchester (K. 
Ramavataram, private communication) were used to derive slowing down times for the target 
backings. For neon ions of initial velocity near 4108 cm/sec, the velocity decreases almost 
exponentially with time in many materials. Consequently the attenuation F of the Doppler 
shifts is given simply by F = a/(a+r) where a is the mean slowing down time and r is the 
mean life of the gamma-ray emitting system (Devons et al. 1955). 

Coincidence gamma-ray spectra (Gove et al. 1961a, Fig. 1(b)) were measured at a number 
of angles in order to determine the small shifts of energy with angle. Gain changes in the 
gamma-ray detection system were monitored by simultaneous observation of the 1.28-Mev 
gamma rays from a Na® source. Least-squares fits were made to the gamma-ray energies as 
a function of cos @ in order to determine F(v/c) for both the backed and unbacked targets 
(see Fig. 1 for an example of these results). The experimental value of F for a given material 


DOPPLER SHIFTS 
IN ALUMINUM 


0 
cos @ 


Fic. 1. Doppler shifts in gamma-ray energy as a function of the angle of emission, 9, relative to the direction 
of motion of the emitting nucleus, which travels in aluminum. The shifts associated with pc spuld ation of the 4.25-Mev 
level in Ne® are shown. Note that the shift is smaller (though finite) for the second gamma-ray in the cascade, 
which shows directly that the 1.63-Mev level has a lifetime near the slowing down time for neon ions in aluminum. 


was obtained from the ratio of two such results, one for the unbacked and one for the backed 
target. Figure 2 shows the determination of the lifetime of the 4.97-Mev level from the value 
of F measured in magnesium. The errors in the lifetimes are asymmetrical as one can see from 
this figure, and in particular the result for the 4.97-Mev level is four standard deviations away 
from zero lifetime. 

In Table I are listed the lifetimes obtained for the 1.63-Mev, 4.25-Mev, and 4.97-Mev 
levels in Ne”, The slowing down time, a, is given for each material used. 

The lifetime value 5.67373 X 107 sec, obtained for the 1.63-Mev level, is in agreement with 
the previous work of Devons e al. (1955) and Lemberg ef al. (1950). 


tissued as A.E.C.L. No. 1282. 
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NEON IONS IN MAGNESIUM 
a=8-9xl0” sec 


oO 
io” io" oC 10" 10° 
T (sec) 


Fic. 2. The relation F = a/(a+r) is plotted against 7 for neon ions in magnesium. Also shown is the measured 
value of F for the 4.97-Mev level in Ne. The band encloses one standard deviation in either direction. 


TABLE I 








Level (Mev) i 4.25 4.97 
E,(Mev) : 2.62 3.34 
Slowing materials Al, Cu Mg, Al 
a(10-% sec) ; is 5.9,2.9 8.9, 5.9 


Texp (10-# sec) 613. 0.76+-72 19 +35 
ifn multipolarity E2 2 El, M2 
M\? g 


M1, E2 ay 
26 4.4X 10-4, 1.6X 10-3, 1.6 X 10-5, 3.3 K 10? 





Table I also gives, for certain assumed multipolarities, the ratio ||? of the measured 
gamma-ray transition probabilities to the Weisskopf estimate which is equal to the ratio 
B(E2)|/B(E2)|s.p. for the E2 transitions. The values of |M|? for the transitions from the 
1.63-Mev and 4.25-Mev levels, which have spin and parity values 2+, 4+ respectively 
(Broude and Gove 1961; Gove et al. 1961a), are seen to be typical of those for known E2 
transitions in light nuclei (Wilkinson 1960). 

The most striking result is the long (19*33 10-8 sec) lifetime of the 4.97-Mev level. The 
3.34-Mev gamma ray from this level has a dipole-to-quadrupole amplitude ratio equal to 
about 0.1 (Broude and Gove 1961; Gove et al. 1961a). From Table I it appears that the values 
of |M |? for this transition are reasonable (Wilkinson 1960) except that for the E2 component 
which is present if the 4.97-Mev level has positive parity. The small value of |M|?, 310-3, 
would make this the most retarded E2 transition yet observed in light nuclei despite the 
presence of two strongly enhanced £2 transitions in the same nucleus. This result for the 
4.97-Mev level suggested a measurement of the polarization of the 3.34-Mev gamma ray in 
the Ne” (pp’y) reaction. This measurement has been made, with the result that the 4.97-Mev 
level has been found to have negative parity (Gove et al. 1961b). This latter result has been 
confirmed in alpha-particle angular distribution measurements by Almqvist et al. (1961), and 
in alpha-gamma angular correlation measurements (Gove ef al. 1961a). 
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LETTERS TO THE EDITOR 


Alpha-Gamma Angular Correlations in the Reaction C(C",a)Ne* 


Angular correlations between gamma rays and alpha particles detected at 0° have been 
measured in the reaction C!(C?,a)Ne™ using beams of 17- and 18-Mev triply charged carbon 
ions from the Chalk River tandem accelerator. The geometry employed was the same as that 
for the lifetime measurements (Clark, Gove, and Litherland 1961) and is shown in an inset 
on Fig. 2. The detectors employed were p-—n silicon junction counters from 5 to 15 mm? in 
area and of resistivity 10,000 ohm cm,f located about 1 centimeter from the target. Two 
Nal(T1) detectors were employed, each 5 in. in diameter and 6 in. long, 6.25 inches from the 
target. Sufficient aluminum absorber was inserted between target and detector to stop the 
carbon ions. One gamma-ray detector was moved in angle with respect to the beam direction 
and the other was fixed at 90° to the beam and served as a monitor. The spectra of gamma 
rays from each detector, in coincidence with a particular alpha group, were displayed on two 
100-channel pulse-height analyzers. 

The upper half of Fig. 1 shows the alpha-particle spectrum measured directly by the zero- 
degree alpha counter. The inset level diagram shows the levels in Ne”; the present experiment 


c?+C'*_-Ne+a, 
18 Mev C” IONS 
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28 , Mev STANDARD 


3°34 Mev 


COUNTS PER CHANNEL 


a6 
CHANNEL NUMBER 


Fic. 1. (a) Alpha-particle groups from the reaction C!2(C!2,a) Ne as observed in silicon p—m junction detector 
at zero degrees to the carbon ion beam. Gates were set on individual groups for the coincidence measurements. 
A partial level scheme for Ne” is also shown. 

(6) The gamma-ray spectrum in coincidence with the alpha-particle group a3 leading to the 4.97-Mev level in 
Ne®, The 1.28-Mev gamma ray from Na” which was used for calibration in lifetime measurements is also present 


in the spectrum. 
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involved measuring the angular distributions of 3.34- and 1.63-Mev gamma rays in coincidence 
with the alpha group a@;, the 2.62- and 1.63-Mev gamma rays in coincidence with a2, and the 
1.63-Mev gamma rays in coincidence with a;. A coincidence spectrum for the first of the above 
cases is shown in the lower half of Fig. 1. In addition to the 1.63- and 3.34-Mev gamma rays a 
1.28-Mev peak from Na” is shown which was employed in lifetime measurements. 

Because the incident particles (C!*), the target (C), and the outgoing particles (He‘*) all 
have zero spin and positive parity the coincident gamma-ray angular correlations measured 
with the alpha detector at 0° are particularly easy to interpret (Litherland and McCallum 1960; 
Litherland and Ferguson 1961) since only the m = 0 substate of the residual nucleus is popu- 
lated. However, if the state of the residual nucleus is 0— it will not be fed in this reaction and 
if it has spin and parity 1+, 2—, 3+, 4—, etc. the alpha particles cannot be emitted at zero 
degrees (Litherland 1961). This latter condition is modified because of the finite size of the 
alpha detector and, in these particular cases, only the m = +1 magnetic substates of the 
level in the residual nucleus are populated. This provides a method for measuring not only the 
spin of a level but also its parity (Litherland 1961). 

The spin of both the 4.97-Mev and the 4.25-Mev levels in Ne” as well as the dipole-quadru- 
pole mixing in the 3.34-Mev radiation had been previously measured (Broude and Gove 1961) 
and the present work was undertaken to check these assignments. The results are summarized 
in Table I where values of a2/ao and a;/ao obtained both experimentally and theoretically are 


TABLE I 


Experimental and theoretical coefficients in the angular correlation of gamma rays 
in coincidence with alpha particles measured at 0° in the C!2(C!?,a) Ne” reaction fitted 
to W(@) = 1+(a2/ao0) P2(cos O-rles ao) Ps(cos 8) 











Level y-ray 
(Mev) (Mev) a2/a0 as/ao 
4.97 1.63 0.23+0.04 —0.77+0.06 Experiment 
0.18 —0.76 Theory 
3.34 0.15+0.04 —0.09 +0.06 Experiment 
0.19 0 Theory 
4.25 1.63 0.39+0.04 —0.45+0.05 Experiment 
0.51 —0.37 Theory 
2.62 0.59 +0.02 —0.62+0.03 Experiment 
0.51 —0.37 Theory 
1.63 1.63 0.58+0.05 —1.40+0.07 Experiment 
71 


0. - “1. 71 Theory 


compared. The experimental points were fitted by a least-squares computer program to a 
Legendre polynomial expansion and the theoretical coefficients in Table I have been corrected 
for the finite size of the gamma-ray detector. In the case of the 3.34-Mev radiation the quadru- 
pole-dipole amplitude ratio of +0.08 obtained by Broude and Gove (1961) was employed. 

The results for 1.63-Mev gamma rays measured in coincidence with alpha particles feeding 
the 4.97-Mev level in Ne” are shown in Fig. 2. The results of this correlation and of the 3.34- 
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Fic. 2. The angular correlation of 1.63-Mev gamma rays in coincidence with alpha particles, feeding the 
4.97-Mev level in Ne, observed at zero degrees as shown in the inset diagram. The dashed and solid curves are 
the theoretical ones for Jr = 2+ and 2— respectively for the level. Note that in this coincidence measurement 
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Mev gamma correlation not only confirm the spin assignment and multipole mixing (assuming 
no E3 radiation) of Broude and Gove (1961) but also establish negative parity for the level. 
The parity has also been established by measurements of the linear polarization of the 3.34-Mev 
gamma ray (Gove et al. 1961) and by direct measurement of the alpha-particle angular dis- 
tributions (Almqvist et al. 1961). 

Although the agreement with theory is not as good the results for the 4.25- and 1.63-Mev 
levels establish unambiguously spin and parity values of 4+ and 2+ respectively for these 
levels. The errors shown on the experimental coefficients do not include any systematic errors 
and are therefore probably underestimates. 

When the gate on the alpha counter was set on the alpha group feeding the 5.63-Mev level 
a weak 4.00-Mev gamma ray was observed in the coincident spectrum. From this measure- 
ment a value of I'y/f! = .04+.02 was obtained and proved that the 5.63 level decays by 
alpha emission. Recent measurements of this ratio give a value of 0.07+0.01 (Kuehner and 


Almqvist 1961). 


ALMOQVIST, E., BROMLEY, D. A., KUEHNER, J. A., and WHALEN, B. 1961. ‘To be published. 
Broupe, C. and Gove, H. E. 1961. Bull. Am. Phys. Soc. Ser. II, 6, 37. 
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The Interpretation of Angular Distributions and Angular Correlations 
from the Reaction C!?(C”,ay)Ne°t 


In recent studies of the C!2(C!2,ay)Ne”® reaction (Gove et al. 1961) it was found that the 
angular correlations of the gamma rays in coincidence with alpha particles, detected in a 
small axially symmetrical counter at zero degrees to the beam of C® ions, agreed with the 
theory given by Litherland and Ferguson (1961) with one exception. It was found that the 
angular correlation of the 1.63-Mev gamma ray following the 3.34-Mev unobserved cascade 
gamma ray from the 4.97-Mev state in Ne®® was in complete disagreement with theory. This 
disagreement has now been explained and has led to an increased understanding of the 
angular distributions of the alpha particles from the C!%(C!,a)Ne*® and similar reactions. 

In the paper by Litherland and Ferguson (1961, hereinafter referred to as LF) it was 
proved that in reactions of the type X(4;,h2)¥*, in which the particle hz is detected in a small 
counter at 0° or 180° to the beam of h; particles, the magnetic quantum numbers of the sub- 
states of Y* do not exceed the sum of the spins of X, 4, and he. In the case of the C(C",a) Ne”? 
reaction one might expect that only the y = 0 magnetic substate, of a state of spin ¢ in the 
residual Ne®°, would be populated. However, the element of the density matrix corresponding 
to the y = 0 magnetic substate, given by equations (21), (28), and (29) of LF, is identically 
zero if the parity of the state of spin ¢ is (—1)**!. This result is a consequence of vanishing 
Clebsch-Gordan coefficientst and is independent of the size of the axially symmetrical counter 
detecting the alpha particle and the complexity of the situation in the compound nucleus. 

Two important conclusions relating to angular correlations and angular distributions follow 
from this result. 

(1) In the case of the C!2(C!?,ay) Ne®® reaction and similar reactions, the angular correlations 
of gamma rays observed in coincidence with an axially symmetrical counter of finite but 
small size situated at zero degrees to the incident beam have no contribution from the y = 0 


tIssued as A.E.C.L. No. 1294. 

tTbe relevant Clebsch-Gordan coefficient is (cy. 1\|bB), where c is the spin of state in Ne, / is the angular 
momentum of the alpha particle, and 0 is the spin of the state in the compound nucleus Mg™. y, A, and @ are the 
projections of the angular momenta ¢, /, and 6 along the quantization axis which in this case is the accelerated 
beam of C!* ions. Because both the target nucleus and bombarding particle have zero spin B is zero. Consequently 
if y is zero then A must also be zero and the sum c+/+6 must be even otherwise the Clebsch-Gordan coefficient 
vanishes. Because both the alpha particle and C!? have even parity and parity is conserved in the reaction it is 
clear that if the state c has parity (—1)°*! then the Clebsch-Gordan coefficient is zero. 
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magnetic substate if the residual of state of spin c has parity (—1)**. The main contribution 
to the angular correlation will be from the y = +1 magnetic substates. The contribution 
from the y = +2 substates will usually be much smaller and be given by an equation similar 
to equation (29) of LF. 

(2) In the study of the angular distributions of alpha particles from such reactions as the 
C2(C2 a) Ne®* reaction no alpha particles will be observed at zero degrees to the beam if the 
state in the residual nucleus of spin c has parity (—1)*+!. Conversely if alpha particles are 
observed at zero degrees to the beam the parity of the state of spin c is (—1)*. It should be 
noted in this case that it is possible for no alpha particles to be emitted at zero degrees to 
the beam due to interference between states in the compound nucleus, This situation is unlikely 
and extremely unlikely to persist at different excitations in the compound nucleus. 

The effects on angular correlations and angular distributions of the above considerations 
have been observed recently. The C!*(C,ay)Ne* angular correlations from the 2— third 
excited state of Ne?® measured by Gove et al. (1961) show that the m = 0 substate of the 
2— excited state is not populated. The angular correlations can be fitted if only the m = +1 
substates are populated. 

The angular distributions of the alpha particles leading to the ground, first and second 
excited states of Ne®® from the C!(C!2,a)Ne?° reaction are in agreement with the known 
positive parity of these states (Almqvist e¢ al. 1961 and Clark et al. 1961). The third excited 
state of Ne? is now known to be a 2— state (Gove et al. 1961) and in this case the angular 
distributions measured as a function of C! energy by Almqvist et al. (1961) show a pronounced 
minimum in the differential cross section at zero degrees to the beam and at all bombarding 
energies. The small but non-zero yield at zero degrees is attributed by Almqvist et al. (1961) 
to the finite size of the counter detecting the alpha particles and to background. 

It is worth emphasizing in conclusion that measurements on heavy ion reactions such as 
those described above or from inelastic alpha scattering from even-even nuclei can decide 
the parity of states whose angular momenta have been determined. For example, the 6.44-Mev 
state in Mg** which has been assigned spin zero (Broude and Gove 1961) was observed in 
the C!2(O!5a)Mg* reaction (Hinds et al. 1961). In this case the observation of the alpha- 
particle group at any angle in the reaction is sufficient to assign positive parity to the state 
at 6.44 Mev in Mg” since the reaction cannot populate 0— levels. 


I would like to acknowledge my thanks to Dr. J. A. Kuehner for a stimulating discussion. 
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Magnetic Spectrometer Measurements of Levels of Ne** 


In connection with recent experiments (see accompanying letters) on the levels of Ne”, 
spectra of alpha particles from C"(C”,a~)Ne*® have been measured using a magnetic spectro- 
meter (Enge 1958). Figure 1 contains spectra measured at 0°, 10°, and 50° in the laboratory 
system and Table I summarizes the excitation energies in Ne*® obtained; the numbers in the 
left-hand column identify the corresponding peaks. The energy resolution obtained was 
AE/E ~ 1% and was limited by the target thickness. The spectrometer energy calibration 
was obtained assuming the excitation corresponding to group 2 to be 4.25 Mev. In general 
probable errors of 0.02 Mev are estimated; however, energy differences are given somewhat 
more accurately. In particular the difference in excitation energy corresponding to groups 
4 and 5 is obtained as 0.17+0.01 Mev. 

These spectra confirm the existence of a level at 5.81 Mev, first observed in the Na*3(p,a) Ne” 
reaction (Temmer 1961). Additional measurements have established spin 1— for the new 
5.80-Mev level as well as 3— for «hat known at 5.64 Mev (Kuehner 1961) and values of the 
ratio Ty/f of 0.003(+0.003) and 0.07(+0.01) respectively (Kuehner and Almqvist 1961). 


*Issued as A.E.C.L. No. 1283. 
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Fic. 1. Spectra of alpha particles from the reaction C!2(C!2,a) Ne*® measured with a magnetic spectrometer are 
presented for angles of 0°, 10°, and 50° in the laboratory system. Included in each spectrum is a scale giving the 
excitation energy in Ne®®, The numbers on the peaks identify them with entries in Table I. A region of the spectrum 
at 0° was obscured due to the presence at the detect >r of a portion of the transmitted beam. 
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TABLE I 


Summary of Ne? level energies in Mev obtained from spectrometer measurements of Fig. 1. The 
corresponding alpha- particle groups are identified by the numbers in the left-hand column. Values of 
the level spins and parities, where known, are given in the right-hand column 








Peak “Weighted Previous 
identification 0° 10° 50° average valuet Jn 
2 4.25* 4.25* 4.25 4+¢ 
3 4.97 4.97 4.97 2-25 
4 5.64 5.64 5.64 5.64 5.63 3—¢ 
5 5.77 5.80 5.83 5.81 1-¢ 
6 6.17 6.17 
7 6.72 6.75 6.74 6.75 0+4 
8 6.85 6.88 6.87 
9 7.05 7.05 7.05 (4—)¢ 
10 7.20 7.19 7.19 7.19 7.19 3-4 
11 7.25 7.25 7.23 0+¢ 
12 7.45 7.47 7.47 7.46 7.46 2+¢ 
13 7.67 7.63 7.65 (6+)¢ 
14 7.86 7.88 7.85 7.86 7.86 2+¢ 
15 7.96 8.00 7.92 7.93 
16 8.52 8.52 
17 8.92 8.91 8.92 


*Energy scale normalization point. 

tTaken from Ajzenberg and Lauritsen (1959). 
“Broude and Gove (1961). 

’Gove et al. (1961). 

©Kuehner (1961). 

4Cameron (1953). 

‘Litherland et al. (1961). 


Preliminary evidence for several new levels is presented. The group of alpha particles 
labelled 9 corresponds to a level at 7.05-Mev excitation and is identified with a level observed 
in earlier (a,y) coincidence studies as 7.02(+0.02)-Mev excitation (Litherland ef al. 1961). 
The absence of this group in the 0° spectrum suggests a spin-parity combination of 7 = 
(—1)/4*! for this level because, as was pointed out by Litherland (1961), only levels of Ne” 
with * = (—1)% can show yield at 0° in the C'!#(C2,a)Ne” reaction. In this connection 
additional spectra taken at 0° but not illustrated have shown ao, a1, and a2 to be present thus 
confirming these levels all to have + parity since their spins J = 0, 2, and 4 (Broude and 
Gove 1961) respectively. 

A shoulder, 13, on the low-energy side of group 12 corresponds to an excitation in Ne” 
of ~ 7.7 Mev and may be identified with a new level near this excitation suggested by y,7 
coincidence measurements (Litherland et al. 1961). The fact that this shoulder is strong at 0° 
requires the corresponding level in Ne” to have parity (—1)/. 

A possible additional new level in Ne” at 7.93-Mev excitation is suggested by the appear- 
ance of a group at all three angles. Its appearance in the 0° spectrum requires the parity to 
be (—1)/. 

Further work will be required on these levels as well as others suggested by these spectra 
to establish their identities. In particular, groups 6, 8, 16, and 17 may be from contaminant 
reactions. It is interesting to note that if a rotational band is associated with the 1 — (5.81 Mev) 
and 3— (7.19 Mev) levels (Litherland et al. 1961a) the spectra show no evidence for a 2— 
member near 6.3 Mev of comparable strength. The absence of 2— and 4— members of this 
band requires a K quantum number of 0 and implies an octupole shape deformation or vibration 
of the nucleus. 

Further measurements using this reaction and the Na*3(p,«) Ne” reaction are planned using 
the instrument in its broad range spectrograph mode of operation instead of the spectrometer 
mode used in this preliminary survey. 
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A 7.02-Mev Level in Ne™* 


A new level in Ne* at 7.02-Mev excitation energy has been observed in a study of coinci- 
dences'between alpha particles and gamma rays from the reaction C!*(C" ay)Ne*®. The new 
level was sought deliberately as a result of a prediction based upon the collective model. 
Recently it has been shown by Clark et a/. (1961) that the cascade electric quadruple transitions 
from the 4.25-Mev (4+) and the 1.63-Mev (2+) levels are both considerably enhanced over 
the single particle estimate. Also the levels at 4.97 and 5.63 Mev have been shown to have 
spins and parities of 2— and 3— respectively (Gove e¢ al. 1961 and Kuehner 1961). The 
presence of a 5— level at 8.84 Mev (McDermott et al. 1960) suggests that, as the collective 
model may be relevant for Ne”, there should be a level in Ne” near 7-Mev excitation with spin 
and parity 4—. Such a level cannot decay by alpha-particle emission to O" because of parity 
conservation and should consequently decay entirely by gamma-ray emission. 

A level at 7.02 Mev has now been located and it has some of the properties expected for a 
4— rotational level based on the 2— level in Ne® at 4.97 Mev. The level was located by 
observing the spectrum of alpha particles, in a 10,000 ohm cm silicon p-m solid state counter, 
in coincidence with the gamma rays from the low-lying levels of Ne*®. The alpha-particle 
counter was situated at 0° to the incident beam of C” ions which bombarded a 50-yg/cm? 
carbon foil. Three 0.4-mg/cm? aluminum foils were used to prevent the C” ions from hitting 
the alpha-particle counter. The gamma rays were observed with the aid of a 5 in. diameter 
by 6 in. long NalI(TI) crystal. 

Figure 1(a) shows the spectrum of alpha particles in coincidence with gamma rays in three 


C 12 se C--Ne?+ os 


e—e |:2 «10° INCH AL 


300 
x—x 6 x 10° INCH AL 
a « 
1-02 “4:25 
200 1-63 


100 


COUNTS PER CHANNEL 





30 40 50 60 70 80 90 
CHANNEL NUMBER 


is 1-63 GAMMA-RAY PULSE 


SPECTRUM 
ENERGIES IN Mev 


(b) 


COINCIDENCE COUNTS 





10 6=6p20— 3900s #0 00s 6—s FO 80 
CHANNEL NUMBER 


Fic. 1. (a) Portions of the spectrum of alpha particles from the reaction C!2(C",a)Ne™® are shown taken in 
coincidence with gamma rays. The effect of introducing an extra 0.4 X 10-4 inches of aluminum absorber is shown. 
The groups are labelled by the excitation energy in Mev of the level in the residual nucleus which they leave 


excited. 
(6) The gamma-ray pulse spectrum taken in coincidence with the alpha-particle group exciting the 7.02-Mev 


level in Ne® is shown. The energies of possible cascade gamma rays from the 7.02-Mev level are shown in Mev. 
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regions of the pulse spectrum. The shift of the peak positions upon insertion of an‘ additional 
0.4-mil aluminum foil in front of the counter proves that the groups are alpha-particle groups. 
The group designated a7,02 corresponds to the formation of Ne” in a new level at 7.02 Mev. 
Figure 1(b) shows the gamma-ray spectrum in coincidence with the alpha-particle group 
corresponding to a7,92 taken at 90° to the incident C” beam and Fig. 2 summarizes the results 
of the analysis of gamma-ray spectra taken at 30°, 60°, and 90°. The energies in Mev shown 
in Fig. 2 correspond to possible transitions from the new level at 7.02 Mev to lower levels in 
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. - ie 2. The energy level diagram of Ne® is shown together with the percentage branching of the new level at 
i ev. 


Ne. The absence of cascade gamma rays from the 5.63-Mev (3—) level, which is fed"by the 
1.39-Mev gamma ray of Fig. 1(b), is due to predominant alpha-particle emission (Gove et al. 
1961; Kuehner and Almqvist 1961). The most recent excitation energy of the new level is 
7.02+.02 Mev. In a subsidiary gamma-gamma coincidence experiment evidence has been 
obtained for another new level in Ne” at about 7.6 Mev which has some of the expected pro- 
perties for a 6+ level in Ne®. Because the evidence is not conclusive the level is shown by 
a dashed line in Fig. 2. 

At present there are no conclusive measurements of the spin and parity of the new level at 
7.02 Mev though the intensities of the gamma rays observed at 30°, 60°, and 90° are consistent 
with the assignment of 4—. Also recent measurements by Almqvist and Kuehner (1961) 
suggest a parity of (—1)/*!. However, the gamma-ray branching ratios shown in Fig. 2 
strongly support the assignment of 4— and the assignment of the level to a rotational band 
based on the 4.97-Mev level. The lack of evidence for an electric dipole transition from the 
new level to the 4.25-Mev level is expected because of the strong inhibition of electric dipole 
radiation in Ne® and the observed strong enhancement of interband electric quadrupole 
transitions (Clark ef al. 1961). It is hoped to measure the spin directly in the near future by 
alpha-gamma correlations (Gove et al. 1961a) and to determine the lifetime of the new state 
by the Doppler shift attenuation method (Devons et al. 1954; Litherland et al. 1961). 
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ANNOUNCEMENTS 





Frequency Measurement of Standard Frequency Transmissions'.? 


Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
as reference standard (with an assumed frequency of 9 192 631 770 c.p.s.). Frequency deviations 
from nominal are quoted in parts per 10!%. A negative sign indicates that the frequency is below 
nominal. 





GBR, 16 ke/s 

May MSF, 8-hour 24-hour WWVB, 
1961 60 ke/s average* averaget 60 ke/s 
—151 —151 —140 —148 
— 169 — 156 — 150 — 146 
— 159 —154 — 153 — 149 
— 153 — 155 —150 —151 
—149 — 145 —148 — 139 
—141 — 137 —138 —139 
— 139 — 135 — 137 — 144 
— 152 —147 — 136 —142 
—145 —146 —134 —149 
—152 —149 —147 —151 
—149 — 150 —147 —154 
— 159 — 156 —149 —146 
— 155 —155 —152 — 156 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. —147 
— 351 —154 N.M. —152 
—149 —152 — 154 —149 
N.M. —151 —150 N.M. 
N.M. — 155 —155 N.M. 
—151 — 153 —154 —154 
—153 —150 —154 —149 
— 153 —152 —150 —149 
—155 —153 — 150 — 150 
N.M. — 153 —152 N.M. 
N.M. —154 —152 —151 
N.M. N.M. N.M. — 153 
— 157 — 154 —154 —151 
N.M. —155 —153 —148 


18 
19 


& 


mew hy 0% 


am wt = 


mt tw 0 
SOND 


Ce G bt 
- oOo 


Average — 152 —151 — 148 —149 
Midmonthly 
mean of WWV — 148 


Note: N.M., no measurement. 
*Time of observations: 10:30 to 18:30 U.T. 
tTime of observations: 15 to 15 hours U.T. 





These measurements will no longer appear regularly in the Canadian Journal of Physics. 
However, individuals or organizations interested may obtain this information on a monthly 
basis by requesting it from the Director, Applied Physics Division, National Research Council, 
Ottawa, Canada. 


RECEIVED JUNE 9, 1961. S. N. KALRA 
DIVISION OF APPLIED PHysics, 

NATIONAL RESEARCH COUNCIL, 

Ottawa, CANADA, 


1[ssued as N.R.C. No. 6395. 
*Cf. Kalra, S. N. 1959. Can. J. Phys. 37, 1328. 
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MANUSCRIPTS 


General.—Manuscripts, in English or French, should be typewritten, double spaced, on 
paper 8} X11 in. The original and one copy are to be submitted. Tables and captions 
for the figures should be placed at the end of the manuscript. Every sheet of the manuscript 
should be numbered. Style, arrangement, spelling, and abbreviations should conform to the 
usage of recent numbers of this journal. Greek letters or unusual signs should be written plainly 
or explained by marginal notes. Characters to be set in boldface type should be indicated 
by a wavy line below each character. Superscripts and subscripts must be legible and carefully 
placed. Manuscripts and illustrations should be carefully checked before they are submitted. 
Authors will be charged for unnecessary deviations from the usual format and for changes 
made in the proof that are considered excessive or unnecessary. 

Abstract.—An abstract of not more than about 200 words, indicating the scope of the work 
and the principal findings, is required, except in Notes. 

References.—References should be listed alphabetically by authors’ names, un- 
numbered, and typed after the text. The form of the citations should be that used in current 
issues of this journal; in references to papers in periodicals, titles should not be given and 
only initial page numbers are required. The names of periodicals should be abbreviated in 
the form given in the most recent List of Periodicals Abstracted by Chemical Abstracts. All 
citations should be checked with the original articles and each one referred to in the text by 
the authors’ names and the year. 

Tables.—Tables should be numbered in roman numerals and each table referred to in the 
text. Titles should always be given but should be brief; column headings should be brief and 
descriptive matter in the tables confined to a minimum. Vertical rules should not be used. 
Numerous small tables should be avoided. 


ILLUSTRATIONS 


General.—All figures (including each figure of the plates) should be numbered con- 
secutively from 1 up, in arabic numerals, and each figure referred to in the text. The author’s 
name, title of the paper, and figure number should be written in the lower left corner of the 
sheets on which the illustrations appear. Captions should not be written on the illustrations. 

Line drawings.— Drawings should be carefully made with India ink on white drawing 
paper, blue tracing linen, or co-ordinate paper ruled in blue only; any co-ordinate lines that 
are to appear in the reproduction should be ruled in black ink. Paper ruled in green, yellow, 
or red should not be used. All lines must be of sufficient thickness to reproduce well. Decimal 
points, periods, and stippled dots must be solid black circles large enough to be reduced if 
necessary. Letters and numerals should be neatly made, preferably with a stencil (do NOT 
use typewriting) and be of such size that the smallest lettering will be not less than 1 mm 
high when the figure is reduced to a suitable size. Many drawings are made too large; originals 
should not be more than 2 or 3 times the size of the desired reproduction. Whenever possible 
two or more drawings should be grouped to reduce the number of cuts required. In such groups 
of drawings, or in large drawings, full use of the space available should be made; the ratio of 
height to width should conform to that of a journal page (43 X7}in.), but allowance must 
be made for the captions. The original drawings and one set of clear copies (e.g. small 
photographs) are to be submitted. 

Photographs.—Prints should be made on glossy paper, with strong contrasts. They 
should be trimmed so that essential features only are shown and mounted carefully, with 
rubber cement, on white cardboard, with no space between those arranged in groups. In 
mounting, full use of the space available should be made. Photographs are to be sub- 
mitted in duplicate; if they are to be reproduced in groups one set should be mounted, 
the duplicate set unmounted. 


REPRINTS 


A total of 100 reprints of each paper, without covers, are supplied free. Additional reprints, 
with or without covers, may be purchased at the time of publication. 

Charges for reprints are based on the number of printed pages, which may be calculated 
approximately by multiplying by 0.6 the number of manuscript pages (double-spaced type- 
written sheets, 8}X11lin.) and including the space occupied by illustrations. Prices and 
instructions for ordering reprints are sent out with the galley proof. 
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